我是R和统计数据的新手，并在R上进行模型复制练习。 该模型测试了50个国家的贸易自由化如何影响能源消耗。 E = a0 + a1i + a2t + a3 * LIBit + nit，其中E表示能量消耗，a0是常数，a1i是国家固定效应，a2t是时间固定效应，LIBit是贸易自由化的虚拟变量。 贸易自由化有三个层次，虚拟0用于“lib之前”，虚拟1用于orignial模型中的“lib期间”和“lib之后”。
当我使用lm（）测试仅具有国家固定效果的模型时，它似乎工作正常。但是当我用它来进行国家和时间的固定效果时，会弹出警告“系数：（1由于奇点而未定义）”。 我也试过plm（），但它明显偏离了真实的结果。 当我使用plm（）时，我似乎没有像原始模型那样获得Constant a0的结果。 请帮忙！