我在混合模型上做一个Uni主题,其中主要统计数据包是Stata . 提供了一个示例来模拟协方差结构,而不适合任何实际的随机效应:

mixed y sex age sex#c.age, || subj:, noconstant residuals(unstructured, t(t)) reml

其中 subject 是群集变量, t 是每个主题4个时间点的 y 值 . 输出包括固定效应参数,并且除了4个方差和6个协方差之外没有随机效应用于剩余误差相关时间结构 .

我试过在网上搜索,但我看不到用R中的 lmer (或其他任何东西)来指定这个的方法 .