adolf.stips at jrc.ec.europa.eu
2011-Dec-12 09:57 UTC
[R] Confidence intervals of gls function?
Dear gls-experts, while reading and testing some examples of the book "introductionary time series analysis with R", I encountered the following fact which puzzles me. Confidence intervals for global temperature time series (P99) computed from general least squares (GLS) to fit the time series. I repeat the example from the book and get the same results: temp.gls=gls(temp ~ time(temp),correlation=corAR1(0.7)) print(confint(temp.gls)) 2.5 % 97.5 % (Intercept) -39.80571914 -28.49658509 time(temp) 0.01442274 0.02011148 But to my surprise the value of the lag 1 acf does not matter: temp.gls=gls(temp ~ time(temp),correlation=corAR1(0.1)) print(confint(temp.gls)) 2.5 % 97.5 % (Intercept) -39.80571914 -28.49658509 time(temp) 0.01442274 0.02011148 I could set it even to 0.0 or -0.1 and get the same result. Why the confidence interval should be independent of the value of the autocorrelation? Many thanks for any hints in this matter, best regards, Adolf stips ------------------------------------------------ Adolf Stips (adolf.stips@jrc.ec.europa.eu) Global Environment Monitoring unit CEC Joint Research Centre, TP 272 I-21027 Ispra, Italy Tel: +39-0332-789876 Fax: +39-0332-789034 I know that I do not know, but even that not for sure! ------------------------------------------------ "The views expressed are purely those of the writer and may not in any circumstances be regarded as stating an official position of the European Commission." [[alternative HTML version deleted]]