我有如下所示的代碼，我試圖在 7 月 - 12 月顯示的繪圖上為 df“l(fā)tyc”中的數(shù)據(jù)畫一條黑線。我的錯(cuò)誤是在 plt.legend 行之前的最后.import warningsimport itertoolsimport numpy as npimport matplotlib.pyplot as pltwarnings.filterwarnings("ignore")plt.style.use('fivethirtyeight')import pandas as pdimport statsmodels.api as smimport matplotlibmatplotlib.rcParams['axes.labelsize'] = 14matplotlib.rcParams['xtick.labelsize'] = 12matplotlib.rcParams['ytick.labelsize'] = 12matplotlib.rcParams['text.color'] = 'k'from sys import exitdf = pd.read_excel("MOSDailyWindSpeed.xlsx")wspdBH1 = df.groupby('Date')[' Simulated WS BH1PI'].sum().reset_index()wspdHOO = df.groupby('Date')[' Simulated WS HOO801'].sum().reset_index()wspdBH1 = wspdBH1.set_index('Date')wspdHOO = wspdHOO.set_index('Date')wspdBH1.indexwspdHOO.indexy = wspdHOO[' Simulated WS HOO801'].resample('MS').mean()#monthly mean --> change site here 'MS' is month start for 'Date' coly['2017':]#look at daily data starting 2017 -view data y.plot(figsize=(15, 6))plt.show()from pylab import rcParamsrcParams['figure.figsize'] = 18, 8decomposition = sm.tsa.seasonal_decompose(y, model='additive')fig = decomposition.plot()plt.show()p = d = q = range(0, 2)pdq = list(itertools.product(p, d, q))seasonal_pdq = [(x[0], x[1], x[2], 12) for x in list(itertools.product(p, d, q))]print('Examples of parameter combinations for Seasonal ARIMA...')print('SARIMAX: {} x {}'.format(pdq[1], seasonal_pdq[1]))print('SARIMAX: {} x {}'.format(pdq[1], seasonal_pdq[2]))print('SARIMAX: {} x {}'.format(pdq[2], seasonal_pdq[3]))print('SARIMAX: {} x {}'.format(pdq[2], seasonal_pdq[4]))for param in pdq:  for param_seasonal in seasonal_pdq:    try:        mod = sm.tsa.statespace.SARIMAX(y,                                        order=param,                                        seasonal_order=param_seasonal,                                        enforce_stationarity=False,                                        enforce_invertibility=False)