R语言计算人群归因危险度PAR的代码

PAR（Population Attributable Risk）是指在人群中，某种因素引起某种疾病的比例，常用于评估某种因素对人群中疾病的贡献。在R语言中，可以使用epitools包中的函数epitab()和epipop()来计算人群归因危险度PAR。

下面是一个简单的例子，假设我们有一个人群数据集popdata，其中包含了某种因素exposure和某种疾病outcome的信息：

library(epitools)

# 创建人群数据集
popdata <- data.frame(exposure = c(0, 0, 1, 1),
                      outcome = c(0, 1, 0, 1))

# 计算人群暴露比例和疾病发生比例
pexposed <- epipop(popdata$exposure)
poutcome <- epipop(popdata$outcome)

# 计算人群归因危险度PAR
par <- (pexposed - pexposed * poutcome) / poutcome

在这个例子中，我们首先使用epipop()函数计算了人群中暴露和发病的比例，然后使用公式PAR = (p1-p0*p2)/p2来计算人群归因危险度PAR。其中p1表示暴露组中患病比例，p0表示未暴露组中患病比例，p2表示整个人群中患病比例。