工业抽样推算操作手册

1、规模以下工业抽样调查推算操作手册（年报和半年报）年报和半年报采用的样本是一致的，在数据处理上的唯一区别是在年报调查中，企业部分多“实收资本”等7个调查指标，下面以年报数据为例来阐述用STATA软件进行总量和方差推算的方法和步骤，在半年报中，只需要在企业部分少推算“实收资本”等7个调查指标即可。一、数据转换用Stata/Transfer将基层表Foxpro数据库转换成Stata数据集。二、目录企业总量和方差推算步骤第一步，给基层表数据集赋予基础权数。通过将基层表数据集和样本数据集进行横向合并来对基层表数据集赋予基础权数。在进行两个数据集横向合并之前，必须确保两个数据集中的观察记录是一一对应关系，

2、换句话说就是，通过变量“qyfrdm”能够实现两个数据集之间的完全匹配。两个数据集的观察记录是否一一对应，可以从横向匹配的结果中看出。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 "基层表数

3、据集和样本数据集有部分企业不对应！"list qymc qyfrdm _merge if _merge!=3clearelse drop _mergesave mlqy*.dta, replace说明：上述命令中包含的“*”表示省代码，下同；将赋予了基础权数的基层表数据集保存成mlqy*.dta是为了使n01_nb.dta不发生任何改变，便于数据管理。第二步，新增一个变量并命名为dws，作为有效样本单位和无效样本单位 有效样本单位和无效样本单位的定义参见规模以下工业抽样调查设计。的标识变量，用于估计企业数。变量dws的赋值方法为：有效样本单位赋值为1，无效样本单位赋值为0。use m

4、lqy*.dta, cleargen dws=1replace dws=0 if reason>=2save mlqy*.dta, replace第三步，根据样本企业回答情况，对样本企业的基础权数进行调整。如果没有无回答样本企业，则不需要进行权数调整；如果最终层内只有少数样本企业无回答，则在最终层内部进行权数调整；如果最终层内绝大部分样本企业无回答，则首先需要进行最终层合并，然后在合并层内进行权数调整。权数调整的范围是最终层（或合并层）内的有效样本。调整权数的具体方法参见规模以下工业抽样调查方案。为了保存权数调整过程，需要新增一个变量（baseweight）保存基础权数，将调整后的权数保

5、存在变量weight中。需要注意的是，只能调整基层表数据集中的权数，不能调整样本数据集中的权数，因为本次调查中的无回答企业，在下次调查中仍然要作为样本企业进行调查，并按照在当次调查中是否回答调查表决定是否进行权数调整。（下同）例一（最终层内有少数样本企业无回答情况举例）：某省发现有1个无回答样本企业，其的基础权数为20，层代码为10611，该层中共抽取了5个样本，假设剩下4个样本企业都是有回答的有效样本企业，则调整因子为5/4，4个有回答的样本企业的权数被调整为20*5/4=25，同时，将那个无回答样本企业的权数调整为“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例二（最终层内有多数样本企业无回答情况举例）：某省发现4个无回答企业，都来自于同一层，层代码为10611，其基础权数都为20，该层中共抽取了5个样本企业，此时，需要将该层中剩下的1个有回答的有效样本企业与临近的层合并后，再进行权数调整。比如有一个临近的层代码为10621，该层中共抽取了5个样本企业，基础权数为10，但其中有1个样本企业消亡，属于无效样本企业。将层

7、代码10611替换成10621，然后在合并层10621内对5（=1+4）个有回答的有效样本企业的权数进行调整，调整因子为，原来属于层10611的1个有回答的有效样本企业的权数调整为20*14/6=46.67，原来属于层10621的4个有回答的有效样本企业的权数调整为10*14/6=23.33。同时，将4个无回答样本企业的权数调整为“0”。由于设置fpc要求fpc标识变量（total_c）在同一层内取值相同，因此，需要对合并层内观察记录的fpc标识变量（total_c）进行调整，调整为合并前两个层的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第四步，在基层数据中设置抽样设计信息。use mlqy*.dta, clearsvyset strata codesvyset pweight wei

9、ghtsvyset fpc total_csave mlqy*.dta, replace说明：是否设置“有限总体校正系数（fpc）”只会影响方差估计值的大小。当每个最终层的抽样比都比较小，即每个最终层的权数都比较大时，是否设置“有限总体校正系数（fpc）”对推算结果的影响不大，但是，当有些最终层的抽样比比较大，特别地，当有全选层，即抽样比等于1，权数等于1 时，需要设置“有限总体校正系数（fpc）”，否则会高估方差。因此，统一要求设置“有限总体校正系数（fpc）”。在STATA中，设置“有限总体校正系数（fpc）”有两种方法，一是将最终层的总体单位数（total_c）设置成“有限总体校正系数（

10、fpc）”，二是将最终层的抽样比设置成“有限总体校正系数（fpc）”，两种设置方法得到的结果是一致的。要求统一采用第一种方法。第五步，进行总量和方差估计（包括估计企业数及其方差）。首先，需要将各调查指标的计量单位转换成万元，a01（全部从业人员期末数）除外；第二，计算各指标的样本汇总值；第三，估计企业数及各调查指标的总量和方差。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），将推算结果保存到一个STATA数据集中，然后将此STATA数据集用

12、Stata/Transfer转换成Foxpro数据库文件，导入到规模以下工业抽样调查数据管理系统中，制表输出。方法二：用EXCEL制表。具体方法为：在STATA结果窗口中选中推算结果表，在EDIT菜单中选用“copy table”功能复制推算结果表，然后打开一个EXCEL文件，执行“粘贴”操作，即将STATA的推算结果拷贝到EXCEL中了，最后在EXCEL文件中计算相对误差（STATA推算结果中没有计算相对误差的参数选项），并编辑制表。举例说明将推算结果保存在一个STATA数据集中的过程。首先，按照综合表的变量结构构造一个数据集；然后，将系统自动保存在内存距阵变量e(est)、e(V_db)等

13、中的推算结果，逐个地存放到数据集的指定位置。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注：当用于推算的数据集中某个或某些估计变量含有缺省值（.）时，内部函数e(V_db)只保存方差估计值

19、，是一个1，n距阵；当数据集中所有估计变量都不含有缺省值时，内部函数e(V_db)保存的是n个估计变量的协方差估计值，是一个n，n距阵。因此，统一要求用于推算的数据集中不能包含缺省值。如果发现有缺省值存在，则可以用以下命令将缺省值替换成0。for var a*: replace X=0 if X=.三、个体经营工业单位总量和方差推算步骤第一步，给基层表数据集赋予基础权数。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第二步，新增一个变量并命名为dws，用于估计个体经营工业单位数。变量dws的赋值方法为：对有调查资料的个体经营工业单位赋值为1，对没有个体经营工业单位的样本PSU观察记录赋值为0。use gt*.dta, cleargen dws=1replace dws=0 if reason=1save gt*.dta, replace第三步，根据样本PSU和个体经营工业单位的回答情况，对个体经营工业单位的基础权数进行调整。如果某个样本PSU内有部分无回答个体经营工业单位

21、，则需要在该样本PSU内部，对其他有回答个体经营工业单位的权数进行调整；如果最终层中只有少数无回答样本PSU 如果某个样本PSU内所有的个体经营工业单位都不回答调查表，则称该样本PSU为无回答样本PSU。，则需要在最终层内对其他有回答样本PSU的权数进行调整；如果某个最终层中绝大多数的样本PSU无回答，则首先需要进行最终层合并，然后在合并层内对其他有回答样本PSU进行权数调整；如果没有无回答个体经营工业单位和样本PSU，则不需要进行权数调整。权数调整的范围是最终层（或合并层）内的有效样本。调整权数的具体方法参见规模以下工业抽样调查方案。“0”。use gt*.dta, cleargen bas

22、eweight=weightreplace weight=baseweight*5/4 if dm=""replace weight=0 if code="" & reason 根据规模以下工业抽样调查数据管理系统中的规定，个体经营工业单位和样本PSU没有填报调查表的原因分为三种：1、该样本PSU中没有个体经营工业单位；2、该样本PSU不回答调查表，为无回答样本PSU；3、该个体经营工业单位不回答调查表，为无回答个体经营工业单位。=3save gt*.dta, replace例二（最终层内有少数样本PSU无回答情况举例）：某省发现有1个无回答样本

23、PSU，其的基础权数为20，层代码为111，该层中共抽取了5个样本PSU，假设剩下4个样本PSU都是有效样本PSU，则调整因子为5/4，4个有效样本PSU的权数，即这4个有效样本PSU内所有个体经营工业单位的权数被调整为20*5/4=25，同时，将那个无回答样本PSU的权数调整为“0”。use gt*.dta, cleargen baseweight=weightreplace weight=baseweight*5/4 if code=111replace weight=0 if code=111 & reason=2save gt*.dta, replace例三（最终层内有较多样本

24、PSU无回答情况举例）：某省发现4个无回答样本PSU，都来自于同一层，层代码为111，其基础权数都为20，该层中共抽取了5个样本PSU，此时，需要将该层中剩下的1个有回答样本PSU与临近的层合并后，再对合并层中有回答样本PSU的权数进行调整。比如有一个临近的层代码为112，该层中共抽取了5个样本PSU，基础权数为10，但其中有1个样本PSU没有个体经营工业单位，属于无效样本PSU。首先将层代码111替换成112，然后在合并层112内对5（=1+4）个有回答的有效样本PSU的权数进行调整，调整因子为，原来属于层111的1个样本PSU的权数调整为20*14/6=46.67，原来属于层112的4个样

25、本PSU的权数调整为10*14/6=23.33。同时，将4个无回答样本PSU的权数调整为“0”。由于设置fpc要求fpc标识变量（total_c）在同一层内取值相同，因此，需要对合并层内观察记录的fpc标识变量（total_c）进行调整，调整为合并前两个层的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第四步，在基层数据中设置抽样设计信息。use gt*.dta, clearsvyset strata codesvyset pweight weightsvyset psu dmsvyset fpc total_csave gt*.dta, replace第五步，进行总量和方差估计（包括估计个体经营工业单位数及其方差）。首先，需要将各调查指标的计量单位转换成万元，a01（全部从业人员期末数）除外；第二，计算各指标的样本

27、汇总值；第三，估计单位数及各调查指标的总量和方差。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*第六步，制综合表参见目录企业部分。四、非目录企业总量和方差推算步骤第一步，给基层表数据集赋予基础权数。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第二步，新增一个变量并命名为dws，用于估计非目录企业数。变量dws的赋值方法为：对有调查资料的非目录企业赋值为1，对没有非目录企业的样本PSU观察记录赋值为0。use fmlqy*.dta, cleargen dws=1replace dws=0 if reas

29、on 根据规模以下工业抽样调查数据管理系统中的规定，非目录企业和样本PSU没有填报调查表的原因分为三种：1、该样本PSU中没有非目录企业；或者，该非目录企业变成无效样本企业，即撤消、被合并、停产且找不到人、转产成为非工业企业、已经纳入规模以上工业统计等； 2、该样本PSU不回答调查表，为无回答样本PSU；3、该非目录企业不回答调查表，为无回答非目录企业。=1save fmlqy *.dta, replace第三步，根据样本PSU和非目录企业的回答情况，对非目录企业的基础权数进行调整。如果某个样本PSU内有部分无回答非目录企业，则需要在该样本PSU内部，对其他有回答非目录企业的权数进行调整；如果

30、最终层中只有少数无回答样本PSU，则需要在最终层内对其他有回答样本PSU的权数进行调整；如果某个最终层中绝大多数的样本PSU无回答，则首先需要进行最终层合并，然后在合并层内对其他有回答样本PSU进行权数调整；如果没有无回答非目录企业和样本PSU，则不需要进行权数调整。调整权数的具体方法参见规模以下工业抽样调查方案。由于非目录企业的无回答权数调整方法与个体经营工业单位相同，因此，可以参照个体经营工业单位子总体的三种无回答情况举例，将其中的“个体经营工业单位”替换成“非目录企业”来理解非目录企业无回答权数调整过程。第四步，用下列命令在基层数据中设置抽样设计信息，然后保存数据集。use fmlqy

31、*.dta, clearsvyset strata codesvyset pweight weightsvyset psu dmsvyset fpc total_csave fmlqy *.dta, replace第五步，进行总量和方差估计（包括估计非目录企业数及其方差）。首先，需要将各调查指标的计量单位转换成万元，a01（全部从业人员期末数）除外；第二，计算各指标的样本汇总值；第三，估计企业数及各调查指标的总量和方差。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*第六步，制综合表。参见目录企业部分。五、企业子总体总量和方差推算步骤第一步，统一PSU标识变量，纵向合并（追加）目录企业和非目录企业基层表数据集。在进行数据集追加之前，需要检查两个数据集中同名变量的类型是否一致，主要是指“字符型”和“数值型”的区别，对于“数值型”内部的细分类不一致无关紧要。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第二步，在基层数据中设置抽样设计信息。use qy *.dta, clearsvyset strata codesvyset pweight weightsvyset psu ppsusvyset fpc total_csave qy *.dta, replace第三步，进行总量和方差估计（包括估计非目录企业数及其方差）。首先，需要将各调查指标的计量单位转换成万元，a01（全部从业人员期末数）除外；第二，计算各指标的样本汇总值；第三，估计企业数及各调查指标的总量和方差。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*第四步，制综合表参见目录企业部分。六、规模以下工业总体总量和方差推算步骤第一步，统一PSU标识变量，纵向合并（追加）企业（目录企业+非目录企业）和个体经营工业单位基层表数据集。在进行数据集追加之前，需要检查两个数据集中同名变量的类型是否一致，主要是指“字符型”和“数值型”的区别，对于“数值型”内部的细分类不一致无关紧要。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第二步，在基层数据中设置抽样设计信息。use zt *.dta, clearsvyset strata codesvyset pweight weightsvyset psu ppsusvyset fpc total_csave zt *.dt

36、a, replace第三步，进行总量和方差估计（包括估计非目录企业数及其方差）。首先，需要将各调查指标的计量单位转换成万元，a01（全部从业人员期末数）除外；第二，计算各指标的样本汇总值；第三，估计企业数及各调查指标的总量和方差。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、制综合表参见目录企业部分。下面是用于推算规模以下工业总体总量和方差的STATA命令程序实例七、用于推算

37、总量和方差的STATA命令程序实例1、目录企业部分clearset more offcd （数据集存放路径）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 "基层表数据集和样本数据集有部分企业不对应！&

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、个体经营工业单位子总体clearset more offcd （数据集存放路径）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、非目录企业clearset more offcd （数据集存放路径）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、工业企业子总体clearset more offcd （数据集存放路径）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