est_intens_weight {highriskzone} | R Documentation |
Estimates the intensity of the point pattern by a kernel method
(See density.ppp
).
est_intens_weight(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 |
data(craterA) #change npixel = 50 to 1000 to get a nicer picture spatstat::spatstat.options(npixel=50) # use only ten observations for fast computation thin.craterA <- craterA[1:10] # weight first 5 observations twice weights <- c(rep(2, 5), rep(1, 5)) int <- est_intens_weight(thin.craterA, weights = weights) plot(int$intensest, main = "pixel image of intensity") plot(craterA$window, main = "contour plot of intensity") contour(int$intensest, add =TRUE)