op <- spatstat.options(npixel=32)
X <- rpoispp(function(x,y) { exp(3+3*x) })
fit <- ppm(X ~x+y)
fit <- ppm(X ~x+y, nd=16)
plot(dffit(fit))
plot(dffit(fit, collapse=TRUE))
spatstat.options(op)
Run the code above in your browser using DataLab