为ppp对象生成集群的spatstat标记

Question

这个问题非常接近here的问题。如果我们想为一个已经存在的点模式生成随机标记，答案是很好的-我们从多元正态分布中提取并与每个点关联。

但是，对于我自己的点模式，我需要生成跟随 // Get the first section header. This form item will hold the message. var item = formItems[i]; // Get clocked type var type = e.response.getItemResponses()[j].getResponse(); 随附的lansing数据集中给出的标记的标记。换句话说，我有一个没有标记的点模式，我想用一个确定的模式模拟标记（例如，为了说明我自己的数据的隔离概念）。我该如何做标记？我了解spatstat与我的数据集之间的点数可能会有所不同，但是我可以缩小窗口或创建更多点。谢谢！

Answer 1

这是在四个不同的矩形中隔离的另一个版本 地区。

library(spatstat)

p <- c(.6,.2,.1,.1)
prob <- rbind(p,
              p[c(4,1:3)],
              p[c(3:4,1:2)],
              p[c(2:4,1)])
X <- unmark(spruces)
labels <- factor(LETTERS[1:4])
subwins <- quadrats(X, 2, 2)
Xsplit <- split(X, subwins)
rslt <- NULL
for(i in seq_along(Xsplit)){
  Y <- Xsplit[[i]]
  marks(Y) <- sample(labels, size = npoints(Y),
                     replace = TRUE, prob = prob[i,])
  rslt <- superimpose(rslt, Y)
}
plot(rslt, main = "", cols = 1:4)
plot(subwins, add = TRUE)

Answer 2

隔离是指一个物种在一个物种中占主导地位的事实。 观察窗口的特定部分。一个极端的例子是 完全基于x坐标。这会产生条 不同类型的点：

library(spatstat)

X <- lansing
Y <- cut(X, X$x, breaks = 6, labels = LETTERS[1:6])
plot(Y, cols = 1:6)

在不了解有关所需隔离类型的更多详细信息的情况下 很难提出更有用的建议。