R/access.R

In spatstat.univar: One-Dimensional Probability Distribution Support for the 'spatstat' Family

#'
#'    access to internal functions, for debugging
#'
#'   kermom()  should agree with kernel.moment()
#' 
#'  Copyright (c) 2023 Adrian Baddeley, Tilman Davies and Martin Hazelton
#'  GNU Public Licence (>= 2.0)

kermom <- function(m, r, kernel="gaussian", mean = 0,
                   sd=1/kernel.factor(kernel)) {
  kernel <- match.kernel(kernel)
  kerncode <- switch(kernel,
                     gaussian=1,
                     rectangular=2,
                     triangular=3,
                     epanechnikov=4,
                     biweight=5,
                     cosine=6,
                     optcosine=7,
                     0)
  df <- data.frame(r=r, mean=mean, sd=sd)
  nr <- nrow(df)
  z <- .C(SK_kermom,
          nx = as.integer(nr),
          x = as.double(df$r),
          mean = as.double(df$mean),
          sd = as.double(df$sd),
          m = as.integer(m),
          kerncode = as.integer(kerncode),
          y = as.double(numeric(nr)),
          errcode = as.integer(integer(1)),
          PACKAGE="spatstat.univar")
  if(z$errcode != 0)
    switch(z$errcode,
           stop("Error in C: negative length"),
           stop("Error in C: unrecognised kernel"),
           stop("Error in C: value of m is not supported"),
           stop(paste("Error in C: error code", z$errcode)))
  return(z$y)
}