线性回归计算器

最小二乘 ŷ=a+bx:斜率/截距/r/R²、残差逐点表、r 显著性 t 检验、散点回归图与预测区间(x₀ 反推)

在线线性回归计算器:粘贴 X:Y 或「X,Y」两列数据(逗号 / 冒号 / 空格 / 分号宽容解析,每行一点), 最小二乘拟合 ŷ = a + bx(b = Sxy/Sxx、a = ȳ − b·x̄),输出斜率、截距、相关系数 r、决定系数 R²、 斜率标准误与 r 显著性 t 检验(t = r√(n−2)/√(1−r²),双侧 p 值基于不完全贝塔精确计算), 逐点残差表(yᵢ、ŷᵢ、残差)与 SVG 散点 + 回归线图; 预测模式输入 x₀ 给 ŷ₀ ± t·s·√(1 + 1/n + (x₀−x̄)²/Sxx) 预测区间并提示外推风险。 Anscombe I 组经典数据(a=3.0001、b=0.5001、r=0.8164)等锚点精确复现。全部计算在浏览器本地完成,不上传任何数据。

R² 高 ≠ 回归合理:先画散点图看形状再谈拟合——Anscombe 四组数据 R² 全是 0.67,但只有第一组真的是直线关系;预测区间只在观测范围内可信,外推需自担风险。

支持「10:8.04」「10,8.04」「10 8.04」「10;8.04」四种分隔,每行取前两个数为 (X, Y);空行与无法解析的行自动跳过。

最小二乘拟合结果(n = 11)

ŷ = 3.0001 + 0.5001·x

斜率 b(SE 0.1179)

0.5001

截距 a

3.0001

相关系数 r

0.8164

决定系数 R²

0.6665

残差标准差 s = 1.2366;b 的 t 统计量 = 4.2415(df = 9), 双侧 p = 0.00217 (** 高度显著)

r 的显著性检验

t = 4.2415

t = r√(n−2)/√(1−r²),df = n−2 = 9

Sxx / Syy / Sxy

110.000 / 41.273

55.010

数据范围

x ∈ [4.000, 14.000]

y ∈ [4.260, 10.840]

散点图与回归线

11.373.4014.603.73

残差逐点表(残差 = yᵢ − ŷᵢ)

#xᵢyᵢŷᵢ残差
110.00008.04008.00100.0390
28.00006.95007.0008-0.0508
313.00007.58009.5013-1.9213
49.00008.81007.50091.3091
511.00008.33008.5011-0.1711
614.00009.960010.0014-0.0414
76.00007.24006.00061.2394
84.00004.26005.0005-0.7405
912.000010.84009.00121.8388
107.00004.82006.5007-1.6807
115.00005.68005.50050.1795

|残差| > 2s 的点标红(可能的离群点);残差之和恒为 0,残差应大致随机散布——若呈弯曲/喇叭形,线性模型不合适。

预测模式(x₀ → ŷ₀ ± 预测区间)

预测区间 = ŷ₀ ± t(0.975, n−2)·s·√(1 + 1/n + (x₀−x̄)²/Sxx), 单个新观测的不确定度(含残差方差);如需「均值响应的置信区间」把根号里的 1 去掉即可(更窄)。

x₀ = 10.0000 的预测

ŷ₀ = 8.0010

95% 预测区间 [ 5.0671, 10.9349 ](SE_pred = 1.2970)

r 与 R² 怎么读

r ∈ [−1, 1] 线性相关方向与强度

R² = r² y 的方差被 x 线性解释的比例

b = r·(sy/sx) 斜率与 r 的关系

R² = 0.67 意思是「y 变异的 67% 能被 x 的线性项解释」。r 只度量线性关系: 强曲线关系(如 y=x²,对称区间)r 可以是 0;一个离群点也能把 r 从 0 拉到 0.9。

R² 的常见误读

  • 相关 ≠ 因果:冰淇淋销量与溺水人数 r 很高,共同原因是气温。
  • R² 低 ≠ 没关系:可能需要加变量、变换(取对数)或用非线性模型。
  • 范围敏感:截取 x 的一段会让 R² 骤变,跨研究比较 R² 要先对齐范围。
#线性回归公式、FAQ 与参考

线性回归计算器说明

最小二乘法找直线 ŷ = a + bx,使残差平方和 Σ(yᵢ−ŷᵢ)² 最小:b = Sxy/Sxx、a = ȳ − b·x̄, 其中 Sxy = Σ(xᵢ−x̄)(yᵢ−ȳ)、Sxx = Σ(xᵢ−x̄)²。r = Sxy/√(Sxx·Syy) 度量线性相关, R² = r² 解释变异比例;r 的显著性用 t = r√(n−2)/√(1−r²)(df = n−2)检验,p 值经不完全贝塔精确计算。 预测区间 SE_pred = s√(1 + 1/n + (x₀−x̄)²/Sxx),x₀ 离均值越远区间越宽——这是回归「两端翘起」的几何来源。 全部计算在浏览器本地完成。

公式汇总

斜率:b = Sxy/Sxx;截距:a = ȳ − b·x̄

相关系数:r = Sxy/√(Sxx·Syy);R² = r²

残差方差:s² = SSE/(n−2),SSE = Syy − b·Sxy = Syy(1−r²)

斜率标准误:SE_b = s/√Sxx;t = b/SE_b = r√(n−2)/√(1−r²)

预测区间:ŷ₀ ± t(0.975, n−2)·s·√(1 + 1/n + (x₀−x̄)²/Sxx)

R² 的含义与常见误读

R² 回答「x 的线性项解释了 y 变异的百分之多少」,它不回答「x 是否导致 y」「模型形式是否正确」, 也不保证预测可靠。四大误读:把相关当因果;R² 高就免检残差(离群点与曲线关系都能造出高 R²); 用 R² 比较不同因变量的模型(不可比);对时间序列数据直接回归(伪回归,r 可以虚高)。 配合本页残差表与散点图一起看,才是完整流程。

常见问题 (FAQ)

r=0.8 强相关吗?

看领域:物理测量里 0.8 很弱,社会科学里算强。更稳的口径是看 R²(0.64 的变异被解释)与散点图形态;样本量小时 |r| 还要除以抽样波动——本页 t 检验给出精确 p 值。

预测区间和置信区间有什么区别?

置信区间描述「x₀ 处平均响应 ŷ」的不确定度(只含参数估计误差,根号里没有 1);预测区间描述「单个新观测 y」的不确定度(还要加一项残差方差 s²,根号里的 1)。所以预测区间恒比置信区间宽,且不随 n→∞ 收窄到零。

数据里有离群点怎么办?

残差表已把 |残差| > 2s 的点标红。先查数据录入错误与实验异常,确认是真实观测后再考虑:稳健回归、加虚拟变量或分别报告「含/不含」两套结果——直接删除需要说明理由。

Anscombe 数据是什么?

F. J. Anscombe (1973) 构造的四组 x,y 数据:均值、回归线(ŷ=3+0.5x)、R²=0.67 几乎完全相同,但散点形态天差地别(直线、曲线、离群点、水平诱饵)。它证明只看回归数字不看图会犯系统错误——本页示例即第一组。

计算数据会上传吗?

不会。数据解析、回归拟合、残差与预测全部在你的浏览器本地计算,不发送任何网络请求,输入数据不被存储或上传。

参考来源

  • Legendre, A. M. (1805) 最小二乘法;Gauss (1809) 正规方程
  • Anscombe, F. J. (1973). Graphs in statistical analysis. The American Statistician 27.
  • r 显著性检验与预测区间:NIST/SEMATECH e-Handbook §4.1.4「Linear Regression」
  • t 分布 p 值经不完全贝塔函数(Lentz 连分式,Numerical Recipes §6.4)计算

相关工具推荐

单位换算器

支持长度、重量、温度、面积、体积等多种单位类型的相互转换

热门 本地

科学计算器

高级数学运算,支持三角函数、对数、指数运算等科学计算,包含内存管理和历史记录

热门 本地

统计学计算器

数据分析工具,支持描述统计、回归分析、概率分布、假设检验等统计计算

热门 本地

方程求解器

一次方程、二次方程、方程组求解,支持图像绘制、牛顿迭代法和数值方法

热门 本地

财务比率计算器

盈利能力、偿债能力、运营能力比率分析,支持图表展示、行业对比和财务健康评分

热门 本地

折扣计算器

计算折扣后价格、节省金额和折扣率,支持折扣率和折数输入

热门 本地

这张卡挂在哪些职业工具箱

相关教程