能告诉你什么
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 中如何操作
- 01
选择广义线性模型
在“统计”模式下打开数据,在“一般线性模型”组中选择“广义线性模型”(GENLIN)。
- 02
选择模型类型
在“模型类型”中选择“分布”和“联接函数”,联接函数列表随分布变化。幂联接需设置“幂”,Tweedie 需设置“Tweedie 幂(1–2)”,负二项需设置“负二项参数 k”(估计值或指定值)。
- 03
设置响应与预测变量
二项分布需设置“响应格式”:“二元(两个类别)”并选择“参考类别”,或“事件数 / 试验次数”并指定“试验次数”变量。然后选择“因变量”并勾选“因子”与“协变量”。
- 04
选择估计与输出
在“估计与输出”中设置“包含截距”、“尺度参数方法”(最大似然估计、固定值(1)、偏差或 Pearson 卡方)、“稳健(三明治)标准误”、“III 类检验统计量”(Wald 或似然比)以及“包含指数参数估计值”。
- 05
按需设置重复测量
勾选“使用广义估计方程”,选择“受试者变量”、可选的“主体内变量”、“工作相关矩阵”(独立、可交换、AR(1)、非结构化)和“协方差矩阵”(稳健估计量或基于模型的估计量),然后点击“运行”。
计算示例
以下结果由 Crosstabs 引擎在构建本页时根据模拟数据实际计算得出,并非手工填写。
数据
120 个模拟客户账户,分属三种套餐(Free、Plus、Team),席位数为 5–50。每月工单数为泊松计数,其对数均值随席位数增加,且 Free 套餐最高;随机种子固定。设置“分布”=泊松、“联接函数”=对数,并勾选“包含指数参数估计值”。
| 项目 | 值 | 自由度 | 值/自由度 |
|---|---|---|---|
| 偏差 | 148.109 | 116 | 1.277 |
| Pearson 卡方 | 135.849 | 116 | 1.171 |
| 对数似然 | −249.115 | ||
| 赤池信息准则 (AIC) | 506.230 | ||
| 贝叶斯信息准则 (BIC) | 517.380 |
| 似然比卡方 | 自由度 | 显著性 |
|---|---|---|
| 57.502 | 3 | <.001 |
| 源 | Wald 卡方 | 自由度 | 显著性 |
|---|---|---|---|
| (截距) | 25.690 | 1 | <.001 |
| plan | 14.207 | 2 | <.001 |
| seats | 33.782 | 1 | <.001 |
| 参数 | B | 标准误 | Wald 卡方 | 显著性 | Exp(B) | Exp(B) 的 95% Wald 置信区间下限 | Exp(B) 的 95% Wald 置信区间上限 |
|---|---|---|---|---|---|---|---|
| (截距) | 0.411 | 0.147 | 7.773 | .005 | 1.509 | 1.130 | 2.014 |
| [plan=Free] | 0.452 | 0.120 | 14.207 | <.001 | 1.571 | 1.242 | 1.987 |
| [plan=Plus] | 0.286 | 0.124 | 5.295 | .021 | 1.331 | 1.043 | 1.698 |
| seats | 0.022 | 0.004 | 33.782 | <.001 | 1.023 | 1.015 | 1.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 检验。
相关过程
- GLM 单变量UNIANOVA
- 二元 Logistic 回归LOGISTIC REGRESSION
- 线性混合模型MIXED
- 概率单位分析PROBIT
- 一般对数线性分析GENLOG
- 复杂抽样一般线性模型CSGLM
- 免费的 SPSS 替代工具:支持范围
用你自己的数据运行广义线性模型
打开 .sav、CSV 或 Excel 文件。所有计算都在浏览器中完成,文件不会离开你的设备。
SPSS 是 IBM Corp. 的商标。Crosstabs 与 IBM 无关联,也未获得 IBM 认可。