我正在使用R中的plm-package处理面板数据 . 除了国家和年份固定效应之外,我还需要包含线性国家特定时间趋势 . 我目前的代码是:

result <- plm (lngdpwdilocal ~ lndn  + trends, data = global_total_dn_uncal, index = c("isonv10","year"), model = "within", effect = "twoways")

“趋势”是一个年变量,可以计算每个国家的年份 . 到目前为止,我得到了相同的结果,就像我没有包含趋势一样 .

以下Stata代码有效(但是,我需要用plm实现它):

xi: xtreg lngdpwdilocal lndn i.year           i.isonv10|year , fe robust

如果它有帮助我很乐意提供更多数据 . 到目前为止,我发现了两个相关的线不幸的是,他们没有让我解决我的问题:

How do I control for unit-specific trend in a panel regression with plm

Fixed effects regression with state specific trends