In calculate_residuals_drift() documentation, first example is for calculate_model_drift() instead for residuals drift.