4# qinzi 你好!我调试了一下你的程序,我猜可能的问题是:
z是指具有结构突变的变量,你在程序中将具有结构突变的变量设为4,也许理论上来讲是可以的(对于这个问题我也不太清楚,我没有做过关于结构突变的研究),但是根据我推测,即使理论上可行,也应该需要比较复杂的计算,对参数设置的要求也会比较复杂。而在实际操作中,一般同一个模型中不会同时这么多具有结构突变的变量(当然不排除特殊情况),一般情况是一次设置1个,其他的变量则设为x, 若存在多个可能具有突变的变量,你可以根据每个变量的经济意义,选取你认为可能会出现结构突变的变量,逐个进行估计。
如在你的程序中,可将q设为1,p设为3,将除z(1维)以外的3个变量都设为x,应该可以得到正确的估计结果,你可以试一下。另外,修改eps1和m也会对结果产生影响,不妨修改试试看。
个人意见,供你参考。