抽樣方案設(shè)計-農(nóng)戶糧食儲量抽樣方案

PAGEPAGE50第二部分農(nóng)戶儲糧損失調(diào)查抽樣方案一、農(nóng)戶儲糧損失調(diào)查抽樣方法1.本項調(diào)查所面臨的客觀條件所研究的對象總體特別巨大（包括全國28個省/自治區(qū)/直轄市）區(qū)的近3億農(nóng)戶），地域特別寬廣，農(nóng)戶儲糧方式和儲糧數(shù)量存在很大差異。序號農(nóng)戶儲糧損失調(diào)查——調(diào)查區(qū)域調(diào)查實施單位技術(shù)支持單位1江蘇、浙江、安徽、福建、江西、湖北、湖南、廣東、廣西、海南、重慶、四川、貴州、云南等14省（自治區(qū)、直轄市）各省（自治區(qū)、直轄市）糧食局成都儲藏所2河北、山西、山東、河南、西藏、陜西、甘肅、青海、寧夏、新疆等10省（自治區(qū)）各?。ㄗ灾螀^(qū)）糧食局河南工大3內(nèi)蒙古、遼寧、吉林、黑龍江等4省（自治區(qū)）各?。ㄗ灾螀^(qū)）糧食局遼寧糧科所2.制定本次調(diào)查抽樣方案的指導思想充分考慮全國各區(qū)域糧食儲藏量的數(shù)量差異，盡可能兼顧樣本調(diào)查戶在各?。ㄗ灾螀^(qū)、直轄市）、各樣本縣（區(qū)、市）中的分散性和均勻分布，從而保證樣本對總體的代表性。3.確定的抽樣方法根據(jù)統(tǒng)計抽樣方法的設(shè)計原理，為確保抽樣調(diào)查方法的科學性和可操作性，通過對各種抽樣方法的比較研究，并結(jié)合農(nóng)戶儲糧實際，確定本調(diào)查采用的抽樣方法為一種優(yōu)化的組合抽樣方式：分層三階段不等概率PPS系統(tǒng)抽樣法及線性系統(tǒng)隨機抽樣法。即：第一階段——采用分層PPS線性系統(tǒng)抽樣法，以七大儲糧生態(tài)區(qū)域（具體劃分見附件2-1）為子總體（層），抽取樣本縣（區(qū)、市）；第二階段——采用線性系統(tǒng)隨機抽樣法，在樣本縣（區(qū)、市）中抽取樣本村；第三階段——采用現(xiàn)場隨機抽樣法，在樣本村中抽取樣本農(nóng)戶；同時——采用多階段抽樣的最優(yōu)分配法對樣本分布影響最大的第一階段的樣本進行最優(yōu)化分配。二、農(nóng)戶儲糧損失調(diào)查抽樣方案本次抽樣調(diào)查中所采用的抽樣技術(shù)，包括三階段抽樣、分層抽樣、不等概率PPS抽樣、線性系統(tǒng)隨機抽樣及現(xiàn)場隨機抽樣法，是我國在大規(guī)模抽樣調(diào)查和市場調(diào)查中常用的抽樣方法。在本次農(nóng)戶儲糧損失抽樣調(diào)查中，根據(jù)調(diào)查對象的性質(zhì)和本次調(diào)查的實際情況，綜合采用上述抽樣方法設(shè)計了一種優(yōu)化的組合抽樣技術(shù)方案如下：第一階段——以七大儲糧生態(tài)區(qū)域為子總體，抽取樣本縣（區(qū)、市）。本階段采用分層PPS線性系統(tǒng)抽樣法（抽樣過程中所涉及到的相關(guān)計算公式及計算過程詳見附件2-2）。1.樣本容量的確定及要達到的精度指標由抽樣技術(shù)理論可知，在本子課題將要進行的大規(guī)模農(nóng)戶儲糧損失調(diào)查中，抽樣數(shù)目主要取決于以下因素：（1）被調(diào)查現(xiàn)象數(shù)量（農(nóng)戶儲糧損失數(shù)量）的波動幅度，即標準差或方差的大??；（2）允許誤差或允許相對誤差率的大??；（3）置信概率或概率把握程度的高低；（4）抽樣技術(shù)方式和方法的不同；（5）調(diào)查對象無回答率的高低（不負責任、有意錯答等非隨機抽樣誤差問題不在此研究）。其中：抽樣極限誤差與抽樣概率度和抽樣平均誤差緊密相關(guān)。抽樣極限誤差與抽樣概率度的值可根據(jù)抽樣調(diào)查的實際條件來設(shè)定。在一定的概率保證程度下（即抽樣概率度一定），要想提高抽樣估計的精度，就必須縮小抽樣極限誤差，就必須通過選擇最佳的抽樣方式、方法和估計量來降低抽樣平均誤差。抽樣極限誤差可采用抽樣相對允許誤差來表示。在其他因素確定的情況下，允許誤差越小，所需樣本容量就越大，需花費的調(diào)查成本就越大。在本次農(nóng)戶儲糧損失抽樣調(diào)查中，按照抽樣調(diào)查的一般規(guī)則，確定抽樣的相對（極限）誤差率為3％；置信概率為95％（對應(yīng)正態(tài)分布的概率度為1.96）。在抽樣數(shù)目的計算過程中，除需要給定允許誤差或允許相對誤差率的數(shù)據(jù)之外，很重要的就是要取得關(guān)于被研究對象的均值、波動幅度、標準差或方差及單位總數(shù)等的數(shù)據(jù)資料，這些數(shù)據(jù)資料信息通常來自現(xiàn)有歷史資料和預(yù)調(diào)查數(shù)據(jù)，而目前所收集的歷史資料多數(shù)屬于典型調(diào)查，數(shù)據(jù)的代表性較差，無法用于推斷總體，因此，在本方案中樣本容量主要取決于試抽樣調(diào)查數(shù)據(jù)波動幅度的大小。為了確定我國農(nóng)戶儲糧損失數(shù)據(jù)的波動幅度，2007年7月份以來，在全國范圍內(nèi)（涵蓋7大儲糧生態(tài)區(qū)域8省/自治區(qū)/直轄市）隨機抽取32個樣本縣（區(qū)、市）對200個農(nóng)戶進行了預(yù)抽樣調(diào)查，取得了191個農(nóng)戶的有效調(diào)查數(shù)據(jù)。本次預(yù)調(diào)查采用按統(tǒng)一調(diào)查表格，深入農(nóng)村一家一戶，采用問卷調(diào)查的方式進行，調(diào)查問卷由調(diào)查人員進行填寫。綜合考慮樣本分布的均勻性和樣本對被研究現(xiàn)象的代表性這兩大因素，同時為了達到抽樣誤差最小的原則，分別對各儲糧生態(tài)區(qū)域所包含的縣級單位數(shù)和各儲糧生態(tài)區(qū)域2006年的糧食產(chǎn)量（通過對預(yù)調(diào)查樣本損失數(shù)據(jù)的分析測算，糧食產(chǎn)量與糧食損失量有極大的正相關(guān)關(guān)系）測算了與其對應(yīng)的尼曼最佳抽樣數(shù)目，并按杰森多指標最優(yōu)平均法進行了樣本在各區(qū)域的最優(yōu)分配。在“糧食儲藏量（極大相關(guān)指標為糧食生產(chǎn)量）越大的地區(qū)，越應(yīng)該成為糧食儲藏損失調(diào)查的重點”的原則下，以各縣（區(qū)、市）2006年的糧食產(chǎn)量為輔助標志（權(quán)數(shù)）進行不等概率（PPS）系統(tǒng)抽樣，根據(jù)相關(guān)因數(shù)計算出應(yīng)該在全國抽取不低于361個縣（區(qū)、市）和1688個農(nóng)戶(在多階段抽樣條件下，按3倍設(shè)計效應(yīng)系數(shù)計算，抽取的農(nóng)戶數(shù)應(yīng)為5064戶，按379個樣本縣（區(qū)、市）計算，平均每個縣（區(qū)、市）應(yīng)抽取不低于13.4戶)，而在實際抽樣過程中，考慮到小數(shù)進位和樣本點重復的影響，通常將樣本擴大5％計算，從而確定初始樣本縣（區(qū)、市）數(shù)為379個縣（區(qū)、市），最終確定的有效抽樣縣（區(qū)、市）數(shù)目為373個縣（區(qū)、市）。2.第一階段樣本的最優(yōu)分配由統(tǒng)計抽樣理論可知：在多階段抽樣中，第一階段抽樣的精度對整個多階段抽樣有極大的影響（實證分析表明第一階段抽樣誤差通常占到整個抽樣誤差的90％~95％）。因此，在本次抽樣方案設(shè)計中，首先對第一階段的樣本數(shù)進行了最優(yōu)化的樣本分配，所采用的方法是尼曼（Neyman）最優(yōu)分配法和多指標綜合的杰森（Jessen）最優(yōu)平均法。（尼曼最優(yōu)分配法是在總體各層中分配樣本容量時，構(gòu)造一個關(guān)于抽樣方差的拉格朗日極值函數(shù)，在一定的約束條件下，使抽樣方差最小的樣本容量分配；杰森最優(yōu)平均法是在處理不同度量的多指標綜合時，對各指標計算出的尼曼最優(yōu)分配數(shù)進行折衷處理的方法）。優(yōu)化處理后具體樣本分配情況如下表：多指標杰森最優(yōu)平均法——第一階段樣本分配情況區(qū)域累計產(chǎn)量總縣數(shù)損失率標準差產(chǎn)量*標準差產(chǎn)量樣本分配比縣*標準差縣樣本分配比應(yīng)抽樣本縣數(shù)最終樣本縣數(shù)118480421191.0962025454.0320.0044130.42400.068413.7986142363226182380.395214354698.630.031494.05760.049315.2913163647692101131.24280443358.820.1758140.34600.073647.26254741723370836230.9621165805507.60.3624599.38830.3143128.229812951419803365931.0252145558240.50.3181607.94360.3188120.69421216217132641670.704115288309.180.0334117.58470.061718.0158197274315321751.24234069962.740.0745217.35000.114035.707836合計或平均46640208520281.242457545531.411907.09421379382具體應(yīng)用方法是：（1）將子總體的全部農(nóng)業(yè)縣級單位按隨機的順序（或無關(guān)標志排隊方式）排列；（2）將各縣（區(qū)、市）的糧食產(chǎn)量數(shù)據(jù)對應(yīng)地列出，并累計起來，形成糧食產(chǎn)量累計數(shù)列；（3）計算得出抽樣距離數(shù)據(jù)k，并得到隨機起點r；（4）通過k與r結(jié)合，抽出全部調(diào)查縣（區(qū)、市）。在確定了各生態(tài)區(qū)域樣本縣（區(qū)、市）數(shù)目的基礎(chǔ)上，對各生態(tài)區(qū)域分別按照無關(guān)標志排隊PPS系統(tǒng)抽樣法抽出了第一階段的樣本單位即樣本調(diào)查縣（區(qū)、市），如下表所示：所屬儲糧生態(tài)區(qū)域28省（自治區(qū)、直轄市）樣本縣糧食產(chǎn)量（噸）1區(qū)：包括西藏、青海、新疆、甘肅和四川等5省/自治區(qū)的119個縣（級行政單位）。共抽14個縣（級行政單位）1西藏自治區(qū)堆龍德慶縣239321西藏自治區(qū)丁青縣232561西藏自治區(qū)扎囊縣214581西藏自治區(qū)日喀則市654921西藏自治區(qū)江孜縣522721西藏自治區(qū)白朗縣392451西藏自治區(qū)波密縣169451青海省尖扎縣149181青海省德令哈市172641甘肅省臨夏縣1078981甘肅省和政縣433971甘肅省迭部縣107771四川省壤塘縣31281四川省巴塘110982區(qū)：包括青海、新疆、甘肅、寧夏、內(nèi)蒙、黑龍江、吉林、河北、山西和陜西等10省/自治區(qū)的238個縣（級行政單位）。共抽16個縣（級行政單位）2青海省門源回族自治縣394902新疆維吾爾自治區(qū)新和縣1082002新疆維吾爾自治區(qū)伽師2057872新疆維吾爾自治區(qū)昭蘇縣453012甘肅省涼州區(qū)6097462寧夏回族自治區(qū)靈武市1777542內(nèi)蒙古土默特右旗6123352內(nèi)蒙古寧城縣5204732內(nèi)蒙古科爾沁左翼后旗7375002內(nèi)蒙古達拉特旗5752392內(nèi)蒙古五原縣3700512內(nèi)蒙古科爾沁右翼前旗5587442黑龍江省泰來4315392吉林省前郭爾羅斯蒙古族自治縣16500002吉林省鎮(zhèn)賚縣5038582河北省萬全縣1209003區(qū)：包括內(nèi)蒙、黑龍江、吉林和遼寧等4省/自治區(qū)的113個縣（級行政單位）。共抽40個縣（級行政單位）3內(nèi)蒙古阿榮旗7505003內(nèi)蒙古扎蘭屯5015003黑龍江省依蘭6473073黑龍江省賓縣6814093黑龍江省木蘭3696323黑龍江省阿城5200313黑龍江省雙城15953383黑龍江省五常17056763黑龍江省龍江11100973黑龍江省依安5821303黑龍江省克山5182873黑龍江省訥河11474253黑龍江省虎林市4308503黑龍江省集賢縣5106473黑龍江省肇州縣6505613黑龍江省嘉蔭縣1242573黑龍江省湯原縣4764273黑龍江省富錦市9501263黑龍江省海林市2209503黑龍江省嫩江縣7426463黑龍江省望奎縣5601903黑龍江省慶安縣6488433黑龍江省肇東市16125653吉林省農(nóng)安縣22905193吉林省九臺市9466573吉林省榆樹市23468003吉林省德惠市15135583吉林省永吉縣3940743吉林省舒蘭市72282733吉林省磐石市5477733吉林省梨樹縣20400453吉林省伊通滿族自治縣7605003吉林省公主嶺市22200523吉林省東豐縣6328643吉林省輝南縣3609013吉林省扶余縣17702283遼寧省康平縣5294913遼寧省寬甸滿族自治縣1101153遼寧省西豐縣3209733遼寧省昌圖縣15011354區(qū)：包括甘肅、寧夏、遼寧、河北、山西、山東、河南、陜西、江蘇、安徽和湖北等11省的623個縣（級行政單位）。共抽128個縣（級行政單位）4甘肅省甘谷縣1368484甘肅省環(huán)縣1018424甘肅省臨洮縣1838354寧夏回族自治區(qū)同心縣1065594遼寧省新民市8106854遼寧省普蘭店市3437754遼寧省海城市5081484遼寧省黑山縣7665844遼寧省蓋州市1579564遼寧省阜新蒙古族自治縣11255604遼寧省大洼縣5311254遼寧省朝陽縣4658974遼寧省北票市3917604河北省正定縣3086204河北省無極縣3498684河北省辛集市5128994河北省新樂市2941214河北省灤南縣3557714河北省遷安市2177144河北省臨漳縣4873674河北省磁縣2831294河北省館陶縣2296584河北省內(nèi)丘縣1484354河北省寧晉縣6038994河北省臨西縣2233544河北省徐水縣3544804河北省安新縣1903834河北省涿州市2933284河北省平泉縣2562414河北省鹽山縣2544924河北省泊頭市3406714河北省永清縣1484894河北省棗強縣2997434河北省景縣5274244河北省深州市4501114山西省長子縣2222234山西省山陰縣1029564山西省介休市756464山西省河津市1133174山西省洪洞縣2557824山西省汾陽市1657114山東省商河縣6091564山東省即墨市5301614山東省平度市13142064山東省桓臺縣3714394山東省沂源縣463664山東省萊州市4725194山東省臨朐縣2978684山東省諸城市6845794山東省高密市6487794山東省魚臺縣2847854山東省梁山縣3404974山東省寧陽縣5051384山東省新泰市4209474山東省五蓮縣2000074山東省郯城縣6191594山東省費縣3397964山東省臨沭縣2920724山東省臨邑縣5786334山東省夏津縣1920824山東省禹城市6710094山東省莘縣6214274山東省冠縣5234164山東省陽信縣3216174山東省鄒平縣5583324山東省成武縣3999804山東省鄆城縣5378234河南省中牟縣3341984河南省杞縣4792474河南省開封縣4492624河南省汝陽縣1509504河南省葉縣4582874河南省安陽縣3761514河南省滑縣10859064河南省浚縣6054084河南省延津縣3384804河南省衛(wèi)輝市2991904河南省武陟縣4556604河南省清豐縣4842744河南省濮陽縣6947654河南省襄城縣5128574河南省郾城區(qū)2516084河南省靈寶市1840474河南省內(nèi)鄉(xiāng)縣2064104河南省唐河縣9426704河南省民權(quán)縣4897004河南省柘城縣5000144河南省夏邑縣7213354河南省永城市9286004河南省扶溝縣3895364河南省西華縣4365004河南省商水縣7393064河南省淮陽縣6605004河南省鹿邑縣6053304河南省西平縣7118704河南省上蔡縣8031654河南省確山縣3545594河南省遂平縣4500294陜西省藍田縣2913054陜西省鳳祥縣2478834陜西省三原縣2306614陜西省武功縣2185924陜西省白水縣797124陜西省宜川縣361754陜西省清澗縣440364陜西省鎮(zhèn)安縣1111834江蘇省銅山縣5330284江蘇省邳州市6065384江蘇省東?？h9345284江蘇省漣水縣6874564江蘇省盱眙縣6592114江蘇省響水縣3336454江蘇省射陽縣7680934江蘇省東臺市6819314江蘇省寶應(yīng)縣6175934江蘇省興化市11698204江蘇省宿豫區(qū)4253244江蘇省泗陽縣4416254安徽省懷遠縣9992564安徽省鳳臺縣493874安徽省鳳陽縣5397004安徽省太和縣6963094安徽省潁上縣7864854安徽省蕭縣6568434安徽省泗縣6791264安徽省渦陽縣9481074安徽省蒙城縣8634504湖北省鄖西縣1649125區(qū)：包括甘肅、河南、陜西、江蘇、安徽、浙江、福建、江西、湖北、湖南、重慶、四川、貴州、廣東和廣西等15省/自治區(qū)/直轄市的593個縣（級行政單位）。共抽120個縣（級行政單位）5河南省光山縣4115885河南省固始縣8854375江蘇省高淳縣1765445江蘇省溧陽市3661415江蘇省太倉市1894835江蘇省如東縣7644645江蘇省通州市4981185江蘇省江都市4499165江蘇省靖江市2733885江蘇省姜堰市5057345安徽省肥西縣4794715安徽省當涂縣3559795安徽省潛山縣2344465安徽省桐城市3815195安徽省定遠縣7836045安徽省天長市6130255安徽省無為縣4054105安徽省霍邱縣10276345安徽省東至縣2345045浙江省蕭山區(qū)2682305浙江省余姚市1954645浙江省嘉善縣1719055浙江省安吉縣1449775浙江省浦江縣846805浙江省江山市2349415浙江省慶元縣571525福建省仙游縣1899765福建省建寧縣1100925福建省長泰縣715055福建省武夷山市1456135福建省武平縣1846185江西省南昌縣7818005江西省安義縣1080905江西省蘆溪縣1124905江西省都昌縣3221855江西省貴溪市3076355江西省寧都縣3146705江西省南康市2359955江西省新干縣2764205江西省萬安縣2044165江西省上高縣2560155江西省豐城市7834775江西省高安市6119565江西省宜黃縣1475885江西省玉山縣1847425江西省鄱陽縣6957645湖北省陽新縣3223185湖北省秭歸縣1159735湖北省襄陽區(qū)8597685湖北省谷城縣2018065湖北省宜城市4628975湖北省沙洋縣5837895湖北省大悟縣2677215湖北省漢川市3843375湖北省監(jiān)利縣8595585湖北省洪湖市5034275湖北省紅安縣2663375湖北省蘄春縣3591155湖北省黃梅縣3613805湖北省廣水市4195425湖北省巴東縣1931865湖北省仙桃市6571335湖南省長沙縣6094485湖南省寧鄉(xiāng)縣8867665湖南省株洲縣3281355湖南省炎陵縣934935湖南省湘潭縣8778995湖南省衡陽縣6662635湖南省衡山縣1979785湖南省耒陽市5213925湖南省新邵縣3025125湖南省洞口縣4947215湖南省武岡市4637895湖南省華容縣4988265湖南省汨羅市4652475湖南省漢壽縣5026125湖南省桃源縣6173395湖南省桑植縣1415625湖南省沅江市3614395湖南省嘉禾縣1546335湖南省祁陽縣6018595湖南省江永縣1458155湖南省新晃自治縣843075湖南省漣源市4241295湖南省花垣縣707925重慶市綦江縣4024075重慶市大足縣4733835重慶市豐都縣3357855重慶市忠縣4323735重慶市奉節(jié)縣4355375重慶市酉陽土家族苗族自治縣3531605重慶市江津市6900435重慶市永川市5033635四川省雙流2804975四川省彭州2637285四川省榮縣4278835四川省合江縣4472625四川省廣漢市2979245四川省三臺縣6902155四川省江油市3050275四川省蒼溪縣3881215四川省威遠縣3308075四川省犍為縣1950025四川省南部縣4981455四川省儀隴縣5068695四川省仁壽縣7908645四川省珙縣1430335四川省鄰水縣4304675四川省宣漢縣5419925四川省大竹縣5584165四川省漢源縣1555145四川省平昌縣3714045四川省樂至縣3768005貴州省正安縣2576315貴州省德江縣1602155貴州省榕江縣977295廣東省陽山縣1405695廣西壯族自治區(qū)陽朔縣1133305廣西壯族自治區(qū)永?？h1466735廣西壯族自治區(qū)宜州市2668986區(qū)：包括四川、云南、貴州和廣西等4省/自治區(qū)的167個縣（級行政單位）。共抽19個縣（級行政單位）6四川省西昌市2508266四川省昭覺941546云南省尋甸回族自治縣1877306云南省會澤縣2768486云南省峨山縣425106云南省大關(guān)縣714266云南省景東縣1287146云南省元謀縣685076云南省祥云縣1529806云南省蘭坪縣733966貴州省盤縣3283786貴州省綏陽縣2722636貴州省仁懷市2562356貴州省興義市2520876貴州省大方縣3182436貴州省織金縣3276856貴州省凱里市972726貴州省獨山縣1117666廣西壯族自治區(qū)鳳山縣376637區(qū)：包括福建、云南、廣東、廣西和海南等5省/自治區(qū)的175個縣（級行政單位）。共抽36個縣（級行政單位）7福建省平和縣1419447云南省瀾滄縣1606047云南省蒙自縣1207557云南省文山縣1273377云南省景洪市1300747廣東省增城市1771407廣東省三水區(qū)406437廣東省鶴山市953037廣東省廉江市4523717廣東省吳川市2095587廣東省化州市3799107廣東省廣寧縣1798937廣東省德慶縣1206697廣東省惠東縣2463977廣東省豐順縣1439927廣東省海豐縣1744447廣東省龍川縣3174257廣東省陽西縣1614237廣東省清新縣1871327廣東省饒平縣1972847廣東省惠來縣2224717廣東省云安縣1360187廣西壯族自治區(qū)隆安縣1310327廣西壯族自治區(qū)橫縣3280377廣西壯族自治區(qū)岑溪市2753107廣西壯族自治區(qū)靈山縣4588807廣西壯族自治區(qū)平南縣3149957廣西壯族自治區(qū)容縣2533567廣西壯族自治區(qū)博白縣4792237廣西壯族自治區(qū)北流市4249727廣西壯族自治區(qū)靖西縣2107897廣西壯族自治區(qū)都安瑤族自治縣1281287廣西壯族自治區(qū)象州縣1784147海南省天等縣1407467海南省文昌市1362947海南省樂東黎族自治縣135387第二階段——在樣本縣（區(qū)、市）中抽取樣本村。第二階段（縣/區(qū)/市抽村）的抽樣方法：采用線性系統(tǒng)隨機抽樣法（也稱無關(guān)標志排隊等距抽樣法）。具體方法分為兩種情況：第一種情況：具備資料可以建立全縣（區(qū)、市）各村的名錄框架（全部村名的一覽表）。在此條件下，第一步，先對全縣（區(qū)、市）各村從1開始依次編碼；第二步，用最大編碼數(shù)除以4（樣本村數(shù)），其結(jié)果稱為抽樣距離，記為k；第三步，在數(shù)碼1—k中作純隨機抽樣，抽出的隨機數(shù)記為r,則編碼為r,r+k,r+2k和r+3k的村即為被抽中的樣本村。第二種情況：不具備資料無法建立全縣（區(qū)、市）各村的名錄框架，只能建立全縣（區(qū)、市）各鄉(xiāng)鎮(zhèn)的名錄框架。在此條件下，就分兩步進行縣（區(qū)、市）抽村的抽樣。第一步：縣（區(qū)、市）抽鄉(xiāng)鎮(zhèn)（城關(guān)鎮(zhèn)除外，下同）。首先，對全縣（區(qū)、市）各鄉(xiāng)鎮(zhèn)從1開始依次編碼；第二，用最大編碼數(shù)除以4（樣本鄉(xiāng)鎮(zhèn)數(shù)），其結(jié)果稱為抽樣距離，記為k；第三，在數(shù)碼1—k中作純隨機抽樣，抽出的隨機數(shù)記為r,則編碼為r,r+k,r+2k和r+3k的鄉(xiāng)鎮(zhèn)即為被抽中的樣本鄉(xiāng)鎮(zhèn)。第二步，分別對每一個被抽中的鄉(xiāng)鎮(zhèn)，建立其所轄各村的名錄框架（全部村名的一覽表）。設(shè)共有k個村，在數(shù)碼1—k中作純隨機抽樣，抽出的隨機數(shù)記為r,則編碼為r的村即為被抽中的樣本村，則總共抽出4個樣本村。第二階段，全國373個樣本縣（區(qū)、市）共抽取1492個樣本村。第三階段——在樣本村中抽取樣本農(nóng)戶。本階段抽樣由調(diào)查人員通過調(diào)查現(xiàn)場隨機抽樣來實施入戶調(diào)查。按照抽樣設(shè)計，第三階段抽樣即村抽農(nóng)戶的數(shù)量為每村隨機抽取4戶農(nóng)戶，全國共抽取樣本農(nóng)戶5968戶進行農(nóng)戶儲糧損失入戶調(diào)查。根據(jù)我國農(nóng)戶儲糧的特點和本次抽樣調(diào)查的具體情況，在進行現(xiàn)場隨機抽樣時，設(shè)計了2種實施方法：第一法——有門牌號的農(nóng)戶（優(yōu)先采用）調(diào)查人員到達第二階段抽樣所抽中的樣本村以后，確認最鄰近街、巷或組的位置，然后從事先已編號的調(diào)查表中隨機抽取1份，則門牌號與所抽出的調(diào)查表編號后兩位數(shù)字相同的農(nóng)戶成為現(xiàn)場入戶調(diào)查的對象；若該組無此門牌號，則可用余數(shù)法——即：用100減去調(diào)查表后兩位數(shù)字的余數(shù)作為所抽中的門牌號；若該農(nóng)戶家中無人可重抽一次。余下的三個調(diào)查戶即為下一個最鄰近街、巷或組的同一個門牌號的農(nóng)戶。第二法——無門牌號的農(nóng)戶調(diào)查人員到達第二階段抽樣所抽中的樣本村以后，從事先已編號的調(diào)查表中隨機抽取1份，以村委會為中心，按右手原則，從村委會右邊第一戶農(nóng)戶起，抽取戶數(shù)與所抽中調(diào)查表編號后兩位數(shù)字相同的農(nóng)戶作為現(xiàn)場入戶調(diào)查的對象。若該村無此農(nóng)戶，則可用余數(shù)法——即：用100減去調(diào)查表后兩位數(shù)字的余數(shù)作為所抽中的農(nóng)戶；若該農(nóng)戶家中無人可重抽一次。余下的三個調(diào)查戶即為此農(nóng)戶右手方向每間隔20戶的農(nóng)戶。注：重抽方法——如果該戶農(nóng)戶沒有適合人選，則選擇向右最鄰近的另一戶農(nóng)戶進行調(diào)查。三、附件附件2-1——7大儲糧生態(tài)區(qū)域的劃分附件2-2——抽樣過程中所涉及到的相關(guān)計算公式及計算過程附件2-3——農(nóng)戶儲糧損失調(diào)查各總體參數(shù)的估計方法附件2-4——各?。ㄗ灾螀^(qū)、直轄市）調(diào)查工作任務(wù)的分配附件2-1——7大儲糧生態(tài)區(qū)域的劃分中國儲糧生態(tài)區(qū)域劃分圖1中國儲糧生態(tài)區(qū)域劃分區(qū)域代碼區(qū)域名稱主要分布范圍1高寒干燥儲糧區(qū)東起橫斷山區(qū)，西抵喀喇昆侖，南至喜馬拉雅，北達阿爾金山—祁連山北麓。包括西藏全部、青海南部及川、滇、新、甘的一部分。海拔4000m2低溫干燥儲糧區(qū)位于我國北部與西北部，包括內(nèi)蒙古的大部分，寧、陜、冀、甘的一部分及新疆的全部。以國界線為西界和北界，東界從大興安嶺的根河河口開始，沿大興安嶺西麓，向南延伸至阿爾山附近，然后向東沿洮兒河谷地跨越大興安嶺至烏蘭浩特以東，再沿大興安嶺東麓南下，經(jīng)突泉、扎魯特旗、開魯至奈曼，南界昆侖山、祁連山北麓至烏鞘嶺，長城、張北、沽源、圍場、阜新。3低溫高濕儲糧區(qū)位于我國的東北部，北、東、東南三界為國界，是我國位置最北、緯度最高的一個區(qū)。本區(qū)西界就是蒙新區(qū)的東界，南界自奈曼至彰武、康平、昌圖、鐵嶺、撫順、寬甸至鴨綠江邊，相當于活動積溫3200℃4中溫干燥儲糧區(qū)西鄰青藏高原，東瀕黃海和渤海，北面與東北地區(qū)及蒙新地區(qū)相接，以秦嶺北麓、伏牛山、淮河為界與華中區(qū)相接，此界相當于活動積溫4500℃或1月份均溫0℃、年降水量5中溫高濕儲糧區(qū)位于秦嶺-淮河與南嶺之間。東至大海，西界以武當山、巫山、武陵山、雪峰山等海拔1000m等高線與西南地區(qū)為界，北界以西峽、方城、淮河，蘇北灌溉總渠一線與華北為界，南界以福清、永春、華安、河源、懷集、梧州、平南、忻城一線與華南相接，此線大致相當于1月份均溫10℃～12℃，活動積溫6中溫低濕儲糧區(qū)位于華北、華中、華南和青藏地區(qū)之間。西南部毗鄰緬甸。本區(qū)以秦嶺為北界，大致相當于1月份均溫0℃等溫線。活動積溫4500℃等值線。東界為華中區(qū)的西界，南界從忻城開始沿百色、那城、文山、開遠、景東、潞西北部、梁河至尖高山一線，大致相當于1月份均溫10℃等溫線，活動積溫60007高溫高濕儲糧區(qū)位于我國最南部，包括大陸和島嶼兩部分。包括海南省以及福建省東南部，廣東和廣西的中南部，云南省的南部和西南部。本地區(qū)絕大部分在北回歸線以南，北界即華中和西南儲糧區(qū)的南界，大致相當于1月份均溫10℃等溫線，活動積溫6000℃～6500附件2-2——抽樣過程中所涉及到的相關(guān)計算公式及計算過程一、單階段純隨機抽樣樣本單位數(shù)的計算由2006年全國各縣（區(qū)、市）的糧食產(chǎn)量統(tǒng)計資料顯示，全國共有產(chǎn)糧縣（區(qū)、市）2028個。根據(jù)在7大儲糧生態(tài)區(qū)域32個樣本縣（區(qū)、市）進行的農(nóng)戶儲糧損失預(yù)調(diào)查數(shù)據(jù)，計算得到32個樣本縣（區(qū)、市）損失率平均數(shù)為2.72425%，標準差為1.24197%，191個農(nóng)戶損失率平均數(shù)為2.77503%，標準差為1.74534%。根據(jù)191個農(nóng)戶的有效調(diào)查數(shù)據(jù)資料，按單階段純隨機抽樣方式進行應(yīng)抽取樣本農(nóng)戶數(shù)量的計算如下：（戶）式中：n——樣本容量，t——正態(tài)分布概率度，s——預(yù)調(diào)查農(nóng)戶損失率標準差，r——相對誤差率，——預(yù)調(diào)查農(nóng)戶的平均損失率。但是，這樣計算得到的1688戶農(nóng)戶，理論上將分布在全國2028個縣（區(qū)、市）級單位中，一步到位的直接抽取需要得到全國所有農(nóng)戶的抽樣框（名錄框架），因此，直接抽取很不現(xiàn)實，且樣本點將非常分散，采用隨機抽樣確定這些農(nóng)戶所需的數(shù)據(jù)資料也難以取得。對于這種超大規(guī)模的抽樣調(diào)查，通常是采用多階段抽樣來實施的。二、多階段隨機抽樣中第一階段樣本單位數(shù)的計算第一階段縣級樣本單位采用分層PPS系統(tǒng)抽樣法進行抽取。首先計算出第一階段的抽樣數(shù)目，然后用最優(yōu)分配方法將其分配到7個儲糧生態(tài)區(qū)域（層）中。分層抽樣樣本容量的計算：（層內(nèi)重復抽樣）（層內(nèi)不重復抽樣）式中：——各區(qū)域預(yù)調(diào)查樣本縣損失率方差的平均數(shù)，——各區(qū)域預(yù)調(diào)查樣本縣的平均損失率，N——第一階段縣級單位總數(shù)。各儲糧生態(tài)區(qū)域預(yù)調(diào)查樣本資料——層內(nèi)樣本方差的計算儲糧區(qū)域預(yù)調(diào)查省/自治區(qū)樣本縣平均損失率（%）樣本縣平均損失率的標準差樣本縣平均損失率的方差1西藏3.63600.0109600.000120002內(nèi)蒙2.40300.0039520.000015603遼寧————0.000076604河南1.73000.0096210.000092565江蘇、四川2.37270.0102520.000105106云南2.42000.0070410.000049587廣西————0.00007660平均——————0.00007660(注：3區(qū)、7區(qū)因預(yù)調(diào)查有效數(shù)據(jù)不足，其平均損失率的方差用其它區(qū)域的平均數(shù)代替。)故有：=0.00007660=2.724%(按預(yù)調(diào)查全部32個縣計算)則：考慮到實際抽樣中可能產(chǎn)生樣本重復，預(yù)先放寬5%計算，則第一階段縣級單位抽樣數(shù)目確定為：三、多階段隨機抽樣中第一階段樣本數(shù)的分配樣本的簡單分配不利于減少抽樣誤差，因此，本次調(diào)查采用尼曼最優(yōu)分配法進行第一階段樣本數(shù)的分配。按尼曼最優(yōu)分配法，第h個區(qū)域應(yīng)分配的抽樣數(shù)目為：（h=1,2,…,L）式中：Wh——總體第h層的層權(quán)，Sh——第h層調(diào)查標志的標準差，n——第一階段分層抽樣的總樣本數(shù)，nh——第h層按最優(yōu)分配法分配的樣本數(shù)。這就需要考慮用什么來計算層權(quán)Wh的問題。通常在使用尼曼最優(yōu)分配法時，采用各層總體單元數(shù)占總體全部單元數(shù)的比重，即來表示層權(quán)。在本次調(diào)查中，采用各區(qū)域的縣級單位數(shù)占全部縣級單位數(shù)和比重來計算層權(quán)。因此，各區(qū)域抽樣數(shù)目就取決于各區(qū)域內(nèi)糧食損失率的波動幅度和各區(qū)域縣級單位數(shù)占比的共同影響。結(jié)合本次預(yù)調(diào)查的實際情況來考慮，這種設(shè)計并不十分理想，因為在各區(qū)域中，糧食的生產(chǎn)量和儲藏量差異非常大，而糧食儲藏量大的區(qū)域，理應(yīng)得到適度偏重的觀察。即，在樣本容量的最優(yōu)分配中，應(yīng)考慮各區(qū)域糧食損失量的占比。但由于糧食損失量的占比數(shù)據(jù)難以取得，因此，本次抽樣選擇了與糧食損失量高度相關(guān)的指標糧食生產(chǎn)量占比進行層權(quán)計算。分別計算出按各區(qū)域縣級單位占比作層權(quán)Wh的最優(yōu)分配數(shù)目和按各區(qū)域糧食生產(chǎn)量占比作層權(quán)Wh的最優(yōu)分配數(shù)目。這樣，問題就轉(zhuǎn)化為多指標抽樣的樣本分配問題。由于按不同指標計算的尼曼最優(yōu)分配數(shù)可能是不同的，因此需要進行多指標的綜合。在本方案中采用杰森最優(yōu)平均分配法進行綜合分配。杰森最優(yōu)平均分配法是指在所考慮的所有目標（指標）中，選取最重要的k個，對每個目標（指標）j，計算最優(yōu)分配的層樣本量nˊjh，然后計算其平均值，則：（h=1,2,…,L）按此方法，最終將379個樣本縣（區(qū)、市）分配到7大儲糧生態(tài)區(qū)域。四、第二、第三階段抽樣數(shù)目的確定在多階段抽樣中，主要有三種方式計算各階段的樣本數(shù)目。第一種方式：按費用函數(shù)分配，以三階段分配加以說明。首先建立費用函數(shù)：式中，c為總費用，c1、c2、c3分別為每調(diào)查一個一階單位、二階單位和三階單位的平均費用，在此基礎(chǔ)上，建立拉格朗日極值函數(shù)并求解得到n和m值，即達到分割三階段樣本單位數(shù)的目的。第二種方式：先計算純隨機抽樣的抽樣數(shù)目，然后估計設(shè)計效應(yīng)系數(shù)由此設(shè)計效應(yīng)系數(shù)可估計復雜抽樣的抽樣數(shù)目。對多階段抽樣的設(shè)計效應(yīng)系數(shù)，各專業(yè)文獻均認為在3倍左右為宜。第三種方式：經(jīng)驗法，在我國農(nóng)產(chǎn)量抽樣調(diào)查中應(yīng)用較多。經(jīng)驗法認為，在多階段抽樣中，第一階段的抽樣誤差占總誤差的極大比例（約90%以上），因此，第一階段的樣本數(shù)目應(yīng)偏大，一般應(yīng)占第一階段總體數(shù)量的20%~40%，若總體數(shù)量較大，則可選較小比例；若總體數(shù)量較小，則應(yīng)選較大的比例。第一種分配方法，要求事先掌握比較完整和準確的費用數(shù)據(jù)，所以實用價值不大。對于本次尚未正式開展的調(diào)查而言，不具備相應(yīng)條件，草率估計也是沒有意義的，因此，重點考慮第二種方法和第三種方法。按第二種方法，本次抽樣設(shè)計效應(yīng)系數(shù)為：（倍）這完全滿足抽樣設(shè)計效應(yīng)系數(shù)的要求。按第三種方法，在本次抽樣設(shè)計中，第一階段樣本數(shù)占總體單位數(shù)的比重為：由于本次調(diào)查的總體特別大，幾乎包含全國所有縣級單位，而第一階段樣本單位占比接近20%，因此，基本上滿足第三種方法的要求。在各主要指標滿足規(guī)范要求的基礎(chǔ)上，為了實際工作操作的簡便，本方案按等數(shù)分配的方式分解了第二、三階段的抽樣數(shù)目，即分別為每個樣本縣（區(qū)、市）隨機抽取了4個樣本村，每個樣本村隨機抽取了4個樣本農(nóng)戶，則確定本次調(diào)查最終的農(nóng)戶樣本數(shù)為5968戶。

附件2-3——農(nóng)戶儲糧損失調(diào)查各總體參數(shù)的估計方法由于本抽樣調(diào)查方案是將全國劃分為7大儲糧生態(tài)區(qū)域即子總體（層）來進行的，因此，首先需要分別對7大儲糧生態(tài)區(qū)域的子總體參數(shù)進行估計，然后才能以所估計出的7個子總體參數(shù)按分層抽樣估計方法對全國總體參數(shù)作出最終的估計。一、各生態(tài)區(qū)域子總體參數(shù)的估計方法設(shè)第L個（L=1,2,3,4,5,6,7）子總體擁有Nl個初級單元（縣），每個初級單元擁有Mi個二級單元（村），每個二級單元又擁有Kij個三級單元（農(nóng)戶）。記各階段的抽樣數(shù)目分別為nl、m、k（本方案中m已確定為4個村、k已確定為4個農(nóng)戶），即在第l個生態(tài)區(qū)域中，抽取nl個縣，每個縣抽取4個村，每個村抽取4個農(nóng)戶。在確定了每個生態(tài)區(qū)域的總縣數(shù)、每個縣的村數(shù)及每個村的農(nóng)戶數(shù)的條件下，理論上可以概括出每一階單元被抽中的概率分別為Zi、Zij、Ziju，它們分別滿足：，，這時，可對總體總和（各區(qū)域中某種糧食的損失量）記為：（l=1,2,3,4,5,6,7）其無偏估計為：其方差為：式中：，的一個無偏估計為：（l=1,2,3,4,5,6,7）式中：設(shè)第l個區(qū)域的農(nóng)戶糧食儲藏量的估計量為，則即為第l個生態(tài)區(qū)域農(nóng)戶儲糧損失率的估計值。二、全國總體參數(shù)的估計方法對于前述各生態(tài)區(qū)域的估計值，還應(yīng)進行匯總分析，以得到全國農(nóng)戶儲糧損失率的總體數(shù)據(jù)。以下是各生態(tài)區(qū)域子總體（層）數(shù)據(jù)匯總的方法。本次抽樣調(diào)查農(nóng)戶儲糧損失量的總估計量為，則：記本次抽樣調(diào)查農(nóng)戶糧食儲藏量的總估計量為，則：則即為全國農(nóng)戶儲糧損失率的估計值，其抽樣方差為：三、（各省增補樣本后）以省為子總體的總體參數(shù)估計方法在以省為子總體的抽樣估計中，由于已不再存在各區(qū)域之間分層優(yōu)化的問題，實際上各省的抽樣已轉(zhuǎn)化為三階段線性系統(tǒng)抽樣法，即第一階段的實質(zhì)已是不重復的簡單隨機抽樣。由于全國樣本是按分區(qū)域PPS抽樣決定的，對省而言，這部分樣本可視為一個簡單隨機子樣本。為了避免全國樣本與分省樣本的差異過大，增補部分仍然按簡單不重復抽樣進行抽取，只不過重復抽中全國樣本后，重復的樣本作為無效樣本重新抽樣，直到最終抽完所必需的省樣本數(shù)。對于每一個省區(qū)而言，其三階段簡單隨機不重復抽樣的總體參數(shù)估計式如下所示。記為第i個縣第j個村第u個樣本農(nóng)戶的某種糧食損失量，則該省區(qū)農(nóng)戶儲糧的平均損失量為：其方差為：其無偏估計為：其中f為各階段的抽樣比。

附件2-4——各省（自治區(qū)、直轄市）調(diào)查工作任務(wù)的分配應(yīng)當說明，上述抽樣方案是基于全國農(nóng)戶總體而設(shè)計的，其調(diào)查結(jié)果對于推斷全國農(nóng)戶儲糧損失總量有較高的置信度，但若用其分布在各省（自治區(qū)、直轄市）的樣本作為該?。ㄗ灾螀^(qū)、直轄市）農(nóng)戶儲糧損失總量的推斷，對大部分?。ㄗ灾螀^(qū)、直轄市）而言，其代表性是不夠充分的。為了增強樣本調(diào)查結(jié)果對各省（自治區(qū)、直轄市）的代表性，滿足各?。ㄗ灾螀^(qū)、直轄市）對本?。ㄗ灾螀^(qū)、直轄市）農(nóng)戶儲糧損失總量進行分析和推斷的需要，在前述抽樣基礎(chǔ)上，另外按照上述隨機抽樣原理，以各?。ㄗ灾螀^(qū)、直轄市）為總體按隨機增補策略對各?。ㄗ灾螀^(qū)、直轄市）的調(diào)查樣本數(shù)進行了增補抽樣。由于沒有對每一個?。ㄗ灾螀^(qū)、直轄市）的農(nóng)戶都進行預(yù)抽樣調(diào)查，因此，采用附件2-2中關(guān)于多階段抽樣中第一階段樣本容量分配的第三種方法進行測算。根據(jù)隨機增補策略抽出了各?。ㄗ灾螀^(qū)、直轄市）根據(jù)本省（自治區(qū)、直轄市）需要所進行調(diào)查應(yīng)該增補的調(diào)查樣本縣（區(qū)、市），有些?。ㄗ灾螀^(qū)、直轄市）沒有進行增補抽樣是因為原抽樣本對該省（自治區(qū)、直轄市）已有充分的代表性。具體處理方法是：對于縣級單位數(shù)量在100以上的?。ㄗ灾螀^(qū)、直轄市），抽取20％的縣級單位數(shù)；對于縣級單位數(shù)量在50~100的?。ㄗ灾螀^(qū)、直轄市），抽取25％的縣級單位數(shù)；對于縣級單位數(shù)量在50以下的省（自治區(qū)、直轄市），抽取30％的縣級單位數(shù)。對于各?。ㄗ灾螀^(qū)、直轄市）在原基礎(chǔ)上增補的樣本（縣/區(qū)/市、村、農(nóng)戶）原則上不作強制性調(diào)查要求，有這方面需求（需要增強本省/自治區(qū)/直轄市）樣本對本?。ㄗ灾螀^(qū)、直轄市）農(nóng)戶儲糧損失總量推斷的代表性）的?。ㄗ灾螀^(qū)、直轄市）可以考慮對增補樣本進行調(diào)查。為了便于各省（自治區(qū)、直轄市）調(diào)查工作任務(wù)的劃分，將各?。ㄗ灾螀^(qū)、直轄市）原抽樣本縣/區(qū)/市（各省/自治區(qū)/直轄市必須調(diào)查）和各省（自治區(qū)、直轄市）增補樣本縣/區(qū)/市（根據(jù)本省/自治區(qū)/直轄市需要進行調(diào)查）分別按省統(tǒng)計列于下表中：序號?。ㄗ灾螀^(qū)、直轄市）原抽樣本縣（必須調(diào)查）增補樣本縣（根據(jù)需要進行調(diào)查）縣名數(shù)量縣名數(shù)量1河北省萬全縣、正定縣、無極縣、辛集市、新樂市、灤南縣、遷安市、臨漳縣、磁縣、館陶縣、內(nèi)丘縣、寧晉縣、臨西縣、徐水縣、安新縣、涿州市、平泉縣、鹽山縣、泊頭市、永清縣、棗強縣、景縣、深州市23深澤縣、玉田縣、肥鄉(xiāng)縣、唐縣、雄縣52山西省長子縣、山陰縣、介休市、河津市、洪洞縣、汾陽市6古交市、靈丘縣、長治縣、屯留縣、沁水縣、陵川縣、聞喜縣、永濟市、五寨縣、偏關(guān)縣、鄉(xiāng)寧縣、隰縣、文水縣、嵐縣143內(nèi)蒙古土默特右旗、寧城縣、科爾沁左翼后旗、達拉特旗、五原縣、科爾沁右翼前旗、阿榮旗、扎蘭屯8土默特左旗、鄂托克前旗、海拉爾、新巴爾虎左旗、烏拉特后旗、卓資縣、察哈爾右翼后旗、豐鎮(zhèn)市84遼寧省康平縣、寬甸滿族自治縣、西豐縣、昌圖縣、新民市、普蘭店市、海城市、黑山縣、蓋州市、阜新蒙古族自治縣、大洼縣、朝陽縣、北票市13沒有增補抽樣/5吉林省前郭爾羅斯蒙古族自治縣、鎮(zhèn)賚縣、農(nóng)安縣、九臺市、榆樹市、德惠市、永吉縣、舒蘭市、磐石市、梨樹縣、伊通滿族自治縣、公主嶺市、東豐縣、輝南縣、扶余縣15沒有增補抽樣/6黑龍江省泰來、依蘭、賓縣、木蘭、阿城、雙城、五常、龍江依安、克山、訥河、虎林市、集賢縣、肇州縣、嘉蔭縣、湯原縣、富錦市、海林市、嫩江縣、望奎縣、慶安縣、肇東市22沒有增補抽樣/7江蘇省銅山縣、邳州市、東海縣、漣水縣、盱眙縣、響水縣、射陽縣、東臺市、寶應(yīng)縣、興化市、宿豫區(qū)、泗陽縣、高淳縣、溧陽市、太倉市、如東縣、通州市、江都市、靖江市、姜堰市20沒有增補抽樣/8浙江省蕭山區(qū)、余姚市、嘉善縣、安