spatstat 中的相对风险估计
Relative risk estimation in spatstat
在计算两个点模式的相对风险估计 (relrisk.ppp
) 时,我 运行 遇到了问题:一个在矩形区域中有四个标记,另一个在圆形区域中有两个标记.
对于具有四个标记的第一个模式,我能够在 large imlist
中获得相对风险和结果对象,每个标记对应 4 个元素。
然而,对于第二个模式,它给出了一个包含 10 个元素的列表,其中第一个矩阵 v
是空的,有 NA
个条目。当创建的点模式对象看起来相同时,我正在为可能出现的错误而烦恼。任何帮助将不胜感激。谢谢
对于您的第一个数据集,结果是图像对象的列表(class im
的四个对象的列表)。对于您的第二个数据集,relrisk.ppp
的结果是单个图像(class im
的对象)。当只有两种可能的点类型(两个可能的标记值)时,这是默认行为。参见 help(relrisk.ppp)
。
在所有情况下,您应该能够 plot
和 print
结果对象。您不需要检查图像的内部数据。
更多解释:当只有两种可能的点类型时,relrisk.ppp
的默认行为是将它们视为case-control 数据,其中属于第一类的点作为对照(如非感染者),属于第二类的点作为病例(如感染者)。强度比(病例除以对照)被估计为图像。
如果您不希望发生这种情况,请设置参数 casecontrol=FALSE
然后 relrisk.ppp
将始终 return 列表图像,每个可能的标记对应一张图像。每个图像都给出了该类型点的空间变化概率。
help(relrisk.ppp)
或者书上都有解释
在计算两个点模式的相对风险估计 (relrisk.ppp
) 时,我 运行 遇到了问题:一个在矩形区域中有四个标记,另一个在圆形区域中有两个标记.
对于具有四个标记的第一个模式,我能够在 large imlist
中获得相对风险和结果对象,每个标记对应 4 个元素。
然而,对于第二个模式,它给出了一个包含 10 个元素的列表,其中第一个矩阵 v
是空的,有 NA
个条目。当创建的点模式对象看起来相同时,我正在为可能出现的错误而烦恼。任何帮助将不胜感激。谢谢
对于您的第一个数据集,结果是图像对象的列表(class im
的四个对象的列表)。对于您的第二个数据集,relrisk.ppp
的结果是单个图像(class im
的对象)。当只有两种可能的点类型(两个可能的标记值)时,这是默认行为。参见 help(relrisk.ppp)
。
在所有情况下,您应该能够 plot
和 print
结果对象。您不需要检查图像的内部数据。
更多解释:当只有两种可能的点类型时,relrisk.ppp
的默认行为是将它们视为case-control 数据,其中属于第一类的点作为对照(如非感染者),属于第二类的点作为病例(如感染者)。强度比(病例除以对照)被估计为图像。
如果您不希望发生这种情况,请设置参数 casecontrol=FALSE
然后 relrisk.ppp
将始终 return 列表图像,每个可能的标记对应一张图像。每个图像都给出了该类型点的空间变化概率。
help(relrisk.ppp)
或者书上都有解释