zr = smoothn(y, robust=True)  # Robust smoothing