У меня проблема с purrr :: map и broom :: tidy при запуске логлинейных GLM в R. По какой-то причине p-значения модели не печатаются при запуске многих моделей, но печатаются с одной моделью.В конце я хотел бы, чтобы несколько моделей печатали значения p для каждой модели, как это происходит в случае одной модели.В приведенном примере используется встроенный набор данных «Титаник» (см. Веб-сайт Уильяма Кинга ).
data(Titanic)
#convert to data frame
T.df <- as.data.frame(Titanic)
head(T.df)
#run glm as loglinear model
model1 <- glm(Freq ~ Sex * Survived, family = poisson, data = T.df)
#print model with tidy--p-values print here
broom::tidy(anova(model1, test = "Chisq"))
#Now run multiple models by class
#Note the models print just fine but without p values
T.df %>%
tidyr::nest(-Class) %>%
dplyr::mutate(
fit = purrr::map(data, ~ anova(glm(Freq ~ Sex * Survived, family = poisson, data = .x)), test="Chisq"),
tidied = purrr::map(fit, broom::tidy)
) %>%
tidyr::unnest(tidied)
Пока я размышляю над этим, как можно помешать метле :: приборке печатать предупреждающие сообщения о нераспознанных столбцах?
Заранее спасибо.