Я недавно обновил пакет MICE и какой-то код, который у меня был раньше, сломался ...
Я не очень часто пользуюсь MICE, поэтому буду крайне признателен за предложения по исправлению
imputed_scabies_model <- with(coded_imputed_scabies,glm(impetigo_active~factor(scabies_infestation)+factor(gender)+factor(agegroups), family = binomial()))
Подгонка модели к вмененному набору данных
imputed_scabies_model_summary <- (summary(pool(imputed_scabies_model)))
Получить сводную статистику
imputed_scabies_model_OR <-exp(cbind(imputed_scabies_model_summary[,"est"],imputed_scabies_model_summary[,"lo 95"],imputed_scabies_model_summary[,"hi 95"]))
Экспонировать оценку, нижний и высокий 95% CI, чтобы получить ИЛИ и 95% CI
imputed_scabies_model_summary <- (cbind(imputed_scabies_model_OR,imputed_scabies_model_summary))
Соберите все воедино - так что у меня будет таблица с ИЛИ, 95% -ным доверительным интервалом, а также исходный результат модели.
Я думаю, проблема в том, что вывод
imputed_scabies_model <- with(coded_imputed_scabies,glm(impetigo_active~factor(scabies_infestation)+factor(gender)+factor(agegroups), family = binomial()))
imputed_scabies_model_summary <- (summary(pool(imputed_scabies_model)))
Используется для получения оценки нижнего и верхнего пределов 95% ДИ в виде отдельных столбцов.
Оценка все еще выходит (теперь она называется оценкой, а не «est» - так легко исправить этот фрагмент кода), но теперь в сводке, скорее всего, есть SE, а не границы CI.
! [Выходные данные imputed_scabies_model_summary <- (итоговая (pool (imputed_scabies_model))))
] <a href="https://i.stack.imgur.com/JmUmH.png" rel="nofollow noreferrer"> 1
Так что, я думаю, мне нужно возвести в степень, а затем сгенерировать CI вручную, но мне было интересно, смогут ли люди, которые используют MICE чаще, дать какое-то руководство.
Надеюсь, мой запрос понятен. Я могу предоставить набор данных, если это будет полезно.