[R] Vectorizing loop

Liaw, Andy andy_liaw at merck.com
Tue Oct 11 14:54:42 CEST 2005


I don't have anything specific to say.  Only the following suggestion:

Try and find out where in the code most of the time is spent.  R has nice
tools for that.  See ?Rprof.  If it's runifpoint() and Kest() that are
taking most of the computing time, you may not be able to do much better.

Andy

> From: Rainer M. Krug
> 
> Sorry
> 
> runifpoint() and Kest are from the package spatstat
> 
> Rainer
> 
> Liaw, Andy wrote:
> > Not unless we know what runifpoint() and Kest() are.  AFAIK 
> these are not
> > part of base R.  If you use functions from add-on packages, 
> please state
> > them so as not to leave others guessing.  (This is in the 
> Posting Guide,
> > which you were asked to read.)
> > 
> > Andy
> > 
> > 
> >>From: Rainer M. Krug
> >>
> >>Hi
> >>
> >>I have the following loop and would like to vectorize it. Any 
> >>ideas if 
> >>it is possible?
> >>
> >>Thanks,
> >>
> >>Rainer
> >>
> >>Tha Loop:
> >>
> >>for (i in 2:Result$NoSims)
> >>{
> >>	ppp <- runifpoint(Result$NoPlants)
> >>	K <- Kest(ppp)
> >>	Result$LSim[i,] <- sqrt(K$iso / pi) - K$r
> >>	CM <- (Result$LSim[i,] * Result$LSim[i,]) / abs(K$r[2] - K$r[1])
> >>	Result$SigCM[i] <- sum(CM, na.rm=TRUE)
> >>	print(i)
> >>	flush.console()
> >>}
> 
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide! 
> http://www.R-project.org/posting-guide.html
> 
>




More information about the R-help mailing list