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.
With independent Gaussian errors, maximizing the likelihood over minimizes the residual sum of squares, so the maximum-likelihood estimators are the same as least squares. The invariance property of maximum likelihood estimation gives
The two estimates have zero covariance by part (ii). Therefore
This is the uncertainty in the estimated expected response. Estimating an additional noisy observation would add its own error variance, but the quantity requested is the mean response .