我使用statsmodels(对其他python选项开放)来运行一些线性回归。我的问题是,我需要回归没有截距和限制系数的范围(0,1),也总和为1。
我试过这样的方法(至少是1的总和):
from statsmodels.formula.api import glm
import pandas as pd
df = pd.DataFrame({'revised_guess':[0.6], "self":[0.55], "alter_1":[0.45], "alter_2":[0.2],"alter_3":[0.8]})
mod = glm("revised_guess ~ self + alter_1 + alter_2 + alter_3 - 1", data=df)
res = mod.fit_constrained(["self + alter_1 + alter_2 + alter_3 = 1"],
start_params=[0.25,0.25,0.25,0.25])
res.summary()
但仍在努力实施“非负”系数约束。