Displaying 1 result from an estimated 1 matches for "fit_lm".
Did you mean:
film
2011 Jun 08
1
Autocorrelation in R
...ick_fix.htm
# "And now for some regression with autocorrelated errors."
# I've tried to follow the example in Pinheiro & Bates (2004), p.
239-244, with no success.
gc_ts = ts(ex32[66:166,"gc"])
yd_ts = ts(ex32[66:166,"yd"])
library(nlme)
trend = time(gc_ts)
fit_lm = lm(gc_ts ~ trend + yd_ts)
acf(resid(fit_lm))
pacf(resid(fit_lm))
gls_ex32_ar1 = gls(gc_ts ~ trend + yd_ts, correlation = corAR1(form=
~yd_ts),method="ML")
summary(gls_ex32_ar1)
_____________________________________________
Dr. Iuri Gavronski
Assistant Professor
Programa de P?s-Grad...