Nothing
#
# reach.R
#
# $Revision: 1.10 $ $Date: 2022/11/03 11:08:33 $
#
## The generic 'reach' is now in spatstat.random
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)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.