在Python中，如何拟合最小分位数b样条回归线？

import numpy as np
import statsmodels.api as sm

# 生成示例数据
np.random.seed(0)
x = np.linspace(0, 10, 100)
y = 2 * np.sin(x) + np.random.normal(size=100)

# 添加常数列
X = sm.add_constant(x)

# 拟合最小分位数b样条回归线
quant_reg = sm.QuantReg(y, X)
res = quant_reg.fit(q=0.5)  # 这里的q表示分位数，0.5表示中位数

# 打印回归结果
print(res.summary())

# 绘制拟合曲线
import matplotlib.pyplot as plt
plt.scatter(x, y, alpha=0.5)
plt.plot(x, res.predict(X), 'r')
plt.show()