我的面板数据框具有25个变量和2000多个观察值。我正在尝试通过使用plm包进行两种固定效果的回归。但是，我想控制这些年来以及各地区之间的差异。但是，运行回归后会弹出以下错误消息：
if（model）z $ model <-mf中的错误： 参数不能解释为逻辑 另外：警告消息： 在lm.fit（x，y，offset = offset，singular.ok = singular.ok，...）中： 多余的参数“索引”，“效果”将被忽略
lm(hyv_aus ~ danger_days_aus+benefit_days_aus+abovemax_aus+belowmin_aus, data=new_final_data, index=c("district", "year"), model = "within", effect="twoways")