运行双向固定效应回归

Question

我正在使用面板数据集。我的模型应如下所示： y =拦截+ x2 + x3 +年份固定效果+区域固定效果+控件+错误

如何在R中运行此双向固定效果模型？有人可以为我编写代码吗？