r - 如何计算R中的相对风险？

Question

如何计算R中的相对风险？

在输入中我有

names(DS)
Tab<-table(DS[,5],DS[,11],DS[,3])

      No Yes
  No   4  16
  Yes  40 168

我是 R 编程语言的新手...

Answer 1

这是一个相当简单的计算；相对风险只是(pos1/total1)/(pos2/total2)哪里pos1是第一组的病例数，pos2在第二组中，和total变量是组总数

但是，您可能对 epitab 感兴趣。 epitools的功能包裹:

## maybe:
##   install.packages("epitools")
library("epitools")

见 ?epicalc对于输入格式，但这应该适用于您的示例:

tab <- matrix(c(4,16,40,168),byrow=TRUE,nrow=2)
epitab(tab,method="riskratio")
## $tab
##           Outcome
## Predictor  Disease1        p0 Disease2        p1 riskratio     lower    upper
##   Exposed1        4 0.2000000       16 0.8000000  1.000000        NA       NA
##   Exposed2       40 0.1923077      168 0.8076923  1.009615 0.8030206 1.269361
##           Outcome
## Predictor  p.value
##   Exposed1      NA
##   Exposed2       1

您当然应该仔细检查结果是否有意义；例如 p0第一组为 4/20=0.2，第二组为 40/208=0.192。风险比为((16/20)/(168/208))=0.8/0.8076=1.009。