pggls(formula, data, subset, na.action, effect = c("individual", "time"), model = c("within","random","pooling","fd"), index = NULL, ...)
"summary"(object, ...)
"print"(x,digits = max(3, getOption("digits") - 2), width = getOption("width"),...)
pggls
,data.frame
,lm
,lm
,"individual"
or "time"
,"within"
, "pooling"
, "random"
or "fd"
,pdata.frame
,c("pggls","panelmodel")
containing:
effect="time"
) covariance of errors, pggls
is a function for the estimation of linear panel models
by general feasible generalized least squares, either with or without
fixed effects. General FGLS is based on a two-step estimation process:
first a model is estimated by OLS (pooling
), fixed effects
(within
) or first differences (fd
), then its residuals are
used to estimate an error covariance matrix for use in a feasible-GLS
analysis.
This framework allows the error covariance structure inside every group (if
effect="individual"
, else symmetric) of observations to be fully
unrestricted and is therefore robust against any type of intragroup
heteroskedasticity and serial correlation. Conversely, this structure is
assumed identical across groups and thus general FGLS estimation is
inefficient under groupwise heteroskedasticity. Note also that this
method requires estimation of $T(T+1)/2$ variance parameters, thus
efficiency requires N > > T (if effect="individual"
, else the
opposite). The model="random"
and model="pooling"
arguments both produce an unrestricted FGLS model as in Wooldridge,
Ch. 10, although the former is deprecated and included only for
retrocompatibility reasons.
If model="within"
(the default) then a FEGLS (fixed
effects GLS, see ibid.) is estimated; if model="fd"
a FDGLS
(first-difference GLS).
Im, K. S. and Ahn, S. C. and Schmidt, P. and Wooldridge, J. M. (1999) Efficient Estimation of Panel Data Models with Strictly Exogenous Explanatory Variables, Journal of Econometrics, 93, 177-201.
Wooldridge, J. M. (2002) Econometric Analysis of Cross Section and Panel Data, MIT Press.
Wooldridge, J. M. (2010) Econometric analysis of cross-section and Panel Data, 2nd ed., MIT Press.
data("Produc", package = "plm")
zz <- pggls(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp, data = Produc, model = "pooling")
summary(zz)
Run the code above in your browser using DataLab