[R] if condition and for cycles

[Frede Aakmann Tøgersen]

Compare the evaluation of

for (i in 1:10) i

to

for (i in 1:10) print(i)

See also http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-is-the-output-not-printed-when-I-source_0028_0029-a-file_003f


Best regards

Frede Aakmann Tøgersen
Scientist


UNIVERSITY OF AARHUS
Faculty of Agricultural Sciences
Dept. of Genetics and Biotechnology
Blichers Allé 20, P.O. BOX 50
DK-8830 Tjele

Phone:   +45 8999 1900
Direct:  +45 8999 1878

E-mail:  FredeA.Togersen at agrsci.dk
Web:	   http://www.agrsci.org				

This email may contain information that is confidential.
Any use or publication of this email without written permission from Faculty of Agricultural Sciences is not allowed.
If you are not the intended recipient, please notify Faculty of Agricultural Sciences immediately and delete this email.

 

 

> -----Oprindelig meddelelse-----
> Fra: r-help-bounces at r-project.org 
> [mailto:r-help-bounces at r-project.org] På vegne af ljantunes
> Sendt: 27. november 2007 19:44
> Til: r-help at r-project.org
> Emne: [R] if condition and for cycles
> 
> 
> 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.ht
ml#a13972845
> Sent from the R help mailing list archive at Nabble.com.
> 
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide 
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>