[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