SPSS 命令: SPSS BAYES REGRESSION

在浏览器中做贝叶斯线性回归

贝叶斯线性回归估计结果变量如何依赖多个预测变量,给出每个系数的完整后验分布,以及模型相对于仅含截距模型的贝叶斯因子。

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

能告诉你什么

Crosstabs 在你的设备上运行与 SPSS BAYES REGRESSION 对应的过程,支持连续协变量和指示符编码的因子、四种 g 先验贝叶斯因子,以及参考先验或共轭先验。

  • 数据在多大程度上支持完整模型优于仅含截距的模型(BF10),并同时给出 R、R² 和调整后 R²。
  • 截距和每个系数的后验众数、后验平均值、后验方差及可信区间。
  • 误差方差的后验分布,直接说明还有多少未被解释的离散程度。
  • 通过共轭先验,了解关于系数大小的先验信息如何改变估计。

适用场景

适合使用

  • 用多个连续或分类预测变量预测一个连续结果,例如根据面积、距离和城区预测租金。
  • 希望系数区间可以按概率解读,或希望纳入先验估计。
  • 需要证明一组预测变量没有作用(BF01 较大),而不只是 F 检验不显著。

改用其他方法

  • 二分类或计数结果:使用二元 Logistic 回归或广义线性模型。
  • 需要每个预测变量的贝叶斯因子或模型平均:面板只报告整个模型的一个贝叶斯因子。
  • 逐步选择变量、残差诊断或共线性统计:使用线性回归。

在 Crosstabs 中如何操作

  1. 01

    打开数据

    打开文件,切换到“统计”模式,在“贝叶斯统计”组中选择“贝叶斯线性回归”。

  2. 02

    设置模型

    选择“因变量(标度)”,在“协变量(标度)”中勾选连续预测变量,在“因子(指示符编码,以最后一个类别为参考)”中勾选分类预测变量。

  3. 03

    选择贝叶斯因子

    在“贝叶斯因子”中选择“JZS(Zellner–Siow)”“Zellner(固定 g,留空则 g = N)”、带参数 a 的“超 g”,或带尺度 s 的“Rouder”。

  4. 04

    选择系数的先验

    保留“参考先验,p(θ, σ²) ∝ 1/σ²”,或选择“共轭正态–逆伽马”,为所有系数设置“形状 a0”“尺度 b0”以及共同的“先验平均值”和“先验方差(× σ²)”。

  5. 05

    运行

    设置“可信区间(%)”,声明工作区权重(如有),然后点击“运行贝叶斯线性回归”。缺失值按列表删除。

计算示例

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

用示例数据试一试

数据

40 套模拟公寓,包括面积(平方米)、到市中心的距离(千米)和月租金,由带噪声的线性模型以固定随机种子生成。

模型摘要
NRR²调整后 R²估算标准误BF10
400.8650.7490.735154.761.718e+9

JZS 方法(Zellner–Siow 柯西先验),模型与仅含截距的模型进行比较。

系数的贝叶斯估计
参数后验众数后验平均值后验方差95% 可信区间下限95% 可信区间上限
(Intercept)491.67491.6713766.945260.45722.90
area14.8414.842.55011.7017.99
distance−33.93−33.9380.698−51.63−16.23

参考先验 p(θ, σ²) ∝ 1/σ²;等尾可信区间。

误差方差的贝叶斯估计
参数后验众数后验平均值95% 可信区间下限95% 可信区间上限
误差方差22723.225320.115919.540089.6

后验为逆伽马(18.5, 443103);为节省宽度,此处省略方差列。

结果解读

两个预测变量解释了租金变异的 74.9%(R² = 0.749),与仅含截距模型相比的 JZS 贝叶斯因子显示支持备择假设的证据为极强(BF10 = 1.718e+9)。面积每增加 1 平方米,租金约增加 14.84(95% 可信区间 11.70 至 17.99);距市中心每远 1 千米,租金变化 −33.93(−51.63 至 −16.23)。在参考先验下,这些区间在数值上与经典置信区间一致,但可以直接按概率解读。

输出表格说明

模型摘要
N、R、R²、调整后 R²、估算标准误,以及模型相对仅含截距模型的 BF10/BF01。
ANOVA
最小二乘拟合得到的回归与残差平方和、自由度、均方、F 值及显著性。
系数的贝叶斯估计
截距、各协变量及各因子指示变量(命名为 因子=水平)的后验众数、平均值、方差及等尾可信区间。
误差方差的贝叶斯估计
σ² 的逆伽马后验,包括众数、平均值、方差及可信区间;图形显示每个系数及方差的后验。

前提假设与注意事项

  • 结果变量与预测变量呈线性关系,误差相互独立、近似正态且方差恒定。
  • 预测变量之间不存在完全共线性;奇异设计会直接报告而不予估计。
  • 因子采用指示符编码并以最后一个类别为参考,因此每个因子系数表示与该类别的差异。

结果如何验证

在含两个协变量和一个三水平因子的参考模型上,JZS 与 Rouder(s = 0.5)因子与 mpmath 30 位精度求积及 R 包 BayesFactor 0.9.12.4.8 的 lmBF 相对误差均在 1e-8 以内,Zellner 与超 g 因子分别在 1e-12 和 1e-8 以内;两种先验下系数后验在 1e-10 以内一致(可信区间端点 1e-9)。频率加权拟合与 BayesFactor 在复制行上的结果在 1e-8 以内一致。 验证证据(英文)

常见问题

可信区间与置信区间有区别吗?

在参考先验下,它们在数值上与经典 t 区间相等,但可直接解读为系数有 95% 的概率落在该区间内。共轭先验会使区间向先验平均值收缩。

为什么贝叶斯因子这么大?

拟合较好且样本量适中时,证据会迅速累积;大于 10 000 的值以科学计数法显示。通常大于 100 即称为极强证据。

它会单独检验每个预测变量吗?

不会。贝叶斯因子比较的是完整模型与仅含截距的模型。可根据各预测变量的可信区间是否包含零来判断,或分别运行两个模型并将其 BF10 相除来比较。

因子以哪个类别为参考?

排序后的最后一个类别,与 SPSS 一致。如需其他基线,请对变量重新编码。

我的数据会被上传吗?

不会。模型在浏览器的 Web Worker 中拟合,数据始终留在你的设备上。

用你自己的数据运行贝叶斯线性回归

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

打开分析工作区

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