拟合混合效应模型示例

研究人员检验四块随机选择的田地上六种紫花苜蓿的产量。为每块田地记录一种紫花苜蓿的产量。

研究人员希望了解紫花苜蓿的品种是否会影响平均产量。该研究人员具有四块可以从中收集数据的田地。但是,研究人员希望可以对不在试验中的田地内紫花苜蓿的生长方式进行建模。因此,研究人员使紫花苜蓿生长所在的田地成为随机因子。研究人员使用混合效应模型来同时评估固定效应和随机效应。

  1. 打开样本数据紫花苜蓿.MTW
  2. 选择统计 > 方差分析 > 混合效应模型 > 拟合混合效应模型
  3. 响应中,输入产量
  4. 随机因子(必需)中,输入田地
  5. 固定因子中,输入品种
  6. 单击图形
  7. 图中的残差中,选择条件标准化
  8. 残差图中,选择四合一
  9. 在每个对话框中单击确定

解释结果

在方差分量表中,田地的 p 值为 0.124。假设检验未显示方差分量不为 0 的证据。误差的方差分量的 p 值为 0.003。因为 p 值小于显著性水平 0.05,所以研究人员可以得出误差的方差分量不为 0 的结论。

固定因子项“品种”的 p 值大约为 0 时,就表明至少一种紫花苜蓿对产量的效应与其他五种类型的效应显著不同。

主效应的系数代表每个水平均值和总体均值之间的差值。例如,与品种 1 相关的紫花苜蓿产量大约比总体均值高 0.385 个单位。此系数的 p 值大约为 0 时,就表明品种 1 对产量的效应与“品种”项的另一个水平效应显著不同。为了确定哪些水平效应在统计意义上相同以及哪些水平效应在统计意义上不同,研究人员计划对项执行多重比较分析。

R2 值表明,模型大约解释产量中 92% 的变异。R-sq(调整)值也较高,值大约为 90.2%。研究人员使用此值比较具有不同预测变量数量的模型。

观测值 1 和 5 是异常观测值,因为它们具有大于 2 的标准化残差。研究人员检查数据以确保这些观测值具有正确的响应值。

正态概率图上的残差大致为一条直线,并且这些点在残差与拟合值图上看似随机散布在 0 周围。

方法

方差估计带约束的极大似然
固定效应的自由度Kenward-Roger

因子信息

因子类型水平数
田地随机41, 2, 3, 4
品种固定61, 2, 3, 4, 5, 6

方差分量

来源变量总和的 %变量标准误Z 值P 值
田地0.07791972.93%0.0675801.1529960.124
误差0.02892427.07%0.0105622.7386130.003
合计0.106843       
-2 对数似然 = 7.736012

固定效应检验

分子自由度分母自由度F 值P 值
品种5.0015.0026.290.000

模型汇总

SR-sqR-sq(调整)AICcBIC
0.17007192.33%90.20%12.5413.52

系数

系数系数标准误自由度T 值P 值
常量3.0945830.1438223.0021.5166920.000
品种         
  10.3854170.07762615.004.9650160.000
  20.1454170.07762615.001.8732870.081
  30.1079170.07762615.001.3902050.185
  4-0.3195830.07762615.00-4.1169380.001
  50.3954170.07762615.005.0938380.000

异常观测值的边际拟合和诊断

观测值产量拟合值残差标准化残差
14.1000003.4800000.6200002.190221R
54.2200003.4900000.7300002.578808R
R  残差大

异常观测值的条件拟合和诊断

观测值产量拟合值残差标准化残差
54.2200003.8953390.3246612.400733R
R  残差大