The distribution function is on , so its median is . The likelihood function is
which decreases over its admissible range. Hence . By the invariance property of maximum likelihood estimation,
By the invariance property of maximum likelihood estimation,
The sampling distribution is , where
has a log-normal distribution with
Writing , its exact fractional bias and variance are
Hence, to leading order,
The maximum-likelihood estimators are the success proportions
The invariance property of maximum likelihood estimation therefore gives
The invariance property of maximum likelihood estimation gives
Since , the ratio has a log-normal distribution. Its exact variance is
and the lowest-order delta method approximation is
For fixed , case (i) has , whereas case (ii) has . Because and is the same in both cases, independent supernova-level variation in case (ii) gives the smaller fractional variance: averaging reduces it, while a shared galaxy fluctuation does not average away.