Hi, suppose I have a time series rt. Why does this code not work? rt <- data n <- length(rt) - 2 at <- vector(length = n) at = rt-(3.75101e-04 + 9.61877e-02 * rt[-1] - 2.48071e-02 * rt[-2]) Thank you, Benjamin