Calculate mean reversion second order differential equation

Question

I would like to calculate mean reversion speeds from second order reactions

    y <-weekly_spread_data

    y.lag <- stats::lag(y, -1)
    delta.y <-diff