spatstat 中的相对风险估计

Relative risk estimation in spatstat

在计算两个点模式的相对风险估计 (relrisk.ppp) 时,我 运行 遇到了问题:一个在矩形区域中有四个标记,另一个在圆形区域中有两个标记.

对于具有四个标记的第一个模式,我能够在 large imlist 中获得相对风险和结果对象,每个标记对应 4 个元素。

然而,对于第二个模式,它给出了一个包含 10 个元素的列表,其中第一个矩阵 v 是空的,有 NA 个条目。当创建的点模式对象看起来相同时,我正在为可能出现的错误而烦恼。任何帮助将不胜感激。谢谢

对于您的第一个数据集,结果是图像对象的列表(class im 的四个对象的列表)。对于您的第二个数据集,relrisk.ppp 的结果是单个图像(class im 的对象)。当只有两种可能的点类型(两个可能的标记值)时,这是默认行为。参见 help(relrisk.ppp)

在所有情况下,您应该能够 plotprint 结果对象。您不需要检查图像的内部数据。

更多解释:当只有两种可能的点类型时,relrisk.ppp的默认行为是将它们视为case-control 数据,其中属于第一类的点作为对照(如非感染者),属于第二类的点作为病例(如感染者)。强度比(病例除以对照)被估计为图像。

如果您不希望发生这种情况,请设置参数 casecontrol=FALSE 然后 relrisk.ppp 将始终 return 列表图像,每个可能的标记对应一张图像。每个图像都给出了该类型点的空间变化概率。

help(relrisk.ppp)或者书上都有解释