https://github.com/cran/spatstat
Tip revision: 2af0235a57146624ae8d5b148d47d8d0ab351213 authored by Adrian Baddeley on 14 May 2007, 05:36:09 UTC
version 1.11-5
version 1.11-5
Tip revision: 2af0235
reach.R
#
# reach.R
#
# $Revision: 1.7 $ $Date: 2007/01/10 05:55:41 $
#
reach <- function(x, ...) {
UseMethod("reach")
}
reach.interact <- function(x, ...) {
verifyclass(x, "interact")
irange <- x$irange
if(is.null(irange))
return(Inf)
if(!is.function(irange))
stop("Internal error - x$irange is not a function")
ir <- irange(x)
if(is.na(ir))
ir <- Inf
return(ir)
}
reach.ppm <- function(x, ..., epsilon=0) {
verifyclass(x, "ppm")
# Poisson case
if(is.poisson.ppm(x))
return(0)
# extract info
inte <- x$interaction
coeffs <- coef(x)
if(newstyle.coeff.handling(inte)) {
# extract only interaction coefficients
Vnames <- x$internal$Vnames
coeffs <- coeffs[Vnames]
}
# apply 'irange' function
irange <- inte$irange
if(is.null(irange))
return(Inf)
ir <- irange(inte, coeffs, epsilon=epsilon)
if(is.na(ir))
ir <- Inf
return(ir)
}