R/bw.scott.R

In spatstat.core: Core Functionality of the 'spatstat' Family

#'
#'   bw.scott.R
#'
#'   Bandwidth selection rule bw.scott for point patterns in any dimension
#'
#'   $Revision: 1.1 $  $Date: 2019/07/22 11:41:41 $

bw.scott <- function(X, isotropic=FALSE, d=NULL) {
  stopifnot(is.ppp(X) || is.lpp(X) || is.pp3(X) || is.ppx(X))
  if(is.null(d)) { d <- spatdim(X, intrinsic=FALSE) } else check.1.integer(d)
  nX <- npoints(X)
  cX <- coords(X, spatial=TRUE, temporal=FALSE, local=FALSE)
  sdX <- apply(cX, 2, sd)
  if(isotropic) {
    #' geometric mean
    sdX <- exp(mean(log(pmax(sdX, .Machine$double.eps))))
  }
  b <- sdX * nX^(-1/(d+4))
  names(b) <- if(isotropic) "sigma" else paste0("sigma.", colnames(cX))
  return(b)
}

bw.scott.iso <- function(X) { bw.scott(X, isotropic=TRUE) }