如何使用plm控制面板回归中的单位特定趋势

Question

我正在使用R中的plmpackage处理面板数据，并且需要指定具有线性状态特定趋势的状态和年份固定效果模型。到目前为止，我有以下代码：

pnl_sb <- pdata.frame(seatbelt.df, index = c("fips","year"))
twyreg <- plm(fatalityrate~sb_useage+primary*sb_useage+secondary*sb_useage+speed65+speed70+ba08+drinkage21+log(income)+age, data = pnl_sb, model = "within", effect = "twoways")

fips是单位变量状态，year是时间变量。我根本不知道如何实现特定于州的趋势。我认为这将假设使用2路固定效果面板数据，将年份作为时间趋势。但是，我有些疑惑。

Answer 1

pnl_sb <- pdata.frame(seatbelt.df, index = c("fips","year"))
twyreg <- plm(fatalityrate~factor(state)*year, data = pnl_sb, model = "within", effect = "twoways")

具体来说：factor(state)确保将使用所有状态。通过将factor(state)乘以year，可以为每个州创建趋势。