spatstat 中对数相对风险的标准误

Question

我的 objective 是呈现空间相对风险的变化。 spatstat 包中的 relrisk 函数可以选择使用 se 参数计算相对风险的标准误差。相对风险通常用对数表示，虽然 relrisk 函数输出的 estimate 可以转换，但 SE 不能。

例如，对于相对风险：

f1 <- spatstat::relrisk(spatstat.data::chorley, relative = TRUE, se = TRUE)
plot(f1$estimate); plot(f1$SE)

但是，对于对数相对风险，不能以与相对风险估计类似的方式转换标准误差。

plot(log(f1$estimate))

有人可以提出计算解决方案吗？或者是否可以将此功能添加到 relrisk 函数中，也许作为逻辑参数 log 来估计对数相对风险及其标准误差？谢谢！

Answer 1

使用delta法，log(R)的标准误差约为se(R)/R，即R的标准误差除以R的估计值。所以最快的解决方案是

plot(with(f1, se/estimate))

Standard error of log relative risk in spatstat