我相信lme4包使用的语法与rstanarm相同，但在拟合斜率上的分组随机效应时，我很难准确地找出不同选项之间的差异。
假设我们有一个带有公式的模型
response ~ fixed + (1 + rand1 | group)字符串
这与其他三个选项之间的假设差异究竟是什么：
一曰：
response ~ fixed + (0 + rand1 | group)型
二：
response ~ fixed + (-1 + rand1 | group)型
三：
response ~ fixed + (1 + rand1 || group)型
我知道here一个答案提到||公式假设group水平的参数偏差之间的协方差矩阵是对角的，但我不清楚关于协方差矩阵的假设如何在第一个和第二个选项之间变化。
1条答案
按热度按时间vnjpjtjt1#
选项1和选项2应给予 * 相同 * 的结果：从
?formula'-'运算符删除指定的术语...它也可以用来删除截距项：当拟合线性模型时，“y ~ x - 1”指定通过原点的线。没有截距的模型也可以指定为'y ~ x +'或'y ~+ x'
(it是真实的，但不是显而易见的，本说明书适用于此）。
如果
rand1是分类的（即因子）因为
(rand1|group)和（例如）(0+rand1|group)都指定了一个一般的正定协方差矩阵（“非结构化”），所以在这两种情况下，你应该得到相同的 * 总体 * 模型拟合（先验的排除效应）。差异在于协方差矩阵的各个元素的含义。使用(rand1|group)，lme 4设置标准 * 治疗对比 *，因此协方差矩阵的第一行/列对应于各组间截距（基线水平值）的方差;随后的行/列对应于基线水平与水平2、3、…对于(0+rand1|group)，行/列是相应组中测量值的方差。如果
rand1连续这个案子的运作方式不同;这里
rand1 | group给出长度为2（{截距，斜率}）的随机效应，协方差矩阵为2x2;用+0或-1抑制截距将产生标量随机效应模型，该模型假设各组间截距恒定，但斜率变化（即，描述斜率组间方差的单个参数估计值）。