楼上的朋友，我看到有文献是这么做的：xtgls ... , igls panels(heteroskedastic)
                                                                        estimates store hetero
                                                                        xtgls ...
                                                                       local df = e(N_g) - 1
                                                                       lrtest hetero . , df(`df')
  我按照以上步骤做，最后stata给出的结果是：does not contain scalar e(ll)，所以就不知道如何办了。你的问题解决了吗？