[R] if condition and for cycles
ljantunes
ljantunes at hotmail.com
Tue Nov 27 19:43:49 CET 2007
I have two questions about R:
1) Does anyone knows why I get the plots but not the summary results of the
linear fitting when I run this code? When I run it outside the if condition
it works fine.
if (op==4) {
t=1:length(time(SeaIce))
fitlinear=lm(SeaIce~t)
windows()
plot.ts(SeaIce)
lines(ts(fitlinear$fitted.values,start=1979,frequency=12))
summary(fitlinear)
fitlinear=ts(fitlinear$fitted.values,start=1979,frequency=12)
c=cos(2*pi*t/12)
s=sin(2*pi*t/12)
fittotal=lm(SeaIce~t+c+s)
summary(fittotal)
windows()
par(mfrow=c(2,1))
plot(t,SeaIce)
lines(fittotal$fit, col=2)
plot(t,fittotal$resid)
windows()
acf(fittotal$resid)
}
2) Why does't this work? I get the message "I'here" ten times but the AIC
results, none.
for (i in 1:10) {
f_arima = arima(SeaIce_STL_Remain,order=c(i,0,0))
AIC(f_arima)
Message("I´here")
}
Thanks in advance,
Luis Antunes
--
View this message in context: http://www.nabble.com/if-condition-and-for-cycles-tf4882403.html#a13972845
Sent from the R help mailing list archive at Nabble.com.
More information about the R-help
mailing list