R/dstrauss.R

In antiphon/rstrauss: Simulates Strauss point process, in 2D and 3D

#' density of a strauss process
#' 
#' Now only for cuboidal window
#' 
#' @param x list with x: coordinate matrix (col-wise), bbox: bounding box (col-wise ranges)
#' @param beta beta parameter
#' @param gamma gamma parameter
#' @param range interaction range
#' @param method how to approximate the constant.
#' @param ... passed to constant approximator.
#' 
#' See \code{\link{approx_strauss_constant}} for the methods
#' 
#' The value is in logarithm scale
#' 
#' @export

dstrauss <- function(x, beta, gamma, range, method="OT", ...) {
  pairs <- sum( sapply(geom(x$x, r=range), length)  )/2
  winV <- prod(apply(x$bbox, 2, diff))
  dim <- ncol(x$bbox)
  # the constant:
  lZ <- approximate_strauss_constant(beta, gamma, range, bbox=x$bbox, method=method, ...)
  #
  nrow(x$x)*log(beta)+pairs*log(gamma) - lZ
}