固定效应回归 - 特定国家的时间趋势

我使用 plm 包中的 plm 对具有国家固定效应的面板数据进行了回归。 Breusch-异教时间测试固定效应表明我的回归不需要年份固定效应。我的回归代码如下：

fe <- plm(
  unemployment ~ fdi + investment + inflation + gdpgrowth + 
    school + Trade + naturalresources,data = westa,model = "within"
)

但是，我想控制“特定国家/地区的时间趋势”。不幸的是，我是 R 的新手，我不知道特定国家/地区的时间趋势所需的代码。我可以使用哪个代码？