计算改进的 kppm 模型的方差协方差矩阵

Calculating variance covariance matrix for improved kppm model

当我尝试改进 kppm 对象的强度估计时,如果我在函数 improve.kppm:

中包含参数 vcov = TRUE,则会出现以下错误
Error in improve.kppm(object, type = type, rmax = rmax, dimyx = dimyx,  : 
  object 'gminus1' not found

如果我不包含参数,该函数会运行,但我无法在改进的 kppm 对象上使用 summary() 函数。我收到与上面相同的错误消息。当我使用 vcov().

时也会发生同样的事情

我用来创建我的 kppm 对象的调用是(为清楚起见减少了协变量的数量):

a05 = kppm(a2005nests ~ nest + nest2, cluster = "Thomas", covariates = fitcov(2))

其中 fitcov(2) 是 returns im 对象列表的函数。这可能是问题所在吗?我注意到,如果我在原始 kppm 调用中使用此函数,则 kppm 对象上的一些 spatstat 函数会抛出错误。通常它说的是 Error: Covariates ‘nest’ and ‘nest2’ were not found.

improve.kppm的逻辑流程有一个错误:如果vcov=TRUEtype != "quasi",变量gminus1没有定义。我们会尽快在spatstat的开发版中修复这个问题。

您是否可能 select type="clik1"type="wclik1" 在最初调用 kppm 时?

目前,您应该能够通过在调用 kppmimprove.kppm.

第二个问题,kppm 找不到协变量,似乎是一个范围界定问题,但我无法在这里重现。如果您能提供一个最小的工作示例,将会有所帮助。