我正试图用scipy.stats.linregress来计算两组数据之间的最小平方意义上的比例因子。然而,尽管输入xi
变量是向量而不是n X 2
矩阵,它还是给了我一个截距。
所以,一个简单的代码如下:
from scipy import stats
from numpy import arrange,array
y = [0, 11, 19, 28, 41, 49, 62, 75, 81]
xi = arange(0,9)
scale, intercept, r_value, p_value, std_err = stats.linregress(xi,y)
运行这个,我得到的刻度是10.383,但我也得到了-0.86的截距。我怎么能告诉它只适合缩放参数,截距应该保持为零。
如果您想要适合一个模型,
y~xi
而不需要截取,您可能需要考虑使用更多面向统计的包,例如statsmodels
:您可以使用
R
独立验证结果。只是现在必须显式地将截距指定为0:引擎盖下的计算很简单:
相关问题 更多 >
编程相关推荐