工業(yè)抽樣推算操作手冊(cè)

1、規(guī)模以下工業(yè)抽樣調(diào)查推算操作手冊(cè)（年報(bào)和半年報(bào)）年報(bào)和半年報(bào)采用的樣本是一致的，在數(shù)據(jù)處理上的唯一區(qū)別是在年報(bào)調(diào)查中，企業(yè)部分多“實(shí)收資本”等7個(gè)調(diào)查指標(biāo)，下面以年報(bào)數(shù)據(jù)為例來(lái)闡述用stata軟件進(jìn)行總量和方差推算的方法和步驟，在半年報(bào)中，只需要在企業(yè)部分少推算“實(shí)收資本”等7個(gè)調(diào)查指標(biāo)即可。一、數(shù)據(jù)轉(zhuǎn)換用stata/transfer將基層表foxpro數(shù)據(jù)庫(kù)轉(zhuǎn)換成stata數(shù)據(jù)集。二、目錄企業(yè)總量和方差推算步驟第一步，給基層表數(shù)據(jù)集賦予基礎(chǔ)權(quán)數(shù)。通過(guò)將基層表數(shù)據(jù)集和樣本數(shù)據(jù)集進(jìn)行橫向合并來(lái)對(duì)基層表數(shù)據(jù)集賦予基礎(chǔ)權(quán)數(shù)。在進(jìn)行兩個(gè)數(shù)據(jù)集橫向合并之前，必須確保兩個(gè)數(shù)據(jù)集中的觀察記錄是一一對(duì)應(yīng)關(guān)系，

2、換句話說(shuō)就是，通過(guò)變量“qyfrdm”能夠?qū)崿F(xiàn)兩個(gè)數(shù)據(jù)集之間的完全匹配。兩個(gè)數(shù)據(jù)集的觀察記錄是否一一對(duì)應(yīng)，可以從橫向匹配的結(jié)果中看出。use n01_nb.dta, clearsort qyfrdmsave n01_nb.dta, replaceuse yb_mlqy*.dta, clearkeep qyfrdm bs code total_c samp_c weightsort qyfrdmmerge qyfrdm using n01_nb.dtacount if _merge!=3local nonqys=r(n)if nonqys0 display 基層表數(shù)據(jù)集和樣本數(shù)據(jù)集有部分企業(yè)不對(duì)

3、應(yīng)！list qymc qyfrdm _merge if _merge!=3clearelse drop _mergesave mlqy*.dta, replace說(shuō)明：上述命令中包含的“*”表示省代碼，下同；將賦予了基礎(chǔ)權(quán)數(shù)的基層表數(shù)據(jù)集保存成mlqy*.dta是為了使n01_nb.dta不發(fā)生任何改變，便于數(shù)據(jù)管理。第二步，新增一個(gè)變量并命名為dws，作為有效樣本單位和無(wú)效樣本單位 有效樣本單位和無(wú)效樣本單位的定義參見(jiàn)規(guī)模以下工業(yè)抽樣調(diào)查設(shè)計(jì)。的標(biāo)識(shí)變量，用于估計(jì)企業(yè)數(shù)。變量dws的賦值方法為：有效樣本單位賦值為1，無(wú)效樣本單位賦值為0。use mlqy*.dta, cleargen dw

4、s=1replace dws=0 if reason=2save mlqy*.dta, replace第三步，根據(jù)樣本企業(yè)回答情況，對(duì)樣本企業(yè)的基礎(chǔ)權(quán)數(shù)進(jìn)行調(diào)整。如果沒(méi)有無(wú)回答樣本企業(yè)，則不需要進(jìn)行權(quán)數(shù)調(diào)整；如果最終層內(nèi)只有少數(shù)樣本企業(yè)無(wú)回答，則在最終層內(nèi)部進(jìn)行權(quán)數(shù)調(diào)整；如果最終層內(nèi)絕大部分樣本企業(yè)無(wú)回答，則首先需要進(jìn)行最終層合并，然后在合并層內(nèi)進(jìn)行權(quán)數(shù)調(diào)整。權(quán)數(shù)調(diào)整的范圍是最終層（或合并層）內(nèi)的有效樣本。調(diào)整權(quán)數(shù)的具體方法參見(jiàn)規(guī)模以下工業(yè)抽樣調(diào)查方案。為了保存權(quán)數(shù)調(diào)整過(guò)程，需要新增一個(gè)變量（baseweight）保存基礎(chǔ)權(quán)數(shù)，將調(diào)整后的權(quán)數(shù)保存在變量weight中。需要注意的是，只能調(diào)整基層

5、表數(shù)據(jù)集中的權(quán)數(shù)，不能調(diào)整樣本數(shù)據(jù)集中的權(quán)數(shù)，因?yàn)楸敬握{(diào)查中的無(wú)回答企業(yè)，在下次調(diào)查中仍然要作為樣本企業(yè)進(jìn)行調(diào)查，并按照在當(dāng)次調(diào)查中是否回答調(diào)查表決定是否進(jìn)行權(quán)數(shù)調(diào)整。（下同）例一（最終層內(nèi)有少數(shù)樣本企業(yè)無(wú)回答情況舉例）：某省發(fā)現(xiàn)有1個(gè)無(wú)回答樣本企業(yè)，其的基礎(chǔ)權(quán)數(shù)為20，層代碼為10611，該層中共抽取了5個(gè)樣本，假設(shè)剩下4個(gè)樣本企業(yè)都是有回答的有效樣本企業(yè)，則調(diào)整因子為5/4，4個(gè)有回答的樣本企業(yè)的權(quán)數(shù)被調(diào)整為20*5/4=25，同時(shí)，將那個(gè)無(wú)回答樣本企業(yè)的權(quán)數(shù)調(diào)整為“0”。use mlqy*.dta, cleargen baseweight=weightreplace weight=bas

6、eweight*5/4 if code=10611replace weight=0 if code=10611 & reason=1save mlqy*.dta, replace例二（最終層內(nèi)有多數(shù)樣本企業(yè)無(wú)回答情況舉例）：某省發(fā)現(xiàn)4個(gè)無(wú)回答企業(yè)，都來(lái)自于同一層，層代碼為10611，其基礎(chǔ)權(quán)數(shù)都為20，該層中共抽取了5個(gè)樣本企業(yè)，此時(shí)，需要將該層中剩下的1個(gè)有回答的有效樣本企業(yè)與臨近的層合并后，再進(jìn)行權(quán)數(shù)調(diào)整。比如有一個(gè)臨近的層代碼為10621，該層中共抽取了5個(gè)樣本企業(yè)，基礎(chǔ)權(quán)數(shù)為10，但其中有1個(gè)樣本企業(yè)消亡，屬于無(wú)效樣本企業(yè)。將層代碼10611替換成10621，然后在合并層10621內(nèi)對(duì)

7、5（=1+4）個(gè)有回答的有效樣本企業(yè)的權(quán)數(shù)進(jìn)行調(diào)整，調(diào)整因子為，原來(lái)屬于層10611的1個(gè)有回答的有效樣本企業(yè)的權(quán)數(shù)調(diào)整為20*14/6=46.67，原來(lái)屬于層10621的4個(gè)有回答的有效樣本企業(yè)的權(quán)數(shù)調(diào)整為10*14/6=23.33。同時(shí)，將4個(gè)無(wú)回答樣本企業(yè)的權(quán)數(shù)調(diào)整為“0”。由于設(shè)置fpc要求fpc標(biāo)識(shí)變量（total_c）在同一層內(nèi)取值相同，因此，需要對(duì)合并層內(nèi)觀察記錄的fpc標(biāo)識(shí)變量（total_c）進(jìn)行調(diào)整，調(diào)整為合并前兩個(gè)層的total_c之和，即20*5+10*5=150。use mlqy*.dta, cleargen baseweight=weightreplace code

8、=10621 if code=10611replace weight=baseweight*14/6 if code=10621 and dws=1replace weight=0 if code=10621 & reason=1gen basetotal=total_creplace total_c=150 if code=10621save mlqy*.dta, replace第四步，在基層數(shù)據(jù)中設(shè)置抽樣設(shè)計(jì)信息。use mlqy*.dta, clearsvyset strata codesvyset pweight weightsvyset fpc total_csave mlqy*.d

9、ta, replace說(shuō)明：是否設(shè)置“有限總體校正系數(shù)（fpc）”只會(huì)影響方差估計(jì)值的大小。當(dāng)每個(gè)最終層的抽樣比都比較小，即每個(gè)最終層的權(quán)數(shù)都比較大時(shí)，是否設(shè)置“有限總體校正系數(shù)（fpc）”對(duì)推算結(jié)果的影響不大，但是，當(dāng)有些最終層的抽樣比比較大，特別地，當(dāng)有全選層，即抽樣比等于1，權(quán)數(shù)等于1 時(shí)，需要設(shè)置“有限總體校正系數(shù)（fpc）”，否則會(huì)高估方差。因此，統(tǒng)一要求設(shè)置“有限總體校正系數(shù)（fpc）”。在stata中，設(shè)置“有限總體校正系數(shù)（fpc）”有兩種方法，一是將最終層的總體單位數(shù)（total_c）設(shè)置成“有限總體校正系數(shù)（fpc）”，二是將最終層的抽樣比設(shè)置成“有限總體校正系數(shù)（fpc）

10、”，兩種設(shè)置方法得到的結(jié)果是一致的。要求統(tǒng)一采用第一種方法。第五步，進(jìn)行總量和方差估計(jì)（包括估計(jì)企業(yè)數(shù)及其方差）。首先，需要將各調(diào)查指標(biāo)的計(jì)量單位轉(zhuǎn)換成萬(wàn)元，a01（全部從業(yè)人員期末數(shù)）除外；第二，計(jì)算各指標(biāo)的樣本匯總值；第三，估計(jì)企業(yè)數(shù)及各調(diào)查指標(biāo)的總量和方差。use mlqy*.dta, clearfor var a*:replace x=x/10replace a01=a01*10egen hz00=sum(dws)egen hz01=sum(a01)egen hz02=sum(a02)egen hz03=sum(a03)egen hz04=sum(a04)egen hz05=sum(a

11、05)egen hz06=sum(a06)egen hz07=sum(a07)egen hz08=sum(a08)egen hz09=sum(a09)egen hz10=sum(a10)egen hz11=sum(a11)egen hz12=sum(a12)egen hz13=sum(a13)egen hz14=sum(a14)egen hz15=sum(a15)svytotal dws a*第六步，制綜合表方法一：首先，編寫(xiě)stata命令文件（do-file），將推算結(jié)果保存到一個(gè)stata數(shù)據(jù)集中，然后將此stata數(shù)據(jù)集用stata/transfer轉(zhuǎn)換成foxpro數(shù)據(jù)庫(kù)文件，導(dǎo)入到規(guī)

12、模以下工業(yè)抽樣調(diào)查數(shù)據(jù)管理系統(tǒng)中，制表輸出。方法二：用excel制表。具體方法為：在stata結(jié)果窗口中選中推算結(jié)果表，在edit菜單中選用“copy table”功能復(fù)制推算結(jié)果表，然后打開(kāi)一個(gè)excel文件，執(zhí)行“粘貼”操作，即將stata的推算結(jié)果拷貝到excel中了，最后在excel文件中計(jì)算相對(duì)誤差（stata推算結(jié)果中沒(méi)有計(jì)算相對(duì)誤差的參數(shù)選項(xiàng)），并編輯制表。舉例說(shuō)明將推算結(jié)果保存在一個(gè)stata數(shù)據(jù)集中的過(guò)程。首先，按照綜合表的變量結(jié)構(gòu)構(gòu)造一個(gè)數(shù)據(jù)集；然后，將系統(tǒng)自動(dòng)保存在內(nèi)存距陣變量e(est)、e(v_db)等中的推算結(jié)果，逐個(gè)地存放到數(shù)據(jù)集的指定位置。use mlqy*.d

13、ta, clearfor var a*:replace x=x/10replace a01=a01*10egen hz00=sum(dws)egen hz01=sum(a01)egen hz15=sum(a15)svytotal dws a*matrix point=e(est)matrix variance=e(v_db)replace smc=（省名稱(chēng)）gen str2 sdm=（省代碼）keep smc sdm hz00 hz01 hz02 hz03 hz04 hz05 hz06 hz07 hz08 hz09 hz10 hz11 hz12 hz13 hz14 hz15order smc

14、sdm hz00 hz01 hz02 hz03 hz04 hz05 hz06 hz07 hz08 hz09 hz10 hz11 hz12 hz13 hz14 hz15gen a00=.gen jfc00=.gen xdw00=.gen a01=.gen jfc01=.gen xdw01=.gen a02=.gen jfc02=.gen xdw02=.gen a03=.gen jfc03=.gen xdw03=.gen a04=.gen jfc04=.gen xdw04=.gen a05=.gen jfc05=.gen xdw05=.gen a06=.gen jfc06=.gen xdw06=.

15、gen a07=.gen jfc07=.gen xdw07=.gen a08=.gen jfc08=.gen xdw08=.gen a09=.gen jfc09=.gen xdw09=.gen a10=.gen jfc10=.gen xdw10=.gen a11=.gen jfc11=.gen xdw11=.gen a12=.gen jfc12=.gen xdw12=.gen a13=.gen jfc13=.gen xdw13=.gen a14=.gen jfc14=.gen xdw14=.gen a15=.gen jfc15=.gen xdw15=.format a00 hz00 %10.0

16、fformat a01 a02 a03 a04 a05 a06 a07 a08 a09 a10 a11 a12 a13 a14 a15 %15.2fformat hz01 hz02 hz03 hz04 hz05 hz06 hz07 hz08 %10.2fformat hz09 hz10 hz11 hz12 hz13 hz14 hz15 %10.2fformat jfc00 jfc01 jfc02 jfc03 jfc04 jfc05 jfc06 jfc07 jfc08 %15.2fformat jfc09 jfc10 jfc11 jfc12 jfc13 jfc14 jfc15 %15.2ffor

17、mat xdw00 xdw01 xdw02 xdw03 xdw04 xdw05 xdw06 xdw07 xdw08 %6.2fformat xdw09 xdw10 xdw11 xdw12 xdw13 xdw14 xdw15 %6.2fkeep in 1forvalues i = 0/15 if i 0 display 基層表數(shù)據(jù)集和樣本數(shù)據(jù)集有部分企業(yè)不對(duì)應(yīng)！ list qymc qyfrdm _merge if _merge!=3clearelse drop _mergesave mlqy*.dta, replacegen dws=1replace dws=0 if reason=2svys

18、et strata codesvyset pweight weightsvyset fpc total_csave mlqy*.dta, replacefor var a*:replace x=x/10replace a01=a01*10egen hz00=sum(dws)egen hz01=sum(a01)egen hz15=sum(a15)svytotal dws a*matrix point=e(est)matrix variance=e(v_db)replace smc=（省名稱(chēng)）gen str2 sdm=（省代碼）keep smc sdm hz00 hz01 hz02 hz03 hz

19、04 hz05 hz06 hz07 hz08 hz09 hz10 hz11 hz12 hz13 hz14 hz15order smc sdm hz00 hz01 hz02 hz03 hz04 hz05 hz06 hz07 hz08 hz09 hz10 hz11 hz12 hz13 hz14 hz15gen a00=.gen jfc00=.gen xdw00=.gen a01=.gen jfc01=.gen xdw01=.gen a15=.gen jfc15=.gen xdw15=.format a00 hz00 %10.0fformat a01 a02 a03 a04 a05 a06 a07

20、a08 a09 a10 a11 a12 a13 a14 a15 %13.2fformat hz01 hz02 hz03 hz04 hz05 hz06 hz07 hz08 %10.2fformat hz09 hz10 hz11 hz12 hz13 hz14 hz15 %10.2fformat jfc00 jfc01 jfc02 jfc03 jfc04 jfc05 jfc06 jfc07 jfc08 %15.2fformat jfc09 jfc10 jfc11 jfc12 jfc13 jfc14 jfc15 %15.2fformat xdw00 xdw01 xdw02 xdw03 xdw04 xd

21、w05 xdw06 xdw07 xdw08 %6.2fformat xdw09 xdw10 xdw11 xdw12 xdw13 xdw14 xdw15 %6.2fkeep in 1forvalues i = 0/15 if i 10 local j=i+1replace a0i=point1,j in 1replace jfc0i=sqrt(variancej,j) in 1replace xdw0i=1.96*jfc0i/a0i*100else local j=i+1replace ai=point1,j in 1replace jfci=sqrt(variancej,j) in 1repl

22、ace xdwi=1.96*jfci/ai*100save zhbml_nb*.dta, replaceclear2、個(gè)體經(jīng)營(yíng)工業(yè)單位子總體clearset more offcd （數(shù)據(jù)集存放路徑）use n04_nb.dta, clearsort dmsave n04_nb.dta, replaceuse yb_gtyjd*.dta, clearkeep dm pseudo code total_c samp_c weightsort dmmerge dm using n04_nb.dtadrop _mergesave gt*.dta, replacegen dws=1replace dws

23、=0 if reason=1svyset strata codesvyset pweight weightsvyset psu dmsvyset fpc total_csave gt*.dta, replacefor var a0*:replace x=x/10replace a01=a01*10egen hz00=sum(dws)egen hz01=sum(a01)egen hz02=sum(a02)egen hz04=sum(a04)svytotal dws a0*matrix point=e(est)matrix variance=e(v_db)replace smc=（省名稱(chēng)）gen

24、str2 sdm=（省代碼）keep smc sdm hz00 hz01 hz02 hz04order smc sdm hz00 hz01 hz02 hz04gen a00=.gen jfc00=.gen xdw00=.gen a01=.gen jfc01=.gen xdw01=.gen a02=.gen jfc02=.gen xdw02=.gen a04=.gen jfc04=.gen xdw04=.format a00 hz00 %10.0fformat a01 a02 a04 %15.2fformat hz01 hz02 hz04 %10.2fformat jfc00 jfc01 jfc

25、02 jfc04 %15.2fformat xdw00 xdw01 xdw02 xdw04 %6.2fkeep in 1forvalues i = 0/3 if i3 local j=i+1 replace a0i=point1,j in 1 replace jfc0i=sqrt(variancej,j) in 1 replace xdw0i=1.96*jfc0i/a0i*100 else local j=i+1 replace a0j=point1,j in 1 replace jfc0j=sqrt(variancej,j) in 1 replace xdw0j=1.96*jfc0j/a0j

26、*100 save zhbgt_nb*.dta, replaceclear3、非目錄企業(yè)clearset more offcd （數(shù)據(jù)集存放路徑）use n02_nb.dta, clearsort dmsave n02_nb.dta, replaceuse yb_gtyjd*.dta, clearkeep dm pseudo code total_c samp_c weightsort dmmerge dm using n02_nb.dtareplace reason=1 if reason=.for var a*: replace x=0 if x=.drop _mergesave fmlq

27、y*.dta, replacegen dws=1replace dws=0 if reason=1svyset strata codesvyset pweight weightsvyset psu dmsvyset fpc total_csave fmlqy*.dta, replacefor var a*:replace x=x/10replace a01=a01*10egen hz00=sum(dws)egen hz01=sum(a01)egen hz15=sum(a15)svytotal dws a*matrix point=e(est)matrix variance=e(v_db)rep

28、lace smc=（省名稱(chēng)）gen str2 sdm=（省代碼）keep smc sdm hz00 hz01 hz02 hz03 hz04 hz05 hz06 hz07 hz08 hz09 hz10 hz11 hz12 hz13 hz14 hz15order smc sdm hz00 hz01 hz02 hz03 hz04 hz05 hz06 hz07 hz08 hz09 hz10 hz11 hz12 hz13 hz14 hz15gen a00=.gen jfc00=.gen xdw00=.gen a01=.gen jfc01=.gen xdw01=.gen a15=.gen jfc15=.g

29、en xdw15=.format a00 hz00 %10.0fformat a01 a02 a03 a04 a05 a06 a07 a08 a09 a10 a11 a12 a13 a14 a15 %13.2fformat hz01 hz02 hz03 hz04 hz05 hz06 hz07 hz08 %10.2fformat hz09 hz10 hz11 hz12 hz13 hz14 hz15 %10.2fformat jfc00 jfc01 jfc02 jfc03 jfc04 jfc05 jfc06 jfc07 jfc08 %15.2fformat jfc09 jfc10 jfc11 jf

30、c12 jfc13 jfc14 jfc15 %15.2fformat xdw00 xdw01 xdw02 xdw03 xdw04 xdw05 xdw06 xdw07 xdw08 %6.2fformat xdw09 xdw10 xdw11 xdw12 xdw13 xdw14 xdw15 %6.2fkeep in 1forvalues i = 0/15 if i 10 local j=i+1replace a0i=point1,j in 1replace jfc0i=sqrt(variancej,j) in 1replace xdw0i=1.96*jfc0i/a0i*100else local j

31、=i+1replace ai=point1,j in 1replace jfci=sqrt(variancej,j) in 1replace xdwi=1.96*jfci/ai*100save zhbfml_nb*.dta, replaceclear4、工業(yè)企業(yè)子總體clearset more offcd （數(shù)據(jù)集存放路徑）use fmlqy*.dta, cleargen str12 ppsu=dmsave temp1.dta, replaceuse mlqy*.dta, cleargen str12 ppsu=qyfrdmappend using temp1.dtasave qy*.dta,

32、 replaceerase temp1.dtasvyset strata codesvyset pweight weightsvyset psu ppsusvyset fpc total_csave qy *.dta, replacefor var a*:replace x=x/10replace a01=a01*10egen hz00=sum(dws)egen hz01=sum(a01)egen hz15=sum(a15)svytotal dws a*matrix point=e(est)matrix variance=e(v_db)replace smc=（省名稱(chēng)）gen str2 sdm

33、=（省代碼）keep smc sdm hz00 hz01 hz02 hz03 hz04 hz05 hz06 hz07 hz08 hz09 hz10 hz11 hz12 hz13 hz14 hz15order smc sdm hz00 hz01 hz02 hz03 hz04 hz05 hz06 hz07 hz08 hz09 hz10 hz11 hz12 hz13 hz14 hz15gen a00=.gen jfc00=.gen xdw00=.gen a01=.gen jfc01=.gen xdw01=.gen a15=.gen jfc15=.gen xdw15=.format a00 hz00

34、%10.0fformat a01 a02 a03 a04 a05 a06 a07 a08 a09 a10 a11 a12 a13 a14 a15 %13.2fformat hz01 hz02 hz03 hz04 hz05 hz06 hz07 hz08 %10.2fformat hz09 hz10 hz11 hz12 hz13 hz14 hz15 %10.2fformat jfc00 jfc01 jfc02 jfc03 jfc04 jfc05 jfc06 jfc07 jfc08 %15.2fformat jfc09 jfc10 jfc11 jfc12 jfc13 jfc14 jfc15 %15.

35、2fformat xdw00 xdw01 xdw02 xdw03 xdw04 xdw05 xdw06 xdw07 xdw08 %6.2fformat xdw09 xdw10 xdw11 xdw12 xdw13 xdw14 xdw15 %6.2fkeep in 1forvalues i = 0/15 if i 10 local j=i+1replace a0i=point1,j in 1replace jfc0i=sqrt(variancej,j) in 1replace xdw0i=1.96*jfc0i/a0i*100else local j=i+1replace ai=point1,j in 1replace jfci=sqrt(variancej,j) in 1replace xdwi=1.96*jfci/ai*100save zhbqy_nb*.dta, replaceclear5、規(guī)模以下工業(yè)總體clearset more offcd （數(shù)據(jù)集存放