У меня есть следующие чувствительность и специфичность для двух диагностических c тестов, которые я применил к одним и тем же людям
library(epiR)
#TEST1
dat1 <- as.table(matrix(c(46,0,2,39), nrow = 2, byrow = TRUE))
colnames(dat1) <- c("Dis+","Dis-")
rownames(dat1) <- c("Test+","Test-")
rval <- epi.tests(dat1, conf.level = 0.95)
print(rval);
Outcome + Outcome - Total
Test + 46 0 46
Test - 2 39 41
Total 48 39 87
Point estimates and 95 % CIs:
---------------------------------------------------------
Apparent prevalence 0.53 (0.42, 0.64)
True prevalence 0.55 (0.44, 0.66)
Sensitivity 0.96 (0.86, 0.99)
Specificity 1.00 (0.91, 1.00)
Positive predictive value 1.00 (0.92, 1.00)
Negative predictive value 0.95 (0.83, 0.99)
Positive likelihood ratio Inf (NaN, Inf)
Negative likelihood ratio 0.04 (0.01, 0.16)
---------------------------------------------------------
#TEST2
dat2 <- as.table(matrix(c(42,0,6,39), nrow = 2, byrow = TRUE))
colnames(dat2) <- c("Dis+","Dis-")
rownames(dat2) <- c("Test+","Test-")
rval <- epi.tests(dat2, conf.level = 0.95)
Outcome + Outcome - Total
Test + 42 0 42
Test - 6 39 45
Total 48 39 87
Point estimates and 95 % CIs:
---------------------------------------------------------
Apparent prevalence 0.48 (0.37, 0.59)
True prevalence 0.55 (0.44, 0.66)
Sensitivity 0.88 (0.75, 0.95)
Specificity 1.00 (0.91, 1.00)
Positive predictive value 1.00 (0.92, 1.00)
Negative predictive value 0.87 (0.73, 0.95)
Positive likelihood ratio Inf (NaN, Inf)
Negative likelihood ratio 0.12 (0.06, 0.26)
---------------------------------------------------------
Я хочу вычислить объединенную глобальную чувствительность и глобальную специфичность с ее доверительными интервалами. Как я могу это сделать?
Спасибо