计量模型，共有4个变量，分别为gap  lnpgdp  （lnpgdp）^2    sbzc ,四个变量均为一阶单整，最优滞后阶数为2阶，存在1个协整关系，VECM模型构建后是稳定的系统，想请教如何继续进行脉冲响应分析和方差分解，有大佬可以提供stata 的运行命令吗？