Add upper/lower limits based on the r2 value

limit_omega0_by_r2(fit_input, factor = 1.5)