关于Matlab coefTest函数中Fixed-effects contrasts的含义与用法问询

固定效应对比（Fixed-effects contrasts）的含义 #

固定效应对比是用于检验线性混合模型中固定效应系数的特定线性组合是否显著不为0的统计方法。它通过构造对比矩阵，将我们关注的组间差异转化为系数的线性假设，以此判断差异是否具有统计学意义。

示例中[0 1 -1]各元素的意义 #

在Shift示例的模型lme = fitlme(tbl,'Rating ~ Shift + (1|Operator)')中，固定效应系数对应如下：

• 第一个系数：截距项（代表夜班（Night）的平均评分，作为基准组）
• 第二个系数：早班（Morning）相对于夜班的效应值
• 第三个系数：晚班（Evening）相对于夜班的效应值

对比矩阵[0 1 -1]对应线性组合：0*截距 + 1*早班效应 + (-1)*晚班效应，也就是早班效应 - 晚班效应。其中：

• 0：表示不纳入截距项的计算
• 1：对应早班的固定效应系数权重
• -1：对应晚班的固定效应系数权重

coefTest(lme,[0 1 -1])检验早班与晚班差异的原理 #

这个调用实际是检验原假设H₀：早班效应 - 晚班效应 = 0，等价于检验早班和晚班的平均评分差异为0。
coefTest会计算该线性组合的t统计量，再结合模型自由度算出对应的p值：

• 若p值小于预设的显著性水平（如0.05），则拒绝原假设，认为早班与晚班的评分存在显著差异；
• 若p值大于显著性水平，则没有足够证据证明两者存在差异。

从均值角度推导：早班平均评分=截距+早班效应，晚班平均评分=截距+晚班效应，两者差值正好等于早班效应 - 晚班效应，因此这个对比矩阵直接对应两组的均值差异检验。

检验晚班与夜班差异的对比矩阵修改 #

晚班与夜班的均值差为：晚班平均评分 - 夜班平均评分 = (截距+晚班效应) - 截距 = 晚班效应。
对应的线性组合为0*截距 + 0*早班效应 + 1*晚班效应，因此对比矩阵应改为[0 0 1]，调用代码为：

pVal = coefTest(lme,[0 0 1])

该调用检验的是原假设H₀：晚班效应=0，即晚班与夜班的评分无显著差异。

内容的提问来源于stack exchange，提问作者Fengting Ji