基于队列数据的RR及其95%CI计算

mytable=xtabs(~XAA+Class,data=mydata)#XAA为列变量名，即你的暴露因素，mydata为你的数据框
p.table=prop.table(mytable, 1)#其中‘1’为按照行计算发生率
chisq.test(mytable)#分析改暴露因素是否在组间存在差异
RR=p.table[2,2]/p.table[1,2]#计算RR值，该公式适合二分类的暴露变量，如果多分类，请自行举一反三
RR#输出值
a=1/mytable[1,1]+1/mytable[1,2]+
  1/mytable[2,1]+1/mytable[2,2]#计算lnRR的方差，具体公式流行病的书上有
CI=1.96*sqrt(a)#lnRR的95%CI
b1=log(RR)+CI
b2=log(RR)-CI
exp(b1)#反对数化
exp(b2)