Poisson: Poisson Point Process Model

Description

Creates an instance of the Poisson point process model which can then be fitted to point pattern data.

Usage

Poisson()

Arguments

Value

An object of class "interact" describing the interpoint interaction structure of the Poisson point process (namely, there are no interactions).

Details

The function ppm, which fits point process models to point pattern data, requires an argument interaction of class "interact" describing the interpoint interaction structure of the model to be fitted. The appropriate description of the Poisson process is provided by the value of the function Poisson.

This works for all types of Poisson processes including multitype and nonstationary Poisson processes.

Examples

Run this code

data(nztrees)
 ppm(nztrees, ~1, Poisson())
 # fit the stationary Poisson process to 'nztrees'
 # no edge correction needed

 data(longleaf)
 <testonly>longleaf <- longleaf[seq(1, longleaf$n, by=50)]</testonly>
 longadult <- longleaf[longleaf$marks >= 30, ]
 longadult <- unmark(longadult)
 ppm(longadult, ~ x, Poisson())
 # fit the nonstationary Poisson process 
 # with intensity lambda(x,y) = exp( a + bx)

 data(lansing)
 # trees marked by species
 <testonly>lansing <- lansing[seq(1,lansing$n, by=30), ]</testonly>
 ppm(lansing, ~ marks, Poisson())
 # fit stationary marked Poisson process
 # with different intensity for each species

ppm(lansing, ~ marks * polynom(x,y,3), Poisson())
 # fit nonstationary marked Poisson process
 # with different log-cubic trend for each species
<testonly># equivalent functionality - smaller dataset
   ppm(amacrine, ~ marks * polynom(x,y,2), Poisson())</testonly>

Run the code above in your browser using DataLab

Data Engineering and BI courses are free this week!