Я хотел бы получить выборки из апостериорного распределения случайных эффектов. Я предполагаю, что это делается внутри функции simulate. Есть ли прямой способ сделать это?
В https://cran.r-project.org/web/packages/lme4/vignettes/lmer.pdf (стр. 29, уравнения 57-59) виньетка дает выражение для задней («условной») ковариационной матрицы , Я думаю, что это можно вычислить, используя переменные, извлеченные из getME, но виньетка говорит, что «на практике обратное в уравнении 57 вычисляется не напрямую, а скорее используется эффективный метод, который использует разреженные структуры. можно легко воспроизвести?