能告诉你什么
Crosstabs 在你的设备上运行与 SPSS BAYES REGRESSION 对应的过程,支持连续协变量和指示符编码的因子、四种 g 先验贝叶斯因子,以及参考先验或共轭先验。
- 数据在多大程度上支持完整模型优于仅含截距的模型(BF10),并同时给出 R、R² 和调整后 R²。
- 截距和每个系数的后验众数、后验平均值、后验方差及可信区间。
- 误差方差的后验分布,直接说明还有多少未被解释的离散程度。
- 通过共轭先验,了解关于系数大小的先验信息如何改变估计。
适用场景
适合使用
- 用多个连续或分类预测变量预测一个连续结果,例如根据面积、距离和城区预测租金。
- 希望系数区间可以按概率解读,或希望纳入先验估计。
- 需要证明一组预测变量没有作用(BF01 较大),而不只是 F 检验不显著。
改用其他方法
- 二分类或计数结果:使用二元 Logistic 回归或广义线性模型。
- 需要每个预测变量的贝叶斯因子或模型平均:面板只报告整个模型的一个贝叶斯因子。
- 逐步选择变量、残差诊断或共线性统计:使用线性回归。
在 Crosstabs 中如何操作
- 01
打开数据
打开文件,切换到“统计”模式,在“贝叶斯统计”组中选择“贝叶斯线性回归”。
- 02
设置模型
选择“因变量(标度)”,在“协变量(标度)”中勾选连续预测变量,在“因子(指示符编码,以最后一个类别为参考)”中勾选分类预测变量。
- 03
选择贝叶斯因子
在“贝叶斯因子”中选择“JZS(Zellner–Siow)”“Zellner(固定 g,留空则 g = N)”、带参数 a 的“超 g”,或带尺度 s 的“Rouder”。
- 04
选择系数的先验
保留“参考先验,p(θ, σ²) ∝ 1/σ²”,或选择“共轭正态–逆伽马”,为所有系数设置“形状 a0”“尺度 b0”以及共同的“先验平均值”和“先验方差(× σ²)”。
- 05
运行
设置“可信区间(%)”,声明工作区权重(如有),然后点击“运行贝叶斯线性回归”。缺失值按列表删除。
计算示例
以下结果由 Crosstabs 引擎在构建本页时根据模拟数据实际计算得出,并非手工填写。
数据
40 套模拟公寓,包括面积(平方米)、到市中心的距离(千米)和月租金,由带噪声的线性模型以固定随机种子生成。
| N | R | R² | 调整后 R² | 估算标准误 | BF10 |
|---|---|---|---|---|---|
| 40 | 0.865 | 0.749 | 0.735 | 154.76 | 1.718e+9 |
JZS 方法(Zellner–Siow 柯西先验),模型与仅含截距的模型进行比较。
| 参数 | 后验众数 | 后验平均值 | 后验方差 | 95% 可信区间下限 | 95% 可信区间上限 |
|---|---|---|---|---|---|
| (Intercept) | 491.67 | 491.67 | 13766.945 | 260.45 | 722.90 |
| area | 14.84 | 14.84 | 2.550 | 11.70 | 17.99 |
| distance | −33.93 | −33.93 | 80.698 | −51.63 | −16.23 |
参考先验 p(θ, σ²) ∝ 1/σ²;等尾可信区间。
| 参数 | 后验众数 | 后验平均值 | 95% 可信区间下限 | 95% 可信区间上限 |
|---|---|---|---|---|
| 误差方差 | 22723.2 | 25320.1 | 15919.5 | 40089.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 中拟合,数据始终留在你的设备上。
相关过程
- 线性回归REGRESSION
- 贝叶斯单因素方差分析BAYES ANOVA
- 贝叶斯 Pearson 相关BAYES CORRELATION
- 权重估计(加权最小二乘)WLS
- GLM 单变量UNIANOVA
- 免费的 SPSS 替代工具:支持范围
用你自己的数据运行贝叶斯线性回归
打开 .sav、CSV 或 Excel 文件。所有计算都在浏览器中完成,文件不会离开你的设备。
SPSS 是 IBM Corp. 的商标。Crosstabs 与 IBM 无关联,也未获得 IBM 认可。