r - 如何将两个数字与置信区间相乘？

Question

如何将两个具有置信区间的数字相乘并得到新置信区间的结果？

我的情况如下。

一种疾病有两种测试，测试 A 和测试 B。

这是测试 A 的混淆矩阵。

data.frame(disease=c(30,70),No_Disease=c(10,90),row.names =c("Test A positive","Test A negative"))

可以使用下面的代码计算正似然值。

library(epiR)
datA <- as.table(matrix(c(30,70,10,90), nrow = 2))
colnames(datA) <- c("Dis+","Dis-")
rownames(datA) <- c("Test+","Test-")
rvalA <- epi.tests(datA, conf.level = 0.95)
print(rvalA)

阳性似然比 3.00 (1.55, 5.80)

测试 B 混淆矩阵。

data.frame(disease=c(60,40),No_Disease=c(40,60),row.names =c("Test B positive","Test B negative"))

检验 B 的正似然值。

datB <- as.table(matrix(c(60,40,40,60), nrow = 2))
colnames(datB) <- c("Dis+","Dis-")
rownames(datB) <- c("Test+","Test-")
rvalB <- epi.tests(datB, conf.level = 0.95)
print(rvalB)

阳性似然比 1.50 (1.12, 2.00)

假设有一个病人A检测阳性，然后B检测阳性。我想将这些正似然比相乘，例如 3 x 1.5，以获得赔率。但我不知道如何在这个乘法之后得到置信区间。如果可能的话，你能用代码教我如何做到这一点吗？先感谢您！

标签： rstatisticsconfidence-interval