如题所述，小弟做一个Vasicek模型的仿真，使用shibor数据，求出模型的参数，参考文献中对方程使用最小二乘法计算：
r(t+1)=a+br(t)+epsilon(t+1),epsilon为误差项。我对R不熟，不知如何计算。求助大家