stats.linregress 中的 r 与 statsmodels 中的 r-squared 相比答案

【问题标题】：r in stats.linregress compared to r-squared in statsmodelsstats.linregress 中的 r 与 statsmodels 中的 r-squared 相比
【发布时间】：2018-08-08 04:40:41
【问题描述】：

我正在开发一个程序来研究某些类星体的星等和红移之间的相关性，我正在使用statsmodels 和scipy.stats.linregress 来计算数据的统计数据； statsmodels 计算 r-squared（以及其他参数），stats.linregress 计算 r（以及其他参数）。

一些示例输出是：

W1 r-squared: 0.855715
W1 r-value  : 0.414026
W2 r-squared: 0.861169
W2 r-value  : 0.517381
W3 r-squared: 0.874051
W3 r-value  : 0.418523
W4 r-squared: 0.856747
W4 r-value  : 0.294094
Visual minus WISE r-squared: 0.87366
Visual minus WISE r-value  : -0.521463

我的问题是，为什么 r 和 r-squared 的值不匹配

（即对于 W1 频段，0.414026**2 != 0.855715）？

我的计算函数代码如下：

def computeStats(x, y, yName):
    from scipy import stats
    import statsmodels.api as sm

    #   Compute model parameters
    model = sm.OLS(y, x, missing= 'drop')
    results = model.fit()
    #   Mask NaN values in both axes
    mask = ~np.isnan(y) & ~np.isnan(x)
    #   Compute fit parameters
    params = stats.linregress(x[mask], y[mask])
    fit = params[0]*x + params[1]
    fitEquation = '$(%s)=(%.4g \pm %.4g) \\times redshift+%.4g$'%(yName,
                params[0],  #   slope
                params[4],  #   stderr in slope
                params[1])  #   y-intercept

    print('%s r-squared: %g'%(name, arrayresults.rsquared))
    print('%s r-value  : %g'%(name, arrayparams[2]))

    return results, params, fit, fitEquation

我是否错误地解释了统计数据？还是这两个模块使用不同的方法计算回归？

【问题讨论】：

标签： python-3.x scipy statsmodels

【解决方案1】：

默认情况下，statsmodels 中的OLS 不包括线性方程中的常数项（即截距）。（常数项对应于design matrix 中的一列。）

要匹配linregress，请像这样创建model：

    model = sm.OLS(y, sm.add_constant(x), missing= 'drop')

【讨论】：