对于面板数据，有理论模型：Yit=αi+θt+β1xit+β2(xit)2+uit , αi是个体固定效应，θt是时间效应，想做回归，R的代码是什么啊？实在是没找到，来此求助。