[学习资料] 非线性似不相关（nlsur）、偏向型技术进步1998-2017的数据和stata的do文件 [推广有奖]

你们很多人都问了这个，首先你们自己看一下这个程序什么意思，很简单的。
就是do文件，program define nlsurab12，这里写的是ab12，这个在创建的do文件里面运行，但是nlsur ad4 @ lnrKY lnwLY lnY Y K L t if dm==11, parameters(b1 b2 b3 b4 b5 b6 b7) initial(b1 0.1 b2 1.37 b3 1 b4 0.0001 b5 1 b6 0.002 b7 1) nequations(3) ifgnls，但是在stata命令行里面运行的是nlsur ad4。这个地方要和前面保持一次，就这里的问题，这么明显，你们看一下程序啊。这个上传的时候，是因为我数据有28个省份，分别做，所以不同省份的ad这个取名不一样，放上去上传的时候没注意。很明显啊这个

点赞支持
感谢分享

有问题，你自己先看一看想一想，我自认为写的很清楚了，确实理解不了，可以直接留言在下面，我看到了会回的

initial(b1 0.1 b2 1.37 b3 1 b4 0.0001 b5 1 b6 0.002 b7 1) 请问这个括号里面的值是怎么确定的呢？需要改吗？已购买，求解答

已购买，也出现这个问题，不知道怎么解决，和前面保持一致是什么意思？能不能给一个能算出和你一样结果的do文件呢，跑一个跑出和您一样结果，才能证明我的操作没问题，不然心里没底。

已经没有那个报错了，还有几个问题请教一下：
1.version=16是什么意思？需要改吗？
2.ab12是对应dm12？
3.initial(b1 0.1 b2 1.37 b3 1 b4 0.0001 b5 1 b6 0.002 b7 1)括号里的是怎么确定的？需要改吗？
4.算出的只是b1-7的值吗？进步偏向指数呢？
5.excel只是均值，您这个命令可以算出每一年的进步偏向指数吗？

这篇论文和下面截图的中文参考论文一定要看。技术进步偏向对产业结构优化的作用机制与实证检验———基于要素技术效率视角——李 爱，盖骁敏。看看了你才会明白为什么要选这样的初始值。version 16代表stata版本不影响。非线性似不相关回归是迭代计算，计算很慢，多的几分钟一个，而且还要调参数，我熟练了算完全部都花了一天半，所以不可能一次性出结果。我是分省份单独计算，所以dm就代表省份，按照统计年鉴上各省的排序从1-28.

技术偏向指数，我用的MATLAB计算的。流程就是参考上面我推荐的李爱的文章里面的公式编写代码计算。
第一步分省份先计算出来那7个待估参数值，第二步带入代码里面，改我程序里面标注的地方。最好你们自己能看懂公式和我程序为什么要那样写，里面有可以调整的地方，自己多看点论文琢磨。不要指望一晚上或者靠别人就能把偏向性技术进步给弄懂，我自己也是花了一个多月时间才最终敲定了这个测算结果。

%计算资本、劳动技术效率和技术进步偏向指数
AAA1=AAA;%t Y K L rKY wLY
[m1,n1]=size(AAA1);
t=10.5;
g=0;
for i=1:20:m1
      lk=i+19;
      PL=AAA1(i:lk,:);
      g=g+1;
      eval(['AS',num2str(g),'=','PL']);
end
A=eval(['AS',num2str(27)]);%27对应的就是dm省份直接改
A1=A;
[m,n]=size(A1);
GK=sum(A1)/m;%求均值


BB1=BB;%(b1 b2 b4 b5 b6 b7)%对应省份计算出来的参数值
aa1=BB1(1,2);
aa=BB(1,2)/(BB(1,2)-1);
bb=(BB(1,2)-1)/BB(1,2);

%资本技术效率
F1=[];
for i=1:m
    a1=A1(i,2)/A1(i,3);
    a2=A1(i,5)/BB1(1,1);
    a3=a2^aa;
    a4=a1*a3;
    F1=[F1;a4];
end
f1=sum(F1)/m

%劳动技术效率
F2=[];
for j=1:m
    qq1=A1(j,2)/A1(j,4);
    qq2=A1(j,6)/(1-BB1(1,1));
    qq3=qq2^aa;
    qq4=qq1*qq3;
    F2=[F2;qq4];
end
f2=sum(F2)/m

%技术进步偏向指数
M=[0];
mk=m-1;
for i=1:mk
    j1=i+1;
ak= F1(j1,1)/F1(i,1);
bk=F2(j1,1)/F2(i,1);
M1=(ak-bk)*bb;
M=[M;M1];
end
BK=[F1,F2,M];
s=sum(M)/mk

如果有用，麻烦评个分、点个赞。确实偏向性技术进步太难了，网上也没有个系统的教程，自己也是弄了很久，希望对你们有用。