使用scipy的stats.linregress和一个固定的y截距

Question

我想知道是否可以使用stats.linregress，除了固定的y-intercept，我目前正在这样做：

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)

如果我仍然能够确定一个R值

，那将是更好的选择

Answer 1

从scipy domcumentation开始，这是不可能的。但您可以使用其他内容，例如statsmodel：

y = [0, 11, 19, 28, 41, 49, 62, 75, 81]
x = arange(0,9)
model = sm.OLS(y,x)
results = model.fit()
print(results.params)