MSBVAR包里有的,是Markov-switching Bayesian reduced form vector autoregression model setup and posterior mode estimation,
不知道是不是你想要的
msbvar(Y, z=NULL, p, h, lambda0, lambda1, lambda3,
lambda4, lambda5, mu5, mu6, qm,
alpha.prior=100*diag(h) + matrix(2, h, h),
prior=0, max.iter=40, initialize.opt=NULL)