У меня есть такие данные:
df <- structure(list(`Orientación dicotómica` = c("Neurogastro", "Neurogastro",
"Neurogastro", "Neurogastro", "No neurogastro", "No neurogastro",
"No neurogastro", "No neurogastro", "No neurogastro"), `Fisiopatología más frecuente variante constipación` = c("Más de una variante",
"Obstrucción del tracto de salida", "Tránsito lento / Inercia",
"Transito normal", "Más de una variante", "Obstrucción del tracto de salida",
"Tránsito lento / Inercia", "Transito normal", "Uso de fármacos"
), n = c(22L, 8L, 12L, 11L, 108L, 12L, 101L, 25L, 1L), Proporcion = c(41.5,
15.1, 22.6, 20.8, 43, 4.8, 40.2, 10, 0.4), ds = c(11.5, 11.5,
11.5, 11.5, 19.6, 19.6, 19.6, 19.6, 19.6), IC25 = c(32, 5.6,
13.1, 11.3, 29.8, -8.4, 27, -3.2, -12.8), IC75 = c(51, 24.6,
32.1, 30.3, 56.2, 18, 53.4, 23.2, 13.6)), row.names = c(NA, -9L
), class = c("tbl_df", "tbl", "data.frame"), na.action = structure(c(`10` = 10L), class = "omit"))
Это выглядит так:
И я пытаюсь запустить chisq-тест с сгруппированными данными, чтобы проанализировать, есть ли статистически значимая разница между нейрогастро и отсутствием нейрогастро в отношении frecuente fisiopatología más frecuente.
Как опубликовано в тесте Chi -Square с сгруппированными данными в dplyr
я пытался
test <- df %>%
group_by(`Fisiopatología más frecuente variante constipación`) %>%
summarise(pval = chisq.test(Proporcion,`Orientación dicotómica`)$p.value)
Но я получаю сообщение об ошибке, которое, как мне кажется, связано с тем, что не существует Uso de fármacos для нейрогастро. Мой подход в порядке? Как я могу запустить тест для остальных групп, кроме отсутствующей?
Спасибо!
tryCatch
для обработки ошибокpval = tryCatch(chisq.test(...), error = function(e) NA)
- person rawr   schedule 25.05.2021summarise()
pval
. x 'x' и 'y' должны иметь как минимум 2 уровня i Входpval
-tryCatch(chisq.test(Proporcion,
Orientación dicotómica)$p.value)
. i Ошибка произошла в группе 5: Fisiopatología más frecuente variante constipación = Uso de fármacos. - person Ale Rey   schedule 26.05.2021error = ...
часть - person rawr   schedule 26.05.2021test <- fisiopato_constipacion_especialidad %>% group_by(
Fisiopatología más frecuente variante constipación) %>% summarise(pval = tryCatch(chisq.test(Proporcion,
Orientación dicotómica)$p.value), error = function(e) NA)
. Результат выше ... - person Ale Rey   schedule 26.05.2021summarise(pval = tryCatch(chisq.test(...)$p.value, error = function(e) NA))
- person rawr   schedule 26.05.2021