This is something that has bothered me for many years about built-in Stata commands for panel data. For small T, large N, if I want an alternative to to xtreg, re to allow for an unrestricted T x T within-unit variance-covariance matrix, I seem to be limited to xtgee. There are some important limitations of xtgee. For one, while one can specify corr(uns) to allow an unrestricted T x T correlation matrix, the variances are assumed constant across time. One can get around this by weighting all variables by the inverse of a time-period specific standard deviation, but this is clunky. Even if one does that, the estimation underlying estimation of GEE has its problems. When I run a simulation with high AR(1) serial correlation (say, rho = 0.75), the number of iterations where xtgee fails can be as high as 20%. If GEE were implemented with a two-step feasible GLS procedure, it would almost always run -- the exception being the extremely rare case where the estimated correlation matrix is not positive definite.
I don't think I've seen a user-written command that implements an unrestricted FGLS estimation for the small T, large N case (and allows for fully robust inference). I'd be happy to learn otherwise. Currently, xtgee is closest but doesn't do enough, and computationally, it's unnecessarily hard.
I don't think I've seen a user-written command that implements an unrestricted FGLS estimation for the small T, large N case (and allows for fully robust inference). I'd be happy to learn otherwise. Currently, xtgee is closest but doesn't do enough, and computationally, it's unnecessarily hard.

Comment