在stata中分析面板数据异方差检验，结果输出“does not contain scalar e(ll)”，以下是过程：

方法1：

xtgls y x, panel(heteroskedastic) igls
estimates store igls
xtgls y x
local df=e(N_g)-1
lrtest igls . , df(`df')

真诚请教，哪里出错了！谢谢！