计算改进的 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=TRUE
和type != "quasi"
,变量gminus1
没有定义。我们会尽快在spatstat
的开发版中修复这个问题。
您是否可能 select type="clik1"
或 type="wclik1"
在最初调用 kppm
时?
目前,您应该能够通过在调用 kppm
和improve.kppm
.
第二个问题,kppm
找不到协变量,似乎是一个范围界定问题,但我无法在这里重现。如果您能提供一个最小的工作示例,将会有所帮助。
当我尝试改进 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=TRUE
和type != "quasi"
,变量gminus1
没有定义。我们会尽快在spatstat
的开发版中修复这个问题。
您是否可能 select type="clik1"
或 type="wclik1"
在最初调用 kppm
时?
目前,您应该能够通过在调用 kppm
和improve.kppm
.
第二个问题,kppm
找不到协变量,似乎是一个范围界定问题,但我无法在这里重现。如果您能提供一个最小的工作示例,将会有所帮助。