[R] P-values in gls

[Wilson, Andrew]

When fitting a simple linear or polynomial regression using lm, R
provides a p-value for the whole model as well as for the individual
coefficients.  When fitting the same models using gls (in order to
correct for autocorrelation), there doesn't seem to be a p-value
provided for the whole model, although LL, AIC and BIC statistics are
provided.  Is it possible to obtain a p-value for the whole model?

Many thanks,
Andrew Wilson