研究HOG与CORN之间的非对称传导关系，用冲量持续门限自回归(MomentumConsistent Threshold Autoregressive,MC-TAR)模型拟合：ecmAsyFit(y=HOG, x=CORN,lag=4, model="mtar",split=FALSE,thresh=-0.420)。
要接着作正向冲击、负向冲击的脉冲响应对比图，如下图Figure 3:

用R语言如何实现，即用哪个程序包的哪个函数命令实现？