Сохранить, было ли предупреждение о конвергенции в lmer, запущенном через afex (R) - PullRequest
0 голосов
/ 25 мая 2020

Я хочу запустить много линейных смешанных моделей (по одной на момент времени временного ряда). Я планирую переставить и исправить множественные сравнения позже. Минимальный пример выводит результат ниже. Есть ли способ также выводить предупреждения? В идеале, сохраняя весь текст предупреждения.

data <- data.frame(timepoint =  rep(1:5,times = 800),
                   subj = c(rep("Subj1",times=1000), rep("Subj2",times=1000),rep("Subj3",times=1000), rep("Subj4",times=1000)),
                   signal = runif(4000),
                   factor_one = sample(rep(c(TRUE,FALSE), 2000)),
                   factor_two = sample(rep(c(TRUE,FALSE), 2000))
                   )

lmeResults <- data %>%
    group_by(timepoint) %>% nest() %>% 
    mutate(lme_out = map(data, ~afex::mixed(signal~factor_one+factor_two+ (1|subj), data = .,method='S',progress=F))) %>% 
    mutate(perf = map(lme_out,"anova_table")) %>%
    mutate(F = map(perf,'F')) %>%
    mutate(p = map(perf,'Pr(>F)')) %>% select(-perf,-data,-lme_out) %>% as.data.frame() 

Результат:

> lmeResults
  timepoint                      F                    p
1         1 0.20837217, 0.02361768 0.6481707, 0.8779006
2         2     1.330282, 1.536380 0.2491032, 0.2155235
3         3 0.73357583, 0.07555716 0.3919840, 0.7834823
4         4   0.0513847, 0.6359361 0.8207293, 0.4254244
5         5   0.1806518, 1.0024798 0.6709284, 0.3170151
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...