Questions tagged [spatstat]

spatstat is an R package for analysing spatial data, mainly Spatial Point Patterns.

spatstat is an R package for analysing spatial data, mainly Spatial Point Patterns. It also supports three-dimensional point patterns, and space-time point patterns in any number of dimensions. The package contains over 1000 functions for plotting spatial data, exploratory data analysis, model-fitting, simulation, spatial sampling, model diagnostics, and formal inference. Data types include point patterns, line segment patterns, spatial windows, pixel images and tessellations.

Resources:

441 questions
0
votes
1 answer

Anova for multiple point patterns not working for Strauss model

I just started getting into spatial analysis and am fitting some models to my data. My main goal is to test for spatial regularity (whether there is inhibition between points). I created my hyperframe for the data below. There are 6 point patterns…
0
votes
1 answer

Plotting the K function in spatstat seems to only consider single data (and not pooled data)

I have obtained an mppm object by fitting a model on several independent datasets using the mppm function from the R package spatstat. I then want to compare the K function for the fitted model with the K function for the observations. However, it…
0
votes
1 answer

How to optimize the r parameter when fitting a single model to several datasets in spatstat?

I would like to fit a single model to several independent datasets in R using the spatstat package. Here, I have 3 independent datasets (ppp objects called NMJ1, NMJ2, and NMJ3), to which I want to fit a common model. The way to go should be to use…
0
votes
1 answer

Is it possible to construct "MultiStrauss" or something like "MultiFiksel" using "Pairwise" function in R package spatstat?

As the title, I am trying to use some multitype processes that can be both regulated and clustered. Thus I'm thinking of using Pairwise to construct interactions with marks but seems the function does not accept marks as input. Is this even…
0
votes
1 answer

How do I create a combined colour plot with legend overlaying two marks in Spatstat?

I wanted to create a combined colour plot of two marks in Spatstat with a legend to show the species as well as the diameters of multiple species in one point process pattern. I started with this plot: set.seed(42) species <- LETTERS[1:16] …
0
votes
1 answer

spatstats: Is there a way to incorporate the GPS accuracy error into the analysis?

I am working on a point-pattern database. I created the spatial analysis using the Linhom function of the Spatstats package on R and it worked perfectly. However, the points have a GPS accuracy error associated. So I would like to plot this error in…
0
votes
2 answers

ppp objects not working with deldir function

I am attempting to use deldir with a ppp object, but receive the following error: Error: $ operator is invalid for atomic vectors This does not work: library(deldir) library(spatstat) points <- ppp(x=c(-77.308703, -77.256582, -77.290600, …
Bryan
  • 1,771
  • 4
  • 17
  • 30
0
votes
1 answer

How to modify the value of nrank when plotting the envelope for an mppm object in spatstat?

My object model is an mppm object obtained by fitting a model on several independent datasets using the mppm function from the R package spatstat. Since the envelope function does not accept inputs of type mppm, I used the following code (see…
0
votes
2 answers

How to interpret result of kppm modelling from spatstat

I am recently using spatstat to do a bit of spatial data analysis. I am now really struggling to interpret the result of kppm models. For example, when I fit a log-gaussian cox model, what does the parameter trend, alpha, sigma2, mu and theta mean?…
L2m
  • 1
  • 1
0
votes
1 answer

Some issues Related to simulation from rLGCP and Point Process Models fitting

I trying two generate points from rLGCP function. I assumed that the presence of these points in the Window is governed by two covaiates ras1 and ras2. Hence I need to comptute log-lambda. rm(list=…
Luc
  • 1
  • 2
0
votes
1 answer

Is there a way to create a list of ppp objects using spatstat in R?

I am trying to create a list of landscape objects, which have been grouped according to a variable called sim. Ideally I want a list with landscape objects, that can be called i.e. [[1]] Planar point pattern: 100 points window: rectangle = [0,…
OpenSauce
  • 354
  • 2
  • 14
0
votes
1 answer

How to stratify a raster object in R

My goal in the current work is to sumulate a database containning n variables : (1)latitude of the location, (2)longitude of the location, (3...n-1) covariates, (n)stratum. In stratum variable, I'm looking for the number of the compartment from…
Luc
  • 1
  • 2
0
votes
1 answer

Monte Carlo test of spatial segregation

I have trouble intrepreting the result which I get from the segregation.test method in spatstat. However, I have three different point patterns A,B,C and I want to prove that C and B are correlating whereas A and B not. You can see the Kernel…
Hashriama
  • 173
  • 11
0
votes
1 answer

Plot an envelope for an mppm object in spatstat

My question is closely related to this previous one: Simulation-based hypothesis testing on spatial point pattern hyperframes using "envelope" function in spatstat I have obtained an mppm object by fitting a model on several independent datasets…
0
votes
1 answer

How to adress error message related to S4 object in R?

I'm with spatstat package in R.I'm working on species distribution models mainly on point process models. My goal is to evaluate the predictive performance of point process models of Gibbs, Log-Gaussian Cox, and many more, point processes. I've…
Luc
  • 1
  • 2