The function f.phiinv() appears to introduce numerical imprecision on the scale of +/- 1e-15. This can introduce values for gamma slightly beyond the interval [0,1], which can be problematic for subsequent analysis. Rounding to the nearest 13th decimal place solves the issue.
The function f.phiinv() appears to introduce numerical imprecision on the scale of +/- 1e-15. This can introduce values for gamma slightly beyond the interval [0,1], which can be problematic for subsequent analysis. Rounding to the nearest 13th decimal place solves the issue.