使用马尔可夫链蒙特卡罗方法,使用贝叶斯方法对缺失数据的多元阈值自回归(MTAR)模型进行参数估计。执行MTAR过程的模拟(mtarsim())、矩阵参数和阈值的估计(mtarns()),使用贝叶斯变量选择(mtarstr())识别自回归阶数,使用Metropolised Carlin and Chib(mtarnumreg)识别状态数,并估计缺失数据,系数和协方差矩阵取决于自回归阶数、阈值和状态数(mtarmissing())。卡尔德隆和涅托(2017)<doi:10.1080/03610926.2014.990758>.