Модельные матрицы в R, для смешанных моделей
Я пишу свой собственный программный пакет для определенного типа смешанных моделей в R и хотел бы простой и эффективный способ формирования матриц проектирования. Фиксированная часть ("X") не проблема, я использую model.matrix, или же sparse.model.matrix из пакета Matrix, если я хочу его разреженным.
Я ищу эквивалентную функцию, желательно в хорошо протестированном пакете, для создания матрицы дизайна случайных эффектов ("Z").
Например, если я подгоняю модель lmer / glmer с использованием пакета lme4, я могу получить нужную матрицу, вызвав getME(fitted_object, "Zt"), но я не хочу подгонять модель только для того, чтобы получить матрицу "Zt" или Z. В идеале я хотел бы позвонить (например) model.matrix.random(formula) и получить его напрямую для любой действительной формулы.
В ответ на обсуждение в комментариях также приветствуются идеи о том, как можно запрограммировать такую функцию, но я надеялся, что уже будет хорошо протестированная реализация. Я попытался просмотреть исходный код lme4, чтобы увидеть, как они это делают, но не смог найти соответствующий код.