# NOT RUN {
ppm(nztrees ~1, Poisson())
# fit the stationary Poisson process to 'nztrees'
# no edge correction needed
lon <- longleaf
# }
# NOT RUN {
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
# }
# NOT RUN {
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())
# }
# NOT RUN {
# fit nonstationary marked Poisson process
# with different log-cubic trend for each species
# }
Run the code above in your browser using DataLab