Fieller's theorem ($(\hat{C}-r\hat{E}{)}^2\le q^2(v_C-2r{v}_{CE}+r^2{v}_E)$)

For jointly approximately normal estimates $(\hat{C},\hat{E})$ with estimated variances $v_C,v_E$ and covariance $v_{CE}$, a confidence interval for the ratio $C/E$, where $E\ne 0$, can be obtained by inverting tests of $C-rE=0$. For each proposed ratio $r$, the contrast has estimate $\hat{C}-r\hat{E}$ and variance $v_C-2r{v}_{CE}+r^2{v}_E$. Retain $r$ when its absolute standardized contrast is at most the chosen critical value $q$. This gives the displayed quadratic inequality, with coefficients $A={\hat{E}}^2-q^2{v}_E$, $B=-2(\hat{C}\hat{E}-q^2{v}_{CE})$, $D={\hat{C}}^2-q^2{v}_C$. When $A>0$ and the discriminant is positive, retain the interval between the roots. When $A<0$, a confidence set can instead be two unbounded rays or the entire real line. Degenerate linear cases are interpreted directly from the inequality. A bounded interval must not be forced when the denominator is poorly separated from zero. Known covariance matrix and a joint normal distribution give exact normal critical values; estimated covariance usually gives approximate coverage, unless an appropriate exact common-scale pivot exists.