est_intens | R Documentation |
Estimates the intensity of the point pattern by a kernel method
(See density.ppp
).
est_intens(ppdata, covmatrix = NULL, weights = NULL)
ppdata |
data of class ppp |
covmatrix |
(Optional) Covariance matrix of the kernel of a normal distribution |
weights |
(Optional) vector of weights attached to each observation |
A list of
intensest |
Estimated intensity (object of class "im", see |
covmatrix |
Covariance matrix. If |
density.ppp
, Hscv
, eval.im
data(craterA)
#change npixel = 50 to 1000 to get a nicer picture
spatstat.geom::spatstat.options(npixel=50)
# use only ten observations for fast computation
thin.craterA <- craterA[1:10]
int <- est_intens(thin.craterA)
# Plot estimated intensity
plot(int$intensest, main = "pixel image of intensity")
plot(craterA$window, main = "contour plot of intensity")
contour(int$intensest, add =TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.