问题描述
我的目标是提出空间相对风险的变化。 relrisk软件包中的spatstat函数可以选择使用se参数来计算相对风险的标准误差。相对风险通常用对数表示，虽然estimate函数输出中的relrisk可以转换，但是SE不能转换。
例如，相对风险：
f1 <- spatstat::relrisk(spatstat.data::chorley,relative = TRUE,se = TRUE)
plot(f1$estimate); plot(f1$SE)
但是，对于对数相对风险，不能以与相对风险估计相似的方式转换标准误差。
plot(log(f1$estimate))
有人可以建议计算解决方案吗？还是有可能将此功能添加到relrisk函数中，也许作为逻辑参数log来估计日志相对风险及其标准误差？谢谢！
解决方法
使用增量法，log(R)的标准误差约为se(R)/R，即R的标准误差除以R的估计值。所以最快的解决方案是
plot(with(f1,se/estimate))