# 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