Very short introduction (these slides).
Live tutorial with examples of analysis.
Questions and discussion.
2 September 2021
Very short introduction (these slides).
Live tutorial with examples of analysis.
Questions and discussion.
spatstat.xxxx
, which the now almost empty package spatstat
then Depends
on.?spatstat
to get an overview and find hidden gems.You can use spatstat
to describe/summarise any given point set with things like
However, spatstat
really focuses on statistical inference for phenomena that generate random locations (point processes).
library(spatstat) set.seed(42) # Reproducibility Xpois <- rpoispp(100, nsim = 3) plot(Xpois, main = "")
Xhc <- rHardcore(beta = 100, R = .05, nsim = 3) plot(Xhc, main = "")
Xthomas <- rThomas(kappa = 5, mu = 20, scale = .1, nsim = 3) plot(Xthomas, main = "")
lambda <- function(x,y){200*(x^2+y^2)} Xinhom <- rpoispp(lambda, nsim = 3) plot(Xinhom, main = "")
Intensity is a first moment property.
Interaction is a higher moment property (inter-point correlation).
They are confounded and without further assumptions it is impossible to separate them in general.
Often a approach like in time series is used:
Explore ?spatstat
which includes lists of commonly (and less commonly) used functions.
Get the book. Unfortunately we don’t have a license to share an online version as many authors have nowadays. Maybe this will change with a second edition. There are three free sample chapters at https://book.spatstat.org/
Ask questions on stackoverflow under the spatstat tag
Report bugs or make feature requests on GitHub. (If possible find the right sub package repo to put the issue under.)