R/quadrattestlpp.R

In spatstat.linnet: Linear Networks Functionality of the 'spatstat' Family

#'
#'  quadrattestlpp.R
#'
#'  Quadrat counting tests for lpp and lppm
#'
#'  $Revision: 1.3 $ $Date: 2026/01/21 06:26:39 $

quadrat.test.lpp <-
  function(X, ...,
           tess=NULL,
           nx=5, ny=nx,
           xbreaks=NULL, ybreaks=NULL,
           alternative = c("two.sided", "regular", "clustered"),
           method = c("Chisq", "MonteCarlo"),
           conditional=TRUE, CR=1,
           lambda=NULL, df.est=NULL,
           nsim=1999)
{
   Xname <- short.deparse(substitute(X))
   method <- match.arg(method)
   alternative <- match.arg(alternative)
   tessinfo <- if(!is.null(tess)) {
                 list(tess=tess)
               } else if(!missing(xbreaks)) {
                 list(xbreaks=xbreaks, ybreaks=ybreaks)
               } else if(!missing(nx)) {
                 list(nx=nx, ny=ny)
               } else list()
   do.call(quadrat.testEngine,
           resolve.defaults(list(quote(X),
                                alternative=alternative,
                                method=method,
                                conditional=conditional,
                                CR=CR,
                                fit=lambda,
                                df.est=df.est,
                                nsim=nsim),
                            tessinfo,
                            list(...), 
                            list(Xname=Xname, fitname="CSR on a network")))
}

quadrat.test.lppm <- 
  function(X, ...,
           tess=NULL,
           nx=5, ny=nx,
           xbreaks=NULL, ybreaks=NULL,
           alternative = c("two.sided", "regular", "clustered"),      
           method=c("Chisq", "MonteCarlo"),
           conditional=TRUE, CR=1, df.est=NULL, 
           nsim=1999)
{
   fitname <- short.deparse(substitute(X))
   dataname <- paste("data from", fitname)
   method <- match.arg(method)
   alternative <- match.arg(alternative)
   if(!is.poisson(X))
    stop("Test is only defined for Poisson point process models")
   if(is.marked(X))
    stop("Sorry, not yet implemented for marked point process models")
   Xdata <- response(X)
   dont.complain.about(Xdata)
   tessinfo <- if(!is.null(tess)) {
                 list(tess=tess)
               } else if(!missing(xbreaks)) {
                 list(xbreaks=xbreaks, ybreaks=ybreaks)
               } else if(!missing(nx)) {
                 list(nx=nx, ny=ny)
               } else list()
   do.call(quadrat.testEngine,
          resolve.defaults(list(quote(Xdata), 
                                alternative=alternative,
                                method=method,
                                conditional=conditional, CR=CR,
                                nsim=nsim, 
                                fit=X,
                                df.est=df.est),
                           tessinfo,
                           list(...),
                           list(Xname=dataname, fitname=fitname)))
}

quadrat.test.linearquadratcount <-
  function(X, ...,
           alternative = c("two.sided", "regular", "clustered"),
           method=c("Chisq", "MonteCarlo"),
           conditional=TRUE, CR=1,
           lambda=NULL, df.est=NULL,
           nsim=1999) {
   trap.extra.arguments(...)
   method <- match.arg(method)
   alternative <- match.arg(alternative)
   quadrat.testEngine(Xcount=X,
                      alternative=alternative,
                      fit=lambda, df.est=df.est,
                      method=method, conditional=conditional, CR=CR, nsim=nsim)
}