Я новичок в RI, использую пакет epicalc для расчета коэффициента шансов.У меня есть эти данные
Реальные значения
values <- matrix(c(5,0,35,2,30,38,40,8,11,7,10,9,35,8,66,72,34,4,31,14), nrow = 2)
values
colnames(values) <- c("Laikipia", "West Pokot", "Garissa", "Isiolo", "Mandera", "Marsabit", "Samburu", "Tana River", "Turkana","Wajir")
rownames(values) <- c("0","1")
values <- as.table(values)
values
cc(outcome=NULL, exposure=NULL, cctable=values)}
Но я получаю бесконечное соотношение шансов, что вполне понятно, так как у меня в Лайкипии одно нулевое значение.Поэтому я попробовал поправку Холдейна-Анскомба и добавил 0,5 в каждом коде ячейки, приведенном ниже, но все равно получаю бесконечные значения
Исправление Холдейна-Анскомба
values1 <- matrix(c(5.5,0.5,35.5,2.5,30.5,38.5,40.5,8.5,11.5,7.5,10.5,9.5,35.5,8.5,66.5,72.5,34.5,4.5,31.5,14.5), nrow = 2)
values1
colnames(values1) <- c("Laikipia", "West Pokot", "Garissa", "Isiolo", "Mandera", "Marsabit", "Samburu", "Tana River", "Turkana","Wajir")
rownames(values1) <- c("O","1")
values <- as.table(values1)
values1
cc(outcome=NULL, exposure=NULL, cctable=values1)
"cc {epicalc}"Расчет коэффициентов и построение графиков
Description
Odds ratio calculation and graphing
Usage
cc(outcome, exposure, decimal = 2, cctable = NULL, graph = TRUE,
original = TRUE, design = "cohort", main, xlab = "auto", ylab,
alpha = .05, fisher.or = FALSE, exact.ci.or = TRUE)
Как решить эту проблему.Спасибо