流行病学计算

---

title: “[R语言] 流行病学计算”
date: 2022-12-06
lastmod: 2022-12-06
draft: false
tags: [“R语言”,“流行病学”,“统计学”]
toc: true
autoCollapseToc: true

主要来自结果数据计算

R × C 计算卡方

compare<-matrix(c(99,48,29,0,
                  2,3,2,0,
                  0,0,0,11),nr = 3)
compare
a=chisq.test(compare)
a$expected
fisher.test(compare)

OR值计算

compare<-matrix(c(99,48,
                  2,3
                  ),nr = 2)
compare
library(epiR)
# cross.sectional cohort.time
epi.2by2(compare, method = "cross.sectional")

标准发病率

⚠️ 调查人群无该年龄层，也不允许删除

library(epitools)
ageadjust.direct(count = Case,  # 调查发病人数
            pop =  pop,     # 调查人数
            stdpop = standard_age_pop) # 年龄标准人口数