啊我图怎么没了

模型设定和做出的图如上

调用第三方包

图是脉冲响应函数（IRF），显示了一个变量（如贸易融资）对另一个变量冲击（如有无融资）的响应。用R、Stata或Python等软件估计IVAR模型。
import statsmodels.api as sm
from statsmodels.tsa.api import VAR
# 假设一个包含相关时间序列的DataFrame 'df'
model = VAR(df)
results = model.fit(4)  # 选择滞后期为4，根据AIC或BIC调整
模型估计后，计算脉冲响应函数（IRF），即跟踪一个变量对另一个变量冲击的反应。
irf = results.irf(10)  # 计算10期的脉冲响应
irf.plot(orth=False)  # 绘制IRF图
诊断测试残差（如自相关、异方差性）

可以麻烦具体一点吗