0

我目前正在尝试使用 Henningsen/Henningsen (2011) 的 R 中的 micEconCES 包进行一些估计。我的问题是我对 R 不是很熟悉，我正在尝试实现我自己的数据集来获得包的估计。他们的论文作者为估计创建了这个数据集。

R> set.seed( 123 ) 
R> cesData <- data.frame(x1 = rchisq(200, 10), x2 = rchisq(200, 10), x3 = rchisq(200, 10), x4 = rchisq(200, 10) ) 
R> cesData$y2 <- cesCalc( xNames = c( "x1", "x2" ), data = cesData, + coef = c( gamma = 1, delta = 0.6, rho = 0.5, nu = 1.1 ) ) 
R> cesData$y2 <- cesData$y2 + 2.5 * rnorm( 200 )
R> cesData$y3 <- cesCalc(xNames = c("x1", "x2", "x3"), data = cesData, coef = c( gamma = 1, delta_1 = 0.7, delta = 0.6, rho_1 = 0.3, rho = 0.5, + nu = 1.1), nested = TRUE ) 
R> cesData$y3 <- cesData$y3 + 1.5 * rnorm(200) 
R> cesData$y4 <- cesCalc(xNames = c("x1", "x2", "x3", "x4"), data = cesData, coef = c(gamma = 1, delta_1 = 0.7, delta_2 = 0.6, delta = 0.5, rho_1 = 0.3, rho_2 = 0.4, rho = 0.5, nu = 1.1), nested = TRUE ) 
R> cesData$y4 <- cesData$y4 + 1.5 * rnorm(200)

第一行设置随机数生成器的“种子”，以便可以使用完全相同的数据集复制这些示例。第二行创建一个包含四个输入变量（称为 x1、x2、x3 和 x4）的数据集，每个变量有 200 个观测值，并且是从具有 10 个自由度的随机 χ2 分布生成的。第三、第五和第七个命令使用 micEconCES 包中包含的函数 cesCalc 计算具有两个、三个和四个输入（分别称为 y2、y3 和 y4）的 CES 函数的确定性输出变量给定 CES 生产函数。现在在我的论文中，我试图在资本和劳动力的两个输入案例的总体水平上估计美国的 CES 函数。所以我所做的是我从世界银行数据库中收集了 1990-2015 年的数据，

作者通过以下方式进行了非线性估计

R> cesNls <- nls( y2 ~ gamma * ( delta * x1^(-rho) + (1 - delta) * x2^(-rho) )^(-phi / rho), + data = cesData, start = c( gamma = 0.5, delta = 0.5, rho = 0.25, phi = 1 ) ) R> print( cesNls )

现在我想为我自己的数据集使用完全相同的东西，称为 Data_Extract_From_World_Development_Indicators。所以我首先做的是

ceslan <- cesCalc( xNames = c( "GrossFixedCapitalFormation", "LaborForce" ), data = Data_Extract_From_World_Development_Indicators, coef = c( gamma = 1, delta = 0.6, rho = 0.5, nu = 1.1 ) ) 所以我复制了

R> cesData$y2 <- cesCalc( xNames = c( "x1", "x2" ), data = cesData, coef = c( gamma = 1, delta = 0.6, rho = 0.5, nu = 1.1 ) )

我所做的只是更改数据集的名称，并将 x1 和 x2 替换为资本和劳动力的两个变量。

之后我尝试进行非线性估计

cesulan <- nls(y2 ~ gamma * (delta * GrossFixedCapitalFormation^(-rho) + (1-delta)*LaborForce^(-rho))^(-phi / rho), data = Data_Extract_From_World_Development_Indicators, start = c(gamma = 0.5, delta = 0.5, rho = 0.25, phi = 1) )

现在这就是我的问题所在：我不知道我的数据集中的 y2 变量是什么。我可以在公式中看到 y2 ~ gamma *... 所以它是针对该术语的其余部分绘制的，但我不知道我需要在那里插入什么样的值。有人有建议吗？

提前致谢

4

1 回答 1

0

在Hennigsen &Hennisgen (2011)中，变量y2是使用函数创建的cesCalc。扰动是为了从理论上测试引入的函数cesEst。这个变量应该是你函数的输出（通常是国内生产总值，但不完全是）。

这必须是数字非负值的系列（数据框的 R 列），大小等于您的其他解释变量x1和x2.

于 2018-02-28T10:41:29.567 回答