Any way you can use a `scipy.stats` function for this? Now that I think of it, we should really include an inverse two-way within-subjects p-value to f-value function (survival function or 1-cdf) that takes 1) the number of conditions in each way, 2) the number of subjects, 3) the p-value, and 4) the condition of interest or conditions of interest (list) and returns the corresponding f-value threshold. We have one for the t-test, we should have one for this.