如题！查资料发现可以使用 statsmodels 库中的函数，但是官方文档的例子是一个一元的。想问下多元咋做？

a = np.random.random(20).reshape(10, 2)
y = np.dot(a, np.array([1, 0.5]).reshape(2, 1))
dat = np.hstack((a, y))
datt = pd.DataFrame(dat, columns = ['x1', 'x2', 'y'])
mod = smf.quantreg('y ~ x1 + x2', datt).fit(q = 0.1)

上面这样确实可以做二元的。但是如果协变量很多，这个公式'y ~ x1 + x2'就没法写成这种形式了
请教，有没有那种像普通最小二乘的，直接输入一个设计矩阵，和反应变量就可以的那种？

你给的这个范例是用的 statsmodels.formula.api 下的函数，接收公式表达式作为参数。

可以看看 statsmodels.api 下的相关函数，直接使用数据集作为输入参数。

    G_will 嗯嗯，谢谢！！我查到了。想问下有没有那种介绍统计模型，统计方法的python参考书啊？