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