SPSS 命令: SPSS GENLIN

在浏览器中拟合广义线性模型(GENLIN)

广义线性模型通过联接函数把结果均值与预测变量联系起来,并允许结果服从正态以外的分布,如计数、右偏的正值费用、比例或是/否事件。

免费 · 无需注册 · 文件留在本机 · 更新于 2026年10月18日

能告诉你什么

Crosstabs 在你的设备上运行与 SPSS GENLIN 对应的过程,提供七种分布、多种联接函数与尺度参数方法、III 类 Wald 或似然比检验、稳健标准误,以及用于重复测量或聚类数据的广义估计方程(GEE)。

  • 整个模型是否优于仅含截距的模型(综合检验中的似然比卡方)。
  • 哪些效应起作用:对每个因子与协变量给出 III 类 Wald 或似然比卡方检验。
  • 通过 Exp(B) 得到结果尺度上的效应:泊松与负二项为率比,二项 Logit 为优势比,对数联接为均值之比。
  • 所选分布是否合适:偏差与 Pearson 卡方除以自由度可提示过度离散,AIC、AICC、BIC 与 CAIC 可用于比较候选模型。
  • 使用 GEE 时,可得到重复测量下的总体平均效应、估计的工作相关矩阵以及 QIC 与 QICC 准则。

适用场景

适合使用

  • 工单数、访问次数、缺陷数等计数(泊松;方差大于均值时用负二项)。
  • 理赔金额、持续时间等右偏正值(伽马或逆高斯,通常配对数联接),或含大量精确零值的金额(Tweedie)。
  • 二元结果或“事件数/试验次数”数据,希望在同一框架下选用 Logit、Probit、对数或互补双对数联接。
  • 同一受试者的重复测量,使用 GEE 并选择可交换、AR(1) 或非结构化工作相关矩阵。

改用其他方法

  • 结果连续且近似正态、预测变量不多:线性回归或 GLM 单变量给出相同估计,输出也更常见。
  • 需要个体层面的随机效应而非总体平均效应:使用线性混合模型。
  • 需要交互项:本面板只拟合因子与协变量的主效应。
  • GEE 与频率权重同时使用:重复测量模型会拒绝频率权重,请先展开数据或清除权重。

在 Crosstabs 中如何操作

  1. 01

    选择广义线性模型

    在“统计”模式下打开数据,在“一般线性模型”组中选择“广义线性模型”(GENLIN)。

  2. 02

    选择模型类型

    在“模型类型”中选择“分布”和“联接函数”,联接函数列表随分布变化。幂联接需设置“幂”,Tweedie 需设置“Tweedie 幂(1–2)”,负二项需设置“负二项参数 k”(估计值或指定值)。

  3. 03

    设置响应与预测变量

    二项分布需设置“响应格式”:“二元(两个类别)”并选择“参考类别”,或“事件数 / 试验次数”并指定“试验次数”变量。然后选择“因变量”并勾选“因子”与“协变量”。

  4. 04

    选择估计与输出

    在“估计与输出”中设置“包含截距”、“尺度参数方法”(最大似然估计、固定值(1)、偏差或 Pearson 卡方)、“稳健(三明治)标准误”、“III 类检验统计量”(Wald 或似然比)以及“包含指数参数估计值”。

  5. 05

    按需设置重复测量

    勾选“使用广义估计方程”,选择“受试者变量”、可选的“主体内变量”、“工作相关矩阵”(独立、可交换、AR(1)、非结构化)和“协方差矩阵”(稳健估计量或基于模型的估计量),然后点击“运行”。

计算示例

以下结果由 Crosstabs 引擎在构建本页时根据模拟数据实际计算得出,并非手工填写。

用示例数据试一试

数据

120 个模拟客户账户,分属三种套餐(Free、Plus、Team),席位数为 5–50。每月工单数为泊松计数,其对数均值随席位数增加,且 Free 套餐最高;随机种子固定。设置“分布”=泊松、“联接函数”=对数,并勾选“包含指数参数估计值”。

拟合优度
项目值自由度值/自由度
偏差148.1091161.277
Pearson 卡方135.8491161.171
对数似然−249.115
赤池信息准则 (AIC)506.230
贝叶斯信息准则 (BIC)517.380
综合检验
似然比卡方自由度显著性
57.5023<.001
模型效应检验
源Wald 卡方自由度显著性
(截距)25.6901<.001
plan14.2072<.001
seats33.7821<.001
参数估计值
参数B标准误Wald 卡方显著性Exp(B)Exp(B) 的 95% Wald 置信区间下限Exp(B) 的 95% Wald 置信区间上限
(截距)0.4110.1477.773.0051.5091.1302.014
[plan=Free]0.4520.12014.207<.0011.5711.2421.987
[plan=Plus]0.2860.1245.295.0211.3311.0431.698
seats0.0220.00433.782<.0011.0231.0151.030

Team 为参考套餐;其冗余行及固定为 1 的(尺度)行已省略。

结果解读

泊松模型相对仅含截距模型的似然比 χ²(3) = 57.50,p < .001。偏差/自由度 = 1.277,Pearson 卡方/自由度 = 1.171,均接近 1,过度离散问题不大。席位每增加 1 个,预期工单数乘以 1.023;套餐整体效应显著(Wald χ²(2) = 14.21,p < .001),同等规模下 Free 账户的工单数是 Team 账户的 1.57 倍。

输出表格说明

模型信息
实际拟合的因变量、概率分布与联接函数;二元响应还会注明建模类别与参考类别。
拟合优度
偏差、标度偏差、Pearson 卡方及其值/自由度,对数似然,以及 AIC、AICC、BIC、CAIC(越小越好)。泊松模型的值/自由度明显大于 1 提示过度离散。
综合检验
比较拟合模型与仅截距模型的似然比卡方。
模型效应检验
每个因子与协变量的 III 类 Wald 或似然比卡方,以及自由度和显著性。
参数估计值
每个参数的 B、标准误、Wald 置信区间与假设检验,尺度参数或负二项 k 行;如请求,还给出 Exp(B) 及其区间。
工作相关矩阵
使用 GEE 时:在所选结构下估计的各测量时点间相关;同时输出相关数据摘要,拟合优度表给出 QIC/QICC。

前提假设与注意事项

  • 方差符合所选分布的均值–方差关系(如泊松方差 = 均值);偏差/自由度远大于 1 时,改用负二项或以 Pearson 卡方估计尺度。
  • 联接函数使预测变量在变换后的均值上呈线性作用;超出分布取值范围的响应(负计数、伽马分布中的 0)会被排除并注明。
  • 观测相互独立;使用 GEE 时要求受试者之间独立,受试者内部允许相关。
  • 尺度参数的最大似然估计仅适用于正态、伽马、逆高斯与 Tweedie 分布。

结果如何验证

GENLIN 与 statsmodels GLM 对照;负二项与 Tweedie 模型另与 R MASS::glm.nb 及 statmod/tweedie 核对。测试涵盖泊松(最大似然、偏差、Pearson 尺度,稳健标准误与频率权重)、伽马、逆高斯、正态、二元 Logit、事件/试验 Probit、估计 k 与固定 k 的负二项以及 Tweedie 1.5:对数似然与信息准则误差在 1e-9 以内,偏差与 Pearson 卡方在 1e-8 以内,稳健标准误在 1e-6 以内。GEE 与 statsmodels GEE 及按 SPSS 矩估计公式实现的 NumPy 版本对照:系数在 1e-8 以内,稳健标准误在 1e-7 以内,工作相关在 1e-8 以内,QIC/QICC 在 1e-7 以内。 验证证据(英文)

常见问题

计数数据用泊松还是负二项?

先拟合泊松模型,查看拟合优度表中偏差与 Pearson 卡方的值/自由度。若明显大于 1,说明计数的变异超过泊松允许的范围;可改用负二项(估计 k)并比较 AIC,或保留泊松但以 Pearson 卡方估计尺度。

如何对率(如每席位工单数)建模?

面板没有偏移量(offset)设置。可在对数联接的泊松模型中把暴露量的对数作为协变量纳入,其系数是估计出的暴露效应而非固定为 1。对于二项比例,请使用“事件数 / 试验次数”。

对数联接下 Exp(B) 表示什么?

它是期望值之比:1.25 表示在其他预测变量不变时,预测变量每增加一个单位(或该因子水平相对参考水平),均值计数(或金额)高 25%。

什么情况下应使用 GEE 而非普通模型?

当同一受试者贡献多行数据时,如重复访问、追踪调查的多个波次、诊所内的多名患者。GEE 保留广义线性模型的均值结构,但使用工作相关矩阵和稳健标准误,使推断考虑到聚类。

III 类检验选 Wald 还是似然比?

似然比检验对每个效应重新拟合去掉该效应的模型,在小样本或接近边界时表现更好;Wald 检验计算更快,是默认选项。使用 GEE 时只提供 Wald 检验。

用你自己的数据运行广义线性模型

打开 .sav、CSV 或 Excel 文件。所有计算都在浏览器中完成,文件不会离开你的设备。

打开分析工作区

SPSS 是 IBM Corp. 的商标。Crosstabs 与 IBM 无关联,也未获得 IBM 认可。