[R] a small bug in spatstat::rmh

Evgueni Parilov parilov at cat.nyu.edu
Wed Jul 7 21:30:06 CEST 2004

Time to time, rmh.default fails to simulate a "lookup"-type process on a 

if(all.equal(diff(r),rep(deltar,nlook-1))) {
        equisp <- 1
        par <- c(beta,nlook,equisp,deltar,rmax,h)
    } else {
        equisp <- 0
        par <- c(beta,nlook,equisp,deltar,rmax,h,r)

According to the manual, all.equal should not be used in if-statement 
directly. This works:
identical(all.equal(diff(r), rep(deltar, nlook - 1)),TRUE)


More information about the R-help mailing list