as.spatstat: Function to covert a spatstat quadrature and spatstat data...

In skiptoniam/qrbp: Quasirandom quardature scheme using dirichlet tesselation for spatial Poisson point process models

Function to covert a spatstat quadrature and spatstat data (covariates)

Description

This function will convert the ppmData to a spatstat object that can be used in the spatstat library. This currenly only works for unmarked point process

Usage

as.spatstat(object, ...)

Arguments

object	A ppmData object
\dots	Ignored

Examples

## Not run: 
library(ppmData)
path <- system.file("extdata", package = "ppmData")
lst <- list.files(path=path,pattern='*.tif',full.names = TRUE)
preds <- rast(lst)
window <- preds[[1]]
presences <- subset(snails,SpeciesID %in% "Tasmaphena sinclairi")
object <- ppmData(npoints = 10000, presences=presences, window = window, covariates = preds)
ss.object <- as.spatstat(object)

## now you should be able to fit a ppm from spatstat
ss.form <- ss.object$Q~1+poly(min_temp_coldest_month,2)+poly(annual_mean_precip,2)
ppmfit <- ppm(ss.form,data=ss.object$covariates,Poisson())
pred <- predict(ppmfit, covariates = ss.object$covariates,type="trend")

## End(Not run)

skiptoniam/qrbp documentation built on May 13, 2023, 2:08 a.m.