ppm(nztrees ~1, Poisson())
# fit the stationary Poisson process to 'nztrees'
# no edge correction needed
lon <- longleaf
longadult <- unmark(subset(lon, marks >= 30))
ppm(longadult ~ x, Poisson())
# fit the nonstationary Poisson process
# with intensity lambda(x,y) = exp( a + bx)
# trees marked by species
lans <- lansing
ppm(lans ~ marks, Poisson())
# fit stationary marked Poisson process
# with different intensity for each species
## Not run:
# ppm(lansing ~ marks * polynom(x,y,3), Poisson())
# ## End(Not run)
# fit nonstationary marked Poisson process
# with different log-cubic trend for each species
Run the code above in your browser using DataLab