这个图是作者自己整理的

一眼看上去，你这个还是VAR模型的回归，因为有滞后期，开关变量就是新加入的变量呗。自己看看STR模型吧。

结果是作者自己整理成的表格

已成功 谢谢！

谢谢大神 ，已解决！能再问下你知道LSTAR参数怎么估计吗？我下面有什么错误吗？总是出现syntax error!
param c(1) 1 c(2) 1 c(3) 1  c(4) 1 c(5) 1 c(6) 1 c(7) 1 c(8) 1 c(9) 1 c(10) 1 c(11) 1 c(12) 1 c(13) 1 c(14) 1 c(15) 1 c(16) 1 c(17) 1 c(18) 1
nls dgdp=c(1)+c(2)*dgdp(-1)+c(3)*dgdp(-2)+c(4)*dgdp(-3)+c(5)*dgdp(-4)+c(6)*dgdp(-5)+c(7)*dral+c(8)*di+[c(9)+c(10)*dgdp(-1)+c(11)*dgdp(-2)+c(12)*dgdp(-3)+c(13)*dgdp(-4)+c(14)*dgdp(-5)+c(15)*dral+c(16)*di]*[1+exp(-c(17)*(dral3-c(18)))]^(-1)