要做计量的式子是:
转换成计量方程是:
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的原始数据如下: