当我尝试提高kppm对象的强度估计时,如果我在函数vcov = TRUE
中包含参数improve.kppm
,则会出现以下错误:
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)
是一个返回im
个对象列表的函数。这可能是问题吗?我注意到,如果我在原始kppm
调用中使用此函数,kppm对象上的一些spatstat函数会抛出错误。通常它会说Error: Covariates ‘nest’ and ‘nest2’ were not found
。
答案 0 :(得分:1)
improve.kppm
的逻辑流程中存在错误:如果vcov=TRUE
和type != "quasi"
,则未定义变量gminus1
。我们会尽快在spatstat
的开发版本中修复此问题。
您是否可以在type="clik1"
的原始来电中选择type="wclik1"
或kppm
?
目前,您可以通过在调用type
和type="quasi"
时省略参数kppm
或明确选择improve.kppm
来避免此错误。< / p>
kppm
无法找到协变量的第二个问题似乎是一个范围问题,但我无法在这里重现它。如果你能提供一个最小的工作实例,它会有所帮助。