需要帮助重新调整python fill_between与数据点

alert ！在这里问这样一个入门级的问题，我感到很尴尬！

嘿，伙计们，我一直在做一个项目，涉及时间序列总Electron 内容数据.我的目标是应用统计分析，找出地震引起的TEC异常.

我跟着this research paper for sliding IQR method.但由于某些原因，我没有得到如本文中所示的结果与给定的公式.所以我决定使用滚动平均值+1.6 STD方法.

问题是，当我使用ax. fill_between(x1 = index，y1 = ub，y2 = lb)方法时，我的置信区间带绘制得比数据点更远一步.请参见所给图，以便更好地理解.

以下是我目前正在做的事情:

df = DAEJ.copy()
window = 12
hourly = df.resample(rule="h").median()
hourly["ma"] = hourly["TEC"].rolling(window=window).mean()
hourly["hour"] = hourly.index.hour
hourly["std_err"] = hourly["TEC"].rolling(window=window).std()
hourly["ub"] = hourly["ma"] + (1.67* hourly["std_err"])
hourly["lb"] = hourly["ma"] - (1.67* hourly["std_err"])
hourly["sig2"] = hourly["TEC"].rolling(window=window).var()
hourly["kur"] = hourly["TEC"].rolling(window=window).kurt()
hourly["pctChange"] = hourly.TEC.pct_change(12, fill_method="bfill")
hourly = hourly.dropna()
dTEC = hourly[(hourly["TEC"] > hourly["ub"])]



fig, ax = plt.subplots(figsize=(12,4))
hourly["TEC"].plot(ax=ax, title="TEC Anomaly", label="Station: DAEJ")
ax.fill_between(x=hourly.index, y1=hourly["ub"], y2=hourly["lb"], color="red", label="Conf Interval", alpha=.4)
ax.legend()

这是我得到的结果！

如图所示，数据和色带没有正确对齐，我知道在计算滚动平均值与窗口为12小时将导致12小时偏移值，但即使在下降第一个值后，我仍然没有得到对齐的数字.