文章里面的annual incidence rates及95%CI如何计算?
可以给你一个示例,因为没有提供太多的有效信息;你可以根据这个示例修改一下
# 安装并加载binom包
if(!require(binom)) {
install.packages("binom")
library(binom)
}
# 输入数据
new_cases <- 100 # 新发病例数
population <- 10000 # 总人口数
# 计算年发病率
annual_incidence_rate <- new_cases / population * 100000
# 计算95%置信区间
confint <- binom.confint(new_cases, population, conf.level = 0.95, methods = "exact")
# 输出结果
print(paste("Annual incidence rate: ", annual_incidence_rate))
print(paste("95% CI: [", confint$lower * 100000, ", ", confint$upper * 100000, "]"))