相关系数在线计算:Pearson、Spearman 与 Kendall
Pearson、Spearman 与 Kendall 三种相关系数一次算出, 给出各自的 p 值、Pearson 的 95% 置信区间与 r²。
三个系数在测不同的东西
Pearson r 测的是线性关系的强度。它用的是原始数值, 所以对离群点很敏感,也只能捕捉直线关系。
Spearman ρ 先把数据换成秩,再算 Pearson。它测的是单调关系—— 只要「一个涨另一个也涨」就算,不要求成直线。
Kendall τ 直接数「同向的数对」与「反向的数对」, 比 Spearman 更稳健,小样本时分布性质更好,但计算量随 n 平方增长。
有并列值时本工具报的是 τb。 同向数减反向数(记作 S)除以什么有两种取法: 除以全部数对 n₀ = n(n−1)/2 得 τa,除以 √[(n₀−n₁)(n₀−n₂)] 得 τb, 其中 n₁、n₂ 分别是 x、y 里并列造成的数对数。 没有并列时两者完全相同;一旦有并列,τa 会系统性偏小—— 例如 x = 1,2,2,3,5,5,5 与 y = 2,3,4,5,5,6,8,τa 给 0.762 而 τb 给 0.868。 R、scipy、SPSS 报的 Kendall τ 都是 τb,所以本工具也用 τb,输出里会写成「τb」提醒你。 p 值同样做了并列校正(Kendall 1945 的方差公式),否则会偏保守。
两个例子说明为什么要一起看
例一,单调但非线性。 x = 1…8,y = 2ˣ: Pearson 给出 0.850,Spearman 与 Kendall 都给出 1.000。 关系是完美单调的,Pearson 却因为它不是直线而低估了。
例二,一个离群点足以让符号反过来。 x = 1,2,…,10 而 y = 10,9,8,7,6,5,4,3,2,50: 除最后一点外 y 随 x 严格递减,但那一个 50 让 Pearson 给出 +0.358,Spearman 给出 −0.455——两者符号相反。 只看 Pearson 会得出完全错误的结论。
所以:三个一起看,差别大就说明数据有故事。 Pearson 远小于 Spearman 通常意味着关系是单调但弯曲的;两者符号相反几乎总是离群点造成的。
r² 的正确读法
r² 是「线性关系能解释的方差比例」,不是「正确率」。 r = 0.7 看起来很高,但 r² = 0.49,意味着还有一半的变异没被解释。 r = 0.5 时 r² 只有 0.25。相关系数的平方关系让人容易高估关联强度。
相关不等于因果,而且这句话还不够
除了「相关不蕴含因果」之外,还有两件事同样重要:
相关系数为 0 不等于无关。 它只说明没有线性(或单调)关系。 U 形关系的 Pearson r 可能正好是 0,但两个变量显然有关。
样本量大时,毫无实际意义的相关也会显著。 n = 1000 时 r = 0.07 就能给出 p < 0.05, 而它解释的方差只有 0.5%。看 p 值之前先看 r 和 r²。
置信区间
Pearson 的区间用 Fisher z 变换算:先把 r 变换到近似正态的尺度, 在那里算对称区间,再变换回来。所以区间关于 r 不对称, r 越接近 ±1 越明显。n < 10 时区间会很宽,这是诚实的反映,不是缺陷。
常见问题
三个系数该报哪一个?
看数据。关系明显是直线且没有极端值,报 Pearson;关系单调但弯曲,或有离群点,报 Spearman 或 Kendall。样本量小(n < 20)时 Kendall 的分布性质更好。关键是:方法要在看到结果之前定好,不能三个都算完挑最好看的那个报。
Pearson 和 Spearman 差很多说明什么?
Spearman 明显更高,通常说明关系是单调但弯曲的(比如指数或对数关系),Pearson 只测直线所以低估了。两者符号相反几乎总是离群点造成的——个别极端值能把 Pearson 拉到相反方向,而秩相关不受影响。这两种情况本页都会主动提示。
r = 0 能说明两个变量无关吗?
不能。它只说明没有线性关系(Spearman 为 0 则是没有单调关系)。U 形关系的 Pearson r 可能正好是 0,但两个变量显然有关。任何相关分析之前都该先画散点图——这是相关系数无法替代的一步。
为什么置信区间关于 r 不对称?
因为用的是 Fisher z 变换:先把 r 变换到近似正态的尺度,在那里算对称区间,再变换回来。r 越接近 ±1,不对称越明显。这是正确行为,不是计算错误。
样本量大时 p 值很小,是不是就说明关联强?
不是。n = 1000 时 r = 0.07 就能给出 p < 0.05,而它解释的方差只有 0.5%。p 值回答「这个关联像不像随机」,r 和 r² 才回答「关联有多强」。先看后者。
相关解读
- 2026-10-26
VitA缺乏肠道富集了口咽菌群,L-丝氨酸补充却让鞘脂病小鼠更糟
三篇文章分别用横断面菌群数据、鞘脂机制研究和可行性RCT,呈现了营养代谢领域同一道坎:发现了中间指标的变化,却还确认不了因果方向。L-丝氨酸在SPTSSA突变小鼠中的逆向效应尤为警示——遗传背景可以把同一种营养素的效果从无害变成致命。
- 2026-10-23
双歧杆菌宿主适应机制与肠道菌群在脑肿瘤和关节炎研究中的新进展
三篇文献从机制到临床勾勒肠道菌群的影响边界:双歧杆菌通过多糖降解和免疫调控影响宿主神经内分泌;肠脑轴连接菌群与胶质瘤,但临床转化有限;一项跨四国队列协议将菌群纳入银屑病关节炎的生物标志物采集框架。机制跑在证据前面。
- 2026-10-22
原花青素与婴儿双歧杆菌改善酒精戒断焦虑抑郁,三组干预效果相当
小鼠AUD戒断模型中,原花青素(50 mg/kg)与婴儿双歧杆菌(1×10⁷ CFU)单用或联用均能改善焦虑抑郁样行为,机制涉及肠屏障修复、外周Treg/Th17再平衡、中枢炎症下调和前额叶皮质髓鞘蛋白上调。三组干预效果相当,联用无额外增益,全部为关联证据。
相关工具
统计检验怎么选:按数据类型、组数与配对关系对照
回答三个问题就能定下用哪种检验:比数值还是比比例、有几组、是否配对。每一行直接链到对应工具。
米氏方程 Km 与 Vmax 计算(非线性拟合 + 三种线性化对照)
输入若干组底物浓度与反应速率,用非线性拟合求 Km 与 Vmax,并把 Lineweaver-Burk、Hanes-Woolf、Eadie-Hofstee 三种线性化的结果并排列出——它们的分歧本身就是数据质量的读数。
IC50 / EC50 四参数拟合计算器(含数据可信度检查)
输入若干组浓度与响应,用四参数 logistic 拟合 IC50/EC50、Hill 斜率与两端平台,并检查浓度有没有铺到平台、固定斜率会让结果偏多少——这两件事比 R² 重要。