用R来编写一个用来估计Logistic混合自回归模型的参数程序

   模型:X_t=S_t*P*X_t-1+(1-S_t)*X_t-1+E_t  ,其中E_t -N(0,a²),P为系数,S_t=0或1 且P(S_t=1|X_t-1)是logistic里面的形式.

我参考的文章是用EM算法来估计的,但由于在M步中有些参数要用迭代算法才能算,我就不知怎样来写了

求教各位大侠帮忙解决一下