handles, labels = ax1.get_legend_handles_labels()
handles = handles[:-len(handles) // 3][1::2]
labels = labels[:-len(handles) // 3][1::2]
ax1.legend(handles=handles, labels=labels, loc='best', shadow=True, numpoints=2
    )
plt.show()
Difficulty setting a legend with statsmodels Python package