作Yt-m*Yt-1=a(Xt-m*Xt-1)+b+随机干扰项  的回归，在R里怎么写程序，Xt  Yt序列都知道的。