具体该如何来操作呢？

我想再问一下蓝色老师，面板固定效应使用r语言怎么输出常数项？
我看到有一个帖子是这样操作的---p1_fe <- plm（formula,data = pgr, model = "within")
                                             p1_fe <- plm(formula,data = pgr, model = "within",type="dmean")
                                              summary(fixef(p1_fe,type="level"))
                                              summary(fixef(p1_fe,type="dmean"))
然后使用后两行命令结果的系数作差得到的。先不论这种方法的可靠性。
我尝试了一下，虽然每个结果都大致相等（还是有一定差别的），但是这样就没有办法得到一个一致的标准差或t值或p值
我查了很多资料没有明确的解释，不知道蓝色老师怎么看这个问题？

两种都没有问题，关键取决于你的常数项的定义。不同的软件定义是不同的。

请把你帖子里面  跪求   的词语删除

好的   其实我也是求解心切  还望蓝老师理解

谢谢蓝色老师的精彩解答！

蓝色老师，我还有个疑问，stata可以做面板分位数回归分析吗？我查了一些资料都表示只能做截面分位数回归分析而不能直接进行面板分位数回归分析，直到现在我也不确定这个问题，后来我就只能求助于r语言了。就这个问题蓝色老师能帮忙解答一下吗？不胜感激！

没有用过面板分位数
下面的链接看看吧，
http://www.rand.org/economics/resources.html

谢谢蓝色老师！
正在努力学习中！