R plm包中的面板数据

Question

我在R中估计面板数据模型时遇到问题。 我想估计各国实际国内生产总值和相对价格水平的变化对旅游业贡献的影响。

如果我使用命令

Y <- cbind(ln_Differences_Contribution)
X <- cbind(ln_price_differences, Differences_ln_gdp)

然后

fixed <- plm(Y~X, data=pdata, model = "within")

我对不同的年份没有影响。

无论如何我可以添加时间变量吗？

Answer 1

如果要控制时间段（“年”）的效果，可以使用effect参数在模型中双向指定时间效果：

plm(y ~ x1 + x2, data = pdata, model = "within", effect = "twoways")

另一种方法是明确将时间变量包含在公式中作为因子，并在问题中指定单向个体模型，假设您的时间变量被称为year：

plm(y ~ x1 + x2 + factor(year), data = pdata, model = "within")

这些年份的显式估计值会对输出产生一定的“污染”（通常，人们对这些估计不感兴趣）。