Python——分位回归

import statsmodels.api as sm
import statsmodels.formula.api as smf
import matplotlib.pyplot as plt
import numpy as np
import pandas as pd
data=sm.datasets.engel.load_pandas().data
#描述统计
fig=plt.figure()
ax=fig.add_subplot(1,1,1)
box=plt.boxplot((data[‘foodexp’],data[‘income’]),
notch=True, patch_artist=True,labels=[‘foodexp’,‘income’])
colors = [‘lightblue’,‘pink’]
for patch, color in zip(box[‘boxes’], colors):
patch.set_facecolor(color)
patch.set_alpha(1)
plt.show()
#分位回归与线性回归
quantiles=np.arange(.05,.9