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

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

2、換句話說就是，通過變量“qyfrdm”能夠?qū)崿F(xiàn)兩個(gè)數(shù)據(jù)集之間的完全匹配。兩個(gè)數(shù)據(jù)集的觀察記錄是否一一對應(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 nonqys'>0 display "基層表數(shù)

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

4、lqy*.dta, cleargen dws=1replace dws=0 if reason>=2save mlqy*.dta, replace第三步，根據(jù)樣本企業(yè)回答情況，對樣本企業(yè)的基礎(chǔ)權(quán)數(shù)進(jìn)行調(diào)整。如果沒有無回答樣本企業(yè)，則不需要進(jìn)行權(quán)數(shù)調(diào)整；如果最終層內(nèi)只有少數(shù)樣本企業(yè)無回答，則在最終層內(nèi)部進(jìn)行權(quán)數(shù)調(diào)整；如果最終層內(nèi)絕大部分樣本企業(yè)無回答，則首先需要進(jìn)行最終層合并，然后在合并層內(nèi)進(jìn)行權(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)整過程，需要新增一個(gè)變量（baseweight）保存基礎(chǔ)權(quán)數(shù)，將調(diào)整后的權(quán)數(shù)保

5、存在變量weight中。需要注意的是，只能調(diào)整基層表數(shù)據(jù)集中的權(quán)數(shù)，不能調(diào)整樣本數(shù)據(jù)集中的權(quán)數(shù)，因?yàn)楸敬握{(diào)查中的無回答企業(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è)無回答情況舉例）：某省發(fā)現(xiàn)有1個(gè)無回答樣本企業(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è)無回答樣本企業(yè)的權(quán)數(shù)調(diào)整為“0”。use mlqy*.dta, cleargen baseweight

6、=weightreplace weight=baseweight*5/4 if code=10611replace weight=0 if code=10611 & reason=1save mlqy*.dta, replace例二（最終層內(nèi)有多數(shù)樣本企業(yè)無回答情況舉例）：某省發(fā)現(xiàn)4個(gè)無回答企業(yè)，都來自于同一層，層代碼為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è)消亡，屬于無效樣本企業(yè)。將層

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

8、baseweight=weightreplace code=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 wei

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

10、fpc）”，二是將最終層的抽樣比設(shè)置成“有限總體校正系數(shù)（fpc）”，兩種設(shè)置方法得到的結(jié)果是一致的。要求統(tǒng)一采用第一種方法。第五步，進(jìn)行總量和方差估計(jì)（包括估計(jì)企業(yè)數(shù)及其方差）。首先，需要將各調(diào)查指標(biāo)的計(jì)量單位轉(zhuǎ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)

11、egen hz04=sum(a04)egen hz05=sum(a05)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é)果保存到一個(gè)STATA數(shù)據(jù)集中，然后將此STATA數(shù)據(jù)集用

12、Stata/Transfer轉(zhuǎn)換成Foxpro數(shù)據(jù)庫文件，導(dǎo)入到規(guī)模以下工業(yè)抽樣調(diào)查數(shù)據(jù)管理系統(tǒng)中，制表輸出。方法二：用EXCEL制表。具體方法為：在STATA結(jié)果窗口中選中推算結(jié)果表，在EDIT菜單中選用“copy table”功能復(fù)制推算結(jié)果表，然后打開一個(gè)EXCEL文件，執(zhí)行“粘貼”操作，即將STATA的推算結(jié)果拷貝到EXCEL中了，最后在EXCEL文件中計(jì)算相對誤差（STATA推算結(jié)果中沒有計(jì)算相對誤差的參數(shù)選項(xiàng)），并編輯制表。舉例說明將推算結(jié)果保存在一個(gè)STATA數(shù)據(jù)集中的過程。首先，按照綜合表的變量結(jié)構(gòu)構(gòu)造一個(gè)數(shù)據(jù)集；然后，將系統(tǒng)自動保存在內(nèi)存距陣變量e(est)、e(V_db)等

13、中的推算結(jié)果，逐個(gè)地存放到數(shù)據(jù)集的指定位置。use mlqy*.dta, 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 hz

14、06 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 a02=.gen jfc02=.gen xdw02=.gen a03=.gen jfc03=.gen xdw03=.gen a04=.gen jfc04=.gen xdw04=.gen a

15、05=.gen jfc05=.gen xdw05=.gen a06=.gen jfc06=.gen xdw06=.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 xd

16、w14=.gen a15=.gen jfc15=.gen xdw15=.format a00 hz00 %10.0fformat 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.2ff

17、ormat jfc09 jfc10 jfc11 jfc12 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(varia

18、ncej',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 zhbml_nb*.dta, replaceclear注：當(dāng)用于推算的數(shù)據(jù)集中某個(gè)或某些估計(jì)變量含有缺省值（.）時(shí)，內(nèi)部函數(shù)e(V_db)只保存方差估計(jì)值

19、，是一個(gè)1，n距陣；當(dāng)數(shù)據(jù)集中所有估計(jì)變量都不含有缺省值時(shí)，內(nèi)部函數(shù)e(V_db)保存的是n個(gè)估計(jì)變量的協(xié)方差估計(jì)值，是一個(gè)n，n距陣。因此，統(tǒng)一要求用于推算的數(shù)據(jù)集中不能包含缺省值。如果發(fā)現(xiàn)有缺省值存在，則可以用以下命令將缺省值替換成0。for var a*: replace X=0 if X=.三、個(gè)體經(jīng)營工業(yè)單位總量和方差推算步驟第一步，給基層表數(shù)據(jù)集賦予基礎(chǔ)權(quán)數(shù)。use n04_nb.dta, clearsort dmsave n04_nb.dta, replaceuse yb_gtyjd*.dta, clearkeep dm pseudo code total_c samp_c we

20、ightsort dmmerge dm using n04_nb.dtadrop _mergesave gt*.dta, replace第二步，新增一個(gè)變量并命名為dws，用于估計(jì)個(gè)體經(jīng)營工業(yè)單位數(shù)。變量dws的賦值方法為：對有調(diào)查資料的個(gè)體經(jīng)營工業(yè)單位賦值為1，對沒有個(gè)體經(jīng)營工業(yè)單位的樣本PSU觀察記錄賦值為0。use gt*.dta, cleargen dws=1replace dws=0 if reason=1save gt*.dta, replace第三步，根據(jù)樣本PSU和個(gè)體經(jīng)營工業(yè)單位的回答情況，對個(gè)體經(jīng)營工業(yè)單位的基礎(chǔ)權(quán)數(shù)進(jìn)行調(diào)整。如果某個(gè)樣本PSU內(nèi)有部分無回答個(gè)體經(jīng)營工業(yè)單位

21、，則需要在該樣本PSU內(nèi)部，對其他有回答個(gè)體經(jīng)營工業(yè)單位的權(quán)數(shù)進(jìn)行調(diào)整；如果最終層中只有少數(shù)無回答樣本PSU 如果某個(gè)樣本PSU內(nèi)所有的個(gè)體經(jīng)營工業(yè)單位都不回答調(diào)查表，則稱該樣本PSU為無回答樣本PSU。，則需要在最終層內(nèi)對其他有回答樣本PSU的權(quán)數(shù)進(jìn)行調(diào)整；如果某個(gè)最終層中絕大多數(shù)的樣本PSU無回答，則首先需要進(jìn)行最終層合并，然后在合并層內(nèi)對其他有回答樣本PSU進(jìn)行權(quán)數(shù)調(diào)整；如果沒有無回答個(gè)體經(jīng)營工業(yè)單位和樣本PSU，則不需要進(jìn)行權(quán)數(shù)調(diào)整。權(quán)數(shù)調(diào)整的范圍是最終層（或合并層）內(nèi)的有效樣本。調(diào)整權(quán)數(shù)的具體方法參見規(guī)模以下工業(yè)抽樣調(diào)查方案?！?”。use gt*.dta, cleargen bas

22、eweight=weightreplace weight=baseweight*5/4 if dm=""replace weight=0 if code="" & reason 根據(jù)規(guī)模以下工業(yè)抽樣調(diào)查數(shù)據(jù)管理系統(tǒng)中的規(guī)定，個(gè)體經(jīng)營工業(yè)單位和樣本PSU沒有填報(bào)調(diào)查表的原因分為三種：1、該樣本PSU中沒有個(gè)體經(jīng)營工業(yè)單位；2、該樣本PSU不回答調(diào)查表，為無回答樣本PSU；3、該個(gè)體經(jīng)營工業(yè)單位不回答調(diào)查表，為無回答個(gè)體經(jīng)營工業(yè)單位。=3save gt*.dta, replace例二（最終層內(nèi)有少數(shù)樣本PSU無回答情況舉例）：某省發(fā)現(xiàn)有1個(gè)無回答樣本

23、PSU，其的基礎(chǔ)權(quán)數(shù)為20，層代碼為111，該層中共抽取了5個(gè)樣本PSU，假設(shè)剩下4個(gè)樣本PSU都是有效樣本PSU，則調(diào)整因子為5/4，4個(gè)有效樣本PSU的權(quán)數(shù)，即這4個(gè)有效樣本PSU內(nèi)所有個(gè)體經(jīng)營工業(yè)單位的權(quán)數(shù)被調(diào)整為20*5/4=25，同時(shí)，將那個(gè)無回答樣本PSU的權(quán)數(shù)調(diào)整為“0”。use gt*.dta, cleargen baseweight=weightreplace weight=baseweight*5/4 if code=111replace weight=0 if code=111 & reason=2save gt*.dta, replace例三（最終層內(nèi)有較多樣本

24、PSU無回答情況舉例）：某省發(fā)現(xiàn)4個(gè)無回答樣本PSU，都來自于同一層，層代碼為111，其基礎(chǔ)權(quán)數(shù)都為20，該層中共抽取了5個(gè)樣本PSU，此時(shí)，需要將該層中剩下的1個(gè)有回答樣本PSU與臨近的層合并后，再對合并層中有回答樣本PSU的權(quán)數(shù)進(jìn)行調(diào)整。比如有一個(gè)臨近的層代碼為112，該層中共抽取了5個(gè)樣本PSU，基礎(chǔ)權(quán)數(shù)為10，但其中有1個(gè)樣本PSU沒有個(gè)體經(jīng)營工業(yè)單位，屬于無效樣本PSU。首先將層代碼111替換成112，然后在合并層112內(nèi)對5（=1+4）個(gè)有回答的有效樣本PSU的權(quán)數(shù)進(jìn)行調(diào)整，調(diào)整因子為，原來屬于層111的1個(gè)樣本PSU的權(quán)數(shù)調(diào)整為20*14/6=46.67，原來屬于層112的4個(gè)樣

25、本PSU的權(quán)數(shù)調(diào)整為10*14/6=23.33。同時(shí)，將4個(gè)無回答樣本PSU的權(quán)數(shù)調(diào)整為“0”。由于設(shè)置fpc要求fpc標(biāo)識變量（total_c）在同一層內(nèi)取值相同，因此，需要對合并層內(nèi)觀察記錄的fpc標(biāo)識變量（total_c）進(jìn)行調(diào)整，調(diào)整為合并前兩個(gè)層的total_c之和，即20*5+10*5=150。use gt*.dta, cleargen baseweight=weightreplace code=112 if code=111replace weight=baseweight*14/6 if code=112 and dws=1replace weight=0 if code=11

26、2 & reason=2gen basetotal=total_creplace total_c=150 if code=112save gt*.dta, replace第四步，在基層數(shù)據(jù)中設(shè)置抽樣設(shè)計(jì)信息。use gt*.dta, clearsvyset strata codesvyset pweight weightsvyset psu dmsvyset fpc total_csave gt*.dta, replace第五步，進(jìn)行總量和方差估計(jì)（包括估計(jì)個(gè)體經(jīng)營工業(yè)單位數(shù)及其方差）。首先，需要將各調(diào)查指標(biāo)的計(jì)量單位轉(zhuǎn)換成萬元，a01（全部從業(yè)人員期末數(shù)）除外；第二，計(jì)算各指標(biāo)的樣本

27、匯總值；第三，估計(jì)單位數(shù)及各調(diào)查指標(biāo)的總量和方差。use gt*.dta, clearfor 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*第六步，制綜合表參見目錄企業(yè)部分。四、非目錄企業(yè)總量和方差推算步驟第一步，給基層表數(shù)據(jù)集賦予基礎(chǔ)權(quán)數(shù)。use n02_nb.dta, clearsort dmsave n02_nb.dta, replaceuse yb_gtyjd*.dta, cle

28、arkeep 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 fmlqy*.dta, replace第二步，新增一個(gè)變量并命名為dws，用于估計(jì)非目錄企業(yè)數(shù)。變量dws的賦值方法為：對有調(diào)查資料的非目錄企業(yè)賦值為1，對沒有非目錄企業(yè)的樣本PSU觀察記錄賦值為0。use fmlqy*.dta, cleargen dws=1replace dws=0 if reas

29、on 根據(jù)規(guī)模以下工業(yè)抽樣調(diào)查數(shù)據(jù)管理系統(tǒng)中的規(guī)定，非目錄企業(yè)和樣本PSU沒有填報(bào)調(diào)查表的原因分為三種：1、該樣本PSU中沒有非目錄企業(yè)；或者，該非目錄企業(yè)變成無效樣本企業(yè)，即撤消、被合并、停產(chǎn)且找不到人、轉(zhuǎn)產(chǎn)成為非工業(yè)企業(yè)、已經(jīng)納入規(guī)模以上工業(yè)統(tǒng)計(jì)等； 2、該樣本PSU不回答調(diào)查表，為無回答樣本PSU；3、該非目錄企業(yè)不回答調(diào)查表，為無回答非目錄企業(yè)。=1save fmlqy *.dta, replace第三步，根據(jù)樣本PSU和非目錄企業(yè)的回答情況，對非目錄企業(yè)的基礎(chǔ)權(quán)數(shù)進(jìn)行調(diào)整。如果某個(gè)樣本PSU內(nèi)有部分無回答非目錄企業(yè)，則需要在該樣本PSU內(nèi)部，對其他有回答非目錄企業(yè)的權(quán)數(shù)進(jìn)行調(diào)整；如果

30、最終層中只有少數(shù)無回答樣本PSU，則需要在最終層內(nèi)對其他有回答樣本PSU的權(quán)數(shù)進(jìn)行調(diào)整；如果某個(gè)最終層中絕大多數(shù)的樣本PSU無回答，則首先需要進(jìn)行最終層合并，然后在合并層內(nèi)對其他有回答樣本PSU進(jìn)行權(quán)數(shù)調(diào)整；如果沒有無回答非目錄企業(yè)和樣本PSU，則不需要進(jìn)行權(quán)數(shù)調(diào)整。調(diào)整權(quán)數(shù)的具體方法參見規(guī)模以下工業(yè)抽樣調(diào)查方案。由于非目錄企業(yè)的無回答權(quán)數(shù)調(diào)整方法與個(gè)體經(jīng)營工業(yè)單位相同，因此，可以參照個(gè)體經(jīng)營工業(yè)單位子總體的三種無回答情況舉例，將其中的“個(gè)體經(jīng)營工業(yè)單位”替換成“非目錄企業(yè)”來理解非目錄企業(yè)無回答權(quán)數(shù)調(diào)整過程。第四步，用下列命令在基層數(shù)據(jù)中設(shè)置抽樣設(shè)計(jì)信息，然后保存數(shù)據(jù)集。use fmlqy

31、*.dta, clearsvyset strata codesvyset pweight weightsvyset psu dmsvyset fpc total_csave fmlqy *.dta, replace第五步，進(jìn)行總量和方差估計(jì)（包括估計(jì)非目錄企業(yè)數(shù)及其方差）。首先，需要將各調(diào)查指標(biāo)的計(jì)量單位轉(zhuǎn)換成萬元，a01（全部從業(yè)人員期末數(shù)）除外；第二，計(jì)算各指標(biāo)的樣本匯總值；第三，估計(jì)企業(yè)數(shù)及各調(diào)查指標(biāo)的總量和方差。use fmlqy*.dta, clearfor var a*:replace X=X/10replace a01=a01*10egen hz00=sum(dws)egen h

32、z01=sum(a01)egen hz15=sum(a15)svytotal dws a*第六步，制綜合表。參見目錄企業(yè)部分。五、企業(yè)子總體總量和方差推算步驟第一步，統(tǒng)一PSU標(biāo)識變量，縱向合并（追加）目錄企業(yè)和非目錄企業(yè)基層表數(shù)據(jù)集。在進(jìn)行數(shù)據(jù)集追加之前，需要檢查兩個(gè)數(shù)據(jù)集中同名變量的類型是否一致，主要是指“字符型”和“數(shù)值型”的區(qū)別，對于“數(shù)值型”內(nèi)部的細(xì)分類不一致無關(guān)緊要。use fmlqy*.dta, cleargen str12 ppsu=dmsave temp1.dta, replaceuse mlqy00.dta, cleargen str12 ppsu=qyfrdmappend

33、 using temp1.dtasave qy*.dta, replaceerase temp1.dta第二步，在基層數(shù)據(jù)中設(shè)置抽樣設(shè)計(jì)信息。use qy *.dta, clearsvyset strata codesvyset pweight weightsvyset psu ppsusvyset fpc total_csave qy *.dta, replace第三步，進(jìn)行總量和方差估計(jì)（包括估計(jì)非目錄企業(yè)數(shù)及其方差）。首先，需要將各調(diào)查指標(biāo)的計(jì)量單位轉(zhuǎn)換成萬元，a01（全部從業(yè)人員期末數(shù)）除外；第二，計(jì)算各指標(biāo)的樣本匯總值；第三，估計(jì)企業(yè)數(shù)及各調(diào)查指標(biāo)的總量和方差。use qy*.dta

34、, clearfor var a*:replace X=X/10replace a01=a01*10egen hz00=sum(dws)egen hz01=sum(a01)egen hz15=sum(a15)svytotal dws a*第四步，制綜合表參見目錄企業(yè)部分。六、規(guī)模以下工業(yè)總體總量和方差推算步驟第一步，統(tǒng)一PSU標(biāo)識變量，縱向合并（追加）企業(yè)（目錄企業(yè)+非目錄企業(yè)）和個(gè)體經(jīng)營工業(yè)單位基層表數(shù)據(jù)集。在進(jìn)行數(shù)據(jù)集追加之前，需要檢查兩個(gè)數(shù)據(jù)集中同名變量的類型是否一致，主要是指“字符型”和“數(shù)值型”的區(qū)別，對于“數(shù)值型”內(nèi)部的細(xì)分類不一致無關(guān)緊要。use gt*.dta, clearge

35、n str12 ppsu=dmdrop a04save temp1.dta, replaceuse qy*.dta, cleardrop a03 a04 a05 a06 a07 a08 a09 a10 a11 a12 a13 a14 a15append using temp1.dtasave zt*.dta, replaceerase temp1.dta第二步，在基層數(shù)據(jù)中設(shè)置抽樣設(shè)計(jì)信息。use zt *.dta, clearsvyset strata codesvyset pweight weightsvyset psu ppsusvyset fpc total_csave zt *.dt

36、a, replace第三步，進(jìn)行總量和方差估計(jì)（包括估計(jì)非目錄企業(yè)數(shù)及其方差）。首先，需要將各調(diào)查指標(biāo)的計(jì)量單位轉(zhuǎn)換成萬元，a01（全部從業(yè)人員期末數(shù)）除外；第二，計(jì)算各指標(biāo)的樣本匯總值；第三，估計(jì)企業(yè)數(shù)及各調(diào)查指標(biāo)的總量和方差。use zt*.dta, clearfor var a0*:replace X=X/10replace a01=a01*10egen hz00=sum(dws)egen hz01=sum(a01)egen hz02=sum(a02)svytotal dws a0*4、制綜合表參見目錄企業(yè)部分。下面是用于推算規(guī)模以下工業(yè)總體總量和方差的STATA命令程序?qū)嵗?、用于推?/p>

37、總量和方差的STATA命令程序?qū)嵗?、目錄企業(yè)部分clearset more offcd （數(shù)據(jù)集存放路徑）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 nonqys>0 display "基層表數(shù)據(jù)集和樣本數(shù)據(jù)集有部分企業(yè)不對應(yīng)！&

38、quot; list qymc qyfrdm _merge if _merge!=3clearelse drop _mergesave mlqy*.dta, replacegen dws=1replace dws=0 if reason>=2svyset 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

39、)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 sdm hz00 hz01 hz02 hz03 hz04 hz05 hz06 hz07 hz08 hz09 hz10 hz11 hz12 hz13 hz14 hz

40、15gen 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 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 jf

41、c00 jfc01 jfc02 jfc03 jfc04 jfc05 jfc06 jfc07 jfc08 %15.2fformat jfc09 jfc10 jfc11 jfc12 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'+1repla

42、ce 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 zhbml_nb*.dta, rep

43、laceclear2、個(gè)體經(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=0 if reason=1svyset strata codesvyset pweight w

44、eightsvyset 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 str2 sdm="（省代碼）"keep smc s

45、dm 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 jfc02 jfc04 %15.2fformat xd

46、w00 xdw01 xdw02 xdw04 %6.2fkeep in 1forvalues i = 0/3 if i'<3 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'=s

47、qrt(variancej',j') in 1 replace xdw0j'=1.96*jfc0j'/a0j'*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.dtarep

48、lace reason=1 if reason=.for var a*: replace X=0 if X=.drop _mergesave fmlqy*.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)ege

49、n 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 sdm hz00 hz01 hz02 hz03 hz04 hz05 hz06 hz07 hz08 hz09 hz10 hz11 hz1

50、2 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 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 %

51、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 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

52、=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 zhbf

53、ml_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, replaceerase temp1.dtasvyset strata codesvyset pweight weightsvyset psu ppsusvyset fpc total_csave q

54、y *.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 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 %10