0

我已经标记了 ppp 对象(shp),其中包括:x 坐标,y 坐标,生长(“成人”、“青少年”、“树苗”)。

我使用齐次泊松过程发现了成人和青少年之间关联的模拟信封(使用配对相关函数)。

aj1<-envelope(shp,pcfcross,nsim = 199,i = "adult", j = "juvenile",savefuns=TRUE) 

现在我想使用同质托马斯过程找到相同的东西。所以我将一个均匀的托马斯过程拟合到一个点模式对象上,使用

fit1 <- kppm(shp ~1, "Thomas",method ="palm").

然后我尝试使用

aj2<-envelope(fit1,pcfcross,nsim = 199,i = "adult", j = "juvenile",savefuns=TRUE) 

但这没有用。如果有人有任何建议或任何替代方法来使用同质托马斯过程找到成人和青少年之间关联(使用配对相关函数)的模拟信封,我将非常感激。

4

3 回答 3

0

不幸的是,我认为不存在任何估计多类型泊松聚类模型的通用软件。我们希望spatstat在不久的将来(但实际上是几个月后)添加这样的功能。我知道有些人在这方面工作过,你可能会向他们索要代码。我认为纸

Jalilian, A.、Guan, Y.、Mateu, J. 和 Waagepetersen, R. (2014) Multivariate product-shot-noise Cox models, Biometrics, 71, 1022-1033。

可能与您非常相关,您可以询问作者是否有任何您可以使用的代码。

于 2016-03-13T09:47:07.433 回答
0

来自Spatial Lecture Note,以下示例是否有帮助

require('spatstat')
data("redwood")
X <- redwood
plot(X)
plot(envelope(X))
fit <- kppm(X, ~1, "Thomas")
plot(fit)
plot(envelope(fit))

托马斯信封示例

于 2016-03-13T08:03:01.867 回答
0

目前kppm不处理标记点模式。当您尝试拟合模型时,您会收到有关“无法处理标记点模式”的错误消息kppm(shp ~ 1, "Thomas", method="palm")

此功能将spatstat在 2016 年的某个时候添加。

于 2016-03-14T05:26:53.440 回答