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

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

2、換句話說就是，通過變量“qyfrdm”能夠?qū)崿F(xiàn)兩個數(shù)據(jù)集之間的完全匹配。兩個數(shù)據(jù)集的觀察記錄是否一一對應，可以從橫向匹配的結(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è)不對

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

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

5、表數(shù)據(jù)集中的權(quán)數(shù)，不能調(diào)整樣本數(shù)據(jù)集中的權(quán)數(shù)，因為本次調(diào)查中的無回答企業(yè)，在下次調(diào)查中仍然要作為樣本企業(yè)進行調(diào)查，并按照在當次調(diào)查中是否回答調(diào)查表決定是否進行權(quán)數(shù)調(diào)整。（下同）例一（最終層內(nèi)有少數(shù)樣本企業(yè)無回答情況舉例）：某省發(fā)現(xiàn)有1個無回答樣本企業(yè)，其的基礎權(quán)數(shù)為20，層代碼為10611，該層中共抽取了5個樣本，假設剩下4個樣本企業(yè)都是有回答的有效樣本企業(yè)，則調(diào)整因子為5/4，4個有回答的樣本企業(yè)的權(quán)數(shù)被調(diào)整為20*5/4=25，同時，將那個無回答樣本企業(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è)無回答情況舉例）：某省發(fā)現(xiàn)4個無回答企業(yè)，都來自于同一層，層代碼為10611，其基礎權(quán)數(shù)都為20，該層中共抽取了5個樣本企業(yè)，此時，需要將該層中剩下的1個有回答的有效樣本企業(yè)與臨近的層合并后，再進行權(quán)數(shù)調(diào)整。比如有一個臨近的層代碼為10621，該層中共抽取了5個樣本企業(yè)，基礎權(quán)數(shù)為10，但其中有1個樣本企業(yè)消亡，屬于無效樣本企業(yè)。將層代碼10611替換成10621，然后在合并層10621內(nèi)對

7、5（=1+4）個有回答的有效樣本企業(yè)的權(quán)數(shù)進行調(diào)整，調(diào)整因子為，原來屬于層10611的1個有回答的有效樣本企業(yè)的權(quán)數(shù)調(diào)整為20*14/6=46.67，原來屬于層10621的4個有回答的有效樣本企業(yè)的權(quán)數(shù)調(diào)整為10*14/6=23.33。同時，將4個無回答樣本企業(yè)的權(quán)數(shù)調(diào)整為“0”。由于設置fpc要求fpc標識變量（total_c）在同一層內(nèi)取值相同，因此，需要對合并層內(nèi)觀察記錄的fpc標識變量（total_c）進行調(diào)整，調(diào)整為合并前兩個層的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ù)中設置抽樣設計信息。use mlqy*.dta, clearsvyset strata codesvyset pweight weightsvyset fpc total_csave mlqy*.d

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

10、”，兩種設置方法得到的結(jié)果是一致的。要求統(tǒng)一采用第一種方法。第五步，進行總量和方差估計（包括估計企業(yè)數(shù)及其方差）。首先，需要將各調(diào)查指標的計量單位轉(zhuǎn)換成萬元，a01（全部從業(yè)人員期末數(shù)）除外；第二，計算各指標的樣本匯總值；第三，估計企業(yè)數(shù)及各調(dià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*第六步，制綜合表方法一：首先，編寫stata命令文件（do-file），將推算結(jié)果保存到一個stata數(shù)據(jù)集中，然后將此stata數(shù)據(jù)集用stata/transfer轉(zhuǎn)換成foxpro數(shù)據(jù)庫文件，導入到規(guī)

12、模以下工業(yè)抽樣調(diào)查數(shù)據(jù)管理系統(tǒng)中，制表輸出。方法二：用excel制表。具體方法為：在stata結(jié)果窗口中選中推算結(jié)果表，在edit菜單中選用“copy table”功能復制推算結(jié)果表，然后打開一個excel文件，執(zhí)行“粘貼”操作，即將stata的推算結(jié)果拷貝到excel中了，最后在excel文件中計算相對誤差（stata推算結(jié)果中沒有計算相對誤差的參數(shù)選項），并編輯制表。舉例說明將推算結(jié)果保存在一個stata數(shù)據(jù)集中的過程。首先，按照綜合表的變量結(jié)構(gòu)構(gòu)造一個數(shù)據(jù)集；然后，將系統(tǒng)自動保存在內(nèi)存距陣變量e(est)、e(v_db)等中的推算結(jié)果，逐個地存放到數(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=（省名稱）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è)不對應！ 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=（省名稱）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、個體經(jī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=（省名稱）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=（省名稱）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=（省名稱）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ù)集存放