大家好，我现在用r语言做面板数据的固定效应，个体固定效应和时间固定效应都有，我的代码是这样的
plm(y~x1+x2+x3+trdtime2+trdtime3+trdtime4+trdtime5+trdtime6+trdtime7+trdtime8+trdtime9+trdtime10+trdtime11+trdtime12+trdtime13+trdtime14+trdtime15+trdtime16+trdtime17+trdtime18+trdtime19+trdtime20+trdtime21+trdtime22+trdtime23+trdtime24+trdtime25+trdtime26,data=mydata,index = c("stkcd", "trdtime"),effect="twoways",model = "within")
对应的stata我是这样写的：
xtreg y x1 x2 x3 trdtime2-trdtime26,fe r
我输入这个之后r跑了好久好久，不应该啊，用stata一下子就出结果了的。而且我发现回归结果和我用stata 的回归结果又细微的差别，为什么呢？
是不是我plm用错了。请各位前辈指教。刚开始用r。因为之后有个回归没有现成的stata的语句，看见有r的包就想来试一下，结果今天刚试了plm就发现做出来的结果和stata不一样啊。