SPSS 命令: SPSS CSORDINAL

复杂调查数据的有序回归

复杂抽样有序回归用累积关联模型对有序调查结果(例如满意度评分)建模,其阈值和系数描述总体。

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

能告诉你什么

Crosstabs 在你的设备上运行与 SPSS CSORDINAL 对应的过程,提供 logit、probit、互补双对数、负双对数和 Cauchit 关联函数、伪似然估计、基于设计的标准误以及基于设计自由度的 Wald 检验。

  • 预测变量是使个案偏向更高还是更低的类别:位置系数为正表示更高类别更可能出现。
  • 相邻类别之间在潜变量尺度上的阈值(切点)。
  • 每个预测变量整体上是否有作用:基于设计自由度的 Wald F 与 Rao–Scott 调整检验。
  • 抽样设计使每个参数的方差改变了多少(设计效应),以及模型拟合程度(伪 R²)。

适用场景

适合使用

  • 在加权整群调查中,以李克特评分、满意度量表或有序分档(收入档、严重程度)为结果变量。
  • 希望每个预测变量只有一个系数,而不是像多项 Logistic 那样每个类别各有一组系数。
  • 比较不同关联函数:潜变量近似正态时用 probit;高类别更常见时用互补双对数;低类别更常见时用负双对数。

改用其他方法

  • 类别没有自然顺序:使用复杂抽样 Logistic 回归(多项)。
  • 需要平行线(比例优势)检验:此处不提供;可改为比较多项模型中各类别的系数。
  • 分点很多、近似连续的评分:复杂抽样一般线性模型可能更易解释。

在 Crosstabs 中如何操作

  1. 01

    打开调查数据

    在工作区打开 .sav、CSV 或 Excel 文件并切换到“统计”模式。文件只在浏览器中读取,不会上传。

  2. 02

    选择过程

    在“复杂抽样”组中选择“复杂抽样有序回归”(CSORDINAL)。

  3. 03

    确认已保存的抽样计划

    面板顶部的“抽样计划”一行会显示权重、分层、聚类和估计方法。若此数据集尚未保存计划,请先打开“抽样计划”(CSPLAN)。计划中的变量不会出现在变量选择框中,工作区权重也不会被使用:复杂抽样始终使用计划中的权重。

  4. 04

    选择因变量和关联函数

    选择“因变量(有序)”,类别按取值排序。在“关联函数”中选择“Logit(SPSS 默认)”“Probit”“互补双对数”“负双对数”或“Cauchit”。

  5. 05

    添加预测变量并运行

    勾选“因子(分类)”和“协变量(数值)”,设置“置信水平”后运行。若估计未收敛,会显示警告。

计算示例

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

用示例数据试一试

数据

一份 120 名受访者的模拟住户调查:3 个层(Metro、Town、Rural),每层 4 个初级抽样单元(PSU),每个 PSU 10 名受访者;抽样权重从 Metro 的约 35 增加到 Rural 的约 110。同一 PSU 内的受访者共享一个随机效应,因此同群个案彼此相似。抽样计划为“权重 + 分层 + 聚类”,采用有放回(WR)估计;所有数值均由固定随机种子生成。

模型效应检验
源Wald Fdf1df2显著性调整的 F显著性(调整 F)
sex9.98619.0129.986.012
age8.54519.0178.545.017

Wald F 为 Fellegi 调整 Wald F,自由度为 (r, ν − r + 1);调整的 F 为相对 SRS 协方差的 Rao–Scott 二阶校正。

参数估计值
项目参数B标准误t显著性设计效应
阈值[satisfaction=1]−0.0970.803−0.121.9072.406
阈值[satisfaction=2]1.1070.6641.668.1301.781
阈值[satisfaction=3]2.2480.5434.136.0031.141
阈值[satisfaction=4]3.4940.7374.742.0011.627
位置[sex=Female]0.7410.2353.160.0120.408
位置[sex=Male]ᵃ0
位置age0.0340.0122.923.0171.124

a. 此参数冗余,因此设为零(最后一个类别为参考)。关联函数:Logit。

结果解读

在累积 logit 模型 P(Y ≤ k) = F(阈值_k − xβ) 中,正系数使个案偏向更高的满意度类别。女性的位置系数为 0.741(p = .012),即女性处于更高类别的优势约为男性的 2.10 倍;年龄每增加一岁,系数增加 0.034(p = .017)。四个阈值是五个类别之间的切点。标准误来自基于设计的三明治估计,Nagelkerke R² 为 0.125。

输出表格说明

模型摘要
仅含阈值模型和最终模型的伪 −2 对数似然、Cox-Snell、Nagelkerke 与 McFadden R²、使用的未加权个案数以及估计的总体规模。
模型效应检验
每个预测变量的 Wald F、Wald 卡方和 Rao–Scott 调整 F,基于设计自由度。
参数估计值
阈值行(比类别数少一个)和预测变量的位置行,各附 B、标准误、置信区间、t、自由度、显著性和设计效应;脚注注明关联函数。

前提假设与注意事项

  • 已保存的计划真实反映了样本的抽取方式:正确的抽样权重、第 1 阶段分层和初级抽样单元(PSU)。计划有误,无论模型多好,标准误都是错的。
  • 模型仅包含主效应(不含交互项)。因子采用 SPSS 式编码,以最后一个类别为冗余参考;任一模型变量缺失的个案按成列方式剔除,但仍以零贡献保留在设计中。
  • 每个预测变量在所有阈值处效应相同(平行线假设;logit 关联下即比例优势)。模型为 P(Y ≤ k) = F(阈值_k − xβ)。
  • 推断基于设计自由度(PSU 数减层数)上的 t 和 F 分布。每层 PSU 很少时自由度很小、区间较宽;只有一个 PSU 的层对方差没有贡献。

结果如何验证

在一个模拟的两阶段分层调查(235 个个案)上,针对三种设计,全部五种关联函数都与独立编写的全精度 NumPy/SciPy 累积关联拟合及三明治方差对照:系数相对误差在 1e-7 以内,标准误和设计效应在 1e-6 以内,−2 对数似然在 1e-10 以内。其中三种关联函数还与 R survey 包的 svyolr 比较:估计值误差在 1e-7 以内;由于 svyolr 使用数值 Hessian,标准误的偏差在 2e-3 以内。 验证证据(英文)

常见问题

如何理解系数的符号?

模型为 P(Y ≤ k) = F(阈值_k − xβ),因此正系数会降低处于各类别及以下的概率——即使个案偏向更高的类别。在 logit 关联下,Exp(B) 是处于更高类别的累积优势比。

应选择哪种关联函数?

Logit 是默认选项,可解释为优势比。Probit 假定潜变量服从正态分布;互补双对数适合集中于高类别的结果,负双对数适合集中于低类别的结果,Cauchit 对极端值更稳健。不确定时可比较各关联函数的伪 R²。

为什么没有截距行?

在累积关联模型中,阈值承担截距的作用,每个切点一个。因此每个因子都以最后一个水平作为冗余参考。

有平行线检验吗?

没有。面板只拟合比例(平行线)模型。若怀疑效应在不同阈值处不同,可拟合多项复杂抽样 Logistic 回归,并比较各类别的系数。

我的调查文件会被上传吗?

不会。工作区在浏览器中读取文件并运行 CSORDINAL,抽样计划也只保存在当前浏览器中。数据不会发送到 Crosstabs 服务器,也无需注册账号。

用你自己的数据运行复杂抽样有序回归

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

打开分析工作区

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