[R] if condition and for cycles

Rolf Turner r.turner at auckland.ac.nz
Tue Nov 27 20:42:03 CET 2007


Output from an R command is printed to the console ***only*** if the  
command
is issued *directly* from the command line.  Otherwise the output is  
invisible.
This includes unassigned output from commands issued ***inside  
functions***.
And, finally, functions include things like ``if'' and for loops.

If you want such output to be visibly displayed you must wrap it in  
an explicit
print statement.

In short you can solve your problems by wrapping your summary  
commands inside
print statements.  E.g. instead of

	summary(fitlinear)

write

	print(summary(fitlinear))


		cheers,

			Rolf Turner

On 28/11/2007, at 7:43 AM, ljantunes wrote:

>
> 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")

	Surely this must be ``message("I'here")'' --- R is case sensitive!!!
	Attention to detail is vital.  Loose lips sink ships! :-)
> }
>
> Thanks in advance,
> Luis Antunes

######################################################################
Attention: 
This e-mail message is privileged and confidential. If you are not the 
intended recipient please delete the message and notify the sender. 
Any views or opinions presented are solely those of the author.

This e-mail has been scanned and cleared by MailMarshal 
www.marshalsoftware.com
######################################################################



More information about the R-help mailing list