This function first computes Lord's CSEM under the binomial model via
csem_lord(ni), where ni = ncol(x). It then rescales the
resulting CSEM curve using the ratio
$$
\sqrt{\frac{1 - \text{KR-20}}{1 - \text{KR-21}}},
$$
where KR-20 and KR-21 are computed from the observed data via
kr20(x) and kr21(x), respectively.