在使用从今年 1 月开始的数据的疫苗每日计数时,我使用 SARIMAX 预测未来几天的计数。我使用以下代码:
import statsmodels.api as sm
model=sm.tsa.statespace.SARIMAX(vax_Data_india1['Daily_count'],order=(1, 1, 1),seasonal_order=(1,1,1,20))
results=model.fit()
vax_Data_india1['forecast']=results.predict(start=71,end=83,dynamic=True)
vax_Data_india1[['Daily_count','forecast']].plot(figsize=(12,8))