R语言 关于使用OLS方法对CES函数进行回归的问题求助

要做计量的式子是：

转换成计量方程是：

y, x1, x2, x3的数据都有，要对该方程进行回归，还需要ρ和γ的数据。假定γ的范围是(0,1)，即x1和x2的替代弹性大于1，呈现出替代关系，份额参数ρ的范围在(0,1)。我想做的是分两步做回归：
（一）定义变量z为

在γ=seq(0.001, 0.999, by=0.001)和ρ=seq(0.001, 0.999, by=0.001)的取值下，求出z在所有(ρ,γ)组合下的值。（共有999*999=998001组这样的z）
（二）对Log y =log A + [(1-θ) /γ] log z+ θ log x3 + u 进行回归。共有998001个回归，然后选出这些回归中拟合优度R2最高的(ρ,γ)组合，变量z斜率参数估计值t值最大的(ρ,γ)组合。

y, x1, x2, x3的原始数据如下：

1. y <- c(-2.996496, -3.010378, -2.429509, -1.740558, -1.838652, -1.805887, -2.137172, -1.860640,-1.927382, -1.970786, -2.088145, -2.112948, -1.669444, -1.117582, -1.157516, -1.481864)
   
2. x1 <- c(0.7522794, 0.7976038, 0.8640433, 0.9832515, 1.0737857, 1.2448860, 1.3956022, 1.4940608, 1.6059855, 1.7718685, 2.0125792, 2.4570253, 2.8883200, 3.4699692, 3.9103071, 4.7074087)
   
3. x2 <- c(1.041848, 1.156326, 1.302007, 1.504695, 1.695061, 1.899581, 2.169386, 2.336143, 2.414350, 2.597997, 2.629217, 2.726548, 2.617122, 2.839649, 2.850892, 2.966809)
   
4. x3 <- c(1.961502, 1.968505, 1.975363, 1.987506, 2.011261, 2.032320, 2.053623, 2.068644, 2.088111, 2.105985, 2.122804, 2.137528, 2.151685, 2.173351, 2.209894, 2.243792)

复制代码

这种计量方程真要好好学学。。
yong79774279629874.jimdo.com，liangliangbai702.jimdo.com，guanxuejiao522334.jimdo.com，suntani57330.jimdo.com，shanheneiyou8.jimdo.com，hui7697801.jimdo.com，ws600801.jimdo.com，luguzonxyce.jimdo.com，hui2095070912.jimdo.com，zhang01806965.jimdo.com，leiz5230.jimdo.com，xianfutj481.jimdo.com，gu67340540099428.jimdo.com，za649107.jimdo.com，zhuo1552223.jimdo.com，zhuomatang2q33.jimdo.com，zhi442573962.jimdo.com，wei050910697.jimdo.com，budit9904.jimdo.com，tongqiong87418.jimdo.com，panchiyin3ow8.jimdo.com，yingyingbat.jimdo.com，qinyou52470.jimdo.com，chengyangshibi.jimdo.com，ixynbnu.jimdo.com，ixynbnu.jimdo.com，chuwu990436.jimdo.com，robin55ym.jimdo.com，zhong3591390.jimdo.com，edrks20.jimdo.com，mezfq62.jimdo.com，lunshai969912037.jimdo.com，orange9036.jimdo.com，souyuane1318.jimdo.com，dunmujia77.jimdo.com，yunguanao460.jimdo.com，chentianzaixian.jimdo.com，ywllx247.jimdo.com，weihanri97417.jimdo.com，jiaomujiaohuan47.jimdo.com，taoyuanxin4551.jimdo.com，huixingc6246.jimdo.com，hanhuanshan403.jimdo.com，buduidui3954.jimdo.com，kuangkongji4042.jimdo.comb28980980renjingy.jimdo.com，qa63597432.jimdo.com，kuding615.jimdo.com，jiaodi8.jimdo.com，183427464.jimdo.com，pomusha01871671.jimdo.com，posiyou43.jimdo.com，fuqiao51069097.jimdo.com，diwozhao6088891054.jimdo.com，womeicu118550631.jimdo.com，womeicu118550631.jimdo.com，chunjiyiba3307.jimdo.com，jiaodi8.jimdo.com，zhiju7151157873.jimdo.com，limeng999999999.jimdo.com，ci33212731180.jimdo.com，12a96b148508.jimdo.com，cisou73715966.jimdo.com，lieyepaolan0.jimdo.com。shishao584.jimdo.comkeliangsong0709290.jimdo.com，

本人已经解决，只能采用数值模拟的方法，设定分配系数和弹性系数的范围，然后进行回归，对于每个（分配系数，弹性系数）组合，都有一个回归结果。

fanghang716 发表于 2015-12-8 23:16
本人已经解决，只能采用数值模拟的方法，设定分配系数和弹性系数的范围，然后进行回归，对于每个（分配系数 ...

lz能不能推荐以下资料，我也正在做，见别人代码中给参数分配数值，就是不知道他其中的原理，希望lz能够帮忙解决一下。

sanyu11 发表于 2017-10-9 11:32
lz能不能推荐以下资料，我也正在做，见别人代码中给参数分配数值，就是不知道他其中的原理，希望lz能够帮 ...

你们好！你们还在做这部分内容么？我也在用CES函数估计个体层面的投入与产出，模型形式和楼主的差不多，但是STATA一直报错R481, can't calculate derivatives。这是什么原因呢？有没有可能是我的数据中含评级变量（5级打分）引起的？

sanyu11 发表于 2017-10-9 11:32
lz能不能推荐以下资料，我也正在做，见别人代码中给参数分配数值，就是不知道他其中的原理，希望lz能够帮 ...

你们好！你们还在做这部分内容么？我也在用CES函数估计个体层面的投入与产出，模型形式和楼主的差不多，但是STATA一直报错R481, can't calculate derivatives。这是什么原因呢？有没有可能是我的数据中含评级变量（5级打分）引起的？