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

时间:2015-11-26 05:03:20

标签: r spatial spatstat

当我尝试提高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

1 个答案:

答案 0 :(得分:1)

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

您是否可以在type="clik1"的原始来电中选择type="wclik1"kppm

目前,您可以通过在调用typetype="quasi"时省略参数kppm或明确选择improve.kppm来避免此错误。< / p>

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