启动“广义线性混合模型”特质

通过以下方式启动“广义线性混合模型”特质：选择分析 > 拟合模型，输入一个或多个 Y 列，然后从特质菜单选择广义线性混合模型。

“拟合模型”启动窗口

您可以指定带有固定效应、随机效应或固定与随机效应组合的模型。还可以指定响应变量的分布。启动窗口中的选项根据您所指定的模型的性质而有所不同。有关“选择列”红色小三角菜单中选项的详细信息，请参见《使用 JMP》中的““列过滤器”菜单”。

使用“广义线性混合模型”特质拟合模型时，您可以允许无界限方差分量。这意味着具有负估计值的方差分量不被报告为零。默认选定该选项。若您关注的是固定效应，则它应保持选定状态，因为将方差估计值的边界定为零将导致固定效应检验发生偏倚。有关“无界限方差分量”选项的详细信息，请参见“负方差”。

图 9.6 选择了“广义线性混合模型”的“拟合模型”启动窗口 

有关对所有特质通用的“拟合模型”窗口各方面的详细信息，请参见“模型规格”。有关“选择列”红色小三角菜单中选项的详细信息，请参见《使用 JMP》中的““列过滤器”菜单”。这里提供“广义线性混合模型”特质特有的信息。

指定分布

在“拟合模型”启动窗口中，将“广义线性混合模型”选作“特质”时即会显示“分布”选项。在此您可以指定 Y 的分布。下面说明了可用的分布。

二项

指定 Y 服从参数为 p 和 n 的二项分布。响应 Y 指示所有试验中具有固定概率 p 的 n 个独立试验中的总成功次数。该分布允许使用样本大小列。若未列出任何列，则假定样本大小为 1。p 的连结函数为 Logit。若选择建模类型为“名义型”的二值响应变量，则只有“二项”是可用的响应分布。

当您选择“二项”作为“分布”时，必须使用以下方式之一指定响应变量。

‒ 未汇总：若您的数据未汇总为事件频数，指定一个二值列作为响应。若该列的建模类型为名义型，那么您可将其中一个水平指定为“目标水平”。默认“目标水平”值是基于水平顺序的两个水平中的较高者。

‒ 已汇总且将样本大小列输入为第二个 Y：若您的数据已汇总为事件（成功）和试验的频数，则按如下顺序将两个连续列指定为 Y：成功次数的计数和试验次数的计数。

Poisson

指定 Y 服从均值为 l 的 Poisson 分布。Poisson 分布通常对给定时间间隔中的事件数建模，经常表示为计数数据。l 的连结函数为对数。即使 Y 采用非整数值，Poisson 回归也是允许的。

“固定效应”选项卡

在“固定效应”选项卡上添加所有固定效应。根据需要使用“添加”、“交叉”、“嵌套”、“宏”和“特性”选项。有关这些选项的详细信息，请参见“模型规格”。请注意，模型中可能没有固定效应。

注意：若连续列涉及随机效应，则该列不会位于中心，即使选择了“模型规格”红色小三角菜单中的“中心多项式”选项也是如此。

图 9.7 显示“固定效应”的完成的“拟合模型”启动窗口 

“随机效应”选项卡

使用“随机效应”选项卡指定传统方差分量模型和随机系数模型。

注意：若连续列涉及随机效应，则该列不会位于中心，即使选择了“模型规格”红色小三角菜单中的“中心多项式”选项也是如此。

方差分量

对于传统方差分量模型，使用“添加”、“交叉”或“嵌套”选项指定如随机区组、整区误差项和子区误差项之类的项。有关这些选项的详细信息，请参见“模型规格”。

图 9.8 显示“随机效应”的已完成“拟合模型”启动窗口 

随机系数

要构造随机系数模型，使用“嵌套随机系数”按钮来创建随机系数组。

1. 从“选择列”列表中选择作为预测变量的连续列。

2. 选择随机效应选项卡，然后点击添加。

3. 在“随机效应”选项卡中选择这些效应。还选择包含随机效应的列，该效应的水平定义各回归模型。

4. 点击嵌套随机系数按钮。

最后一步创建在随机效应的各个水平内有相关性的随机截距和随机斜率效应。由于对象间的变异性，对象嵌套在其他效应内。若您相信截距可能对所有组是固定的，可以选择截距[<组>]&随机系数(1)，然后点击删除。

您可以按该方式定义多个随机系数组，如在层次线性模型中那样。当斜率和截距系数具有随机批次效应和处理*随机批次效应时，这可能是必要的。在以下层次线性模型中这也可能是必要的：学习成绩具有随机学生效应和随机学校效应，并且学生嵌套在学校内。

数据格式

“拟合模型”平台的“广义线性混合模型”特质要求所有响应测量值都包含在一个响应列中。重复测量数据有时记录在多个列中，其中，每行都是一个对象，重复的测量值记录在单独的响应列中。必须在运行“广义线性混合模型”特质之前堆叠该格式的数据。Cholesterol.jmp 和 Cholesterol Stacked.jmp 样本数据表分别演示了宽格式和堆叠格式。请注意，宽表中的每行都对应于堆叠表中患者的一个水平。