scipy.stats.linregress返回一个对应于斜率的p值,但不返回截取的p值。考虑文档中的以下示例:
>>> from scipy import stats
>>> import numpy as np
>>> x = np.random.random(10)
>>> y = np.random.random(10)
>>> slope, intercept, r_value, p_value, std_err = stats.linregress(x,y)
>>> p_value
0.40795314163864016
根据文档,p-value
是“假设检验的双侧p值,其零假设是斜率为零。”我希望得到相同的统计数据,但针对的是截距而不是斜率。
statsmodels.regression.linear_model.OLS返回开箱即用的两个系数的p值:
>>> import numpy as np
>>> import statsmodels.api as sm
>>> X = sm.add_constant(x)
>>> model = sm.OLS(y,X)
>>> results = model.fit()
>>> results.pvalues
array([ 0.00297559, 0.40795314])
仅使用scipy,我如何获得截取的p值(0.40795314163864016)?
发布于 2017-11-11 21:22:48
来自SciPy.org文档:https://docs.scipy.org/doc/scipy-.14.0/reference/generated/scipy.stats.linregress.html
在:打印“r平方:”,r_value**2
输出:R平方: 0.15286643777
对于其他参数,请尝试:
print ('Intercept is: ', (intercept))
print ('Slope is: ', (slope))
print ('R-Value is: ', (r_value))
print ('Std Error is: ', (std_err))
print ('p-value is: ', (p_value))
https://stackoverflow.com/questions/28753502
复制