[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