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.
Continuous differentiability gives the first-order expansion
The assumed convergence in distribution implies , so after multiplication by the remainder is . The Slutsky theorem gives the multivariate delta method
For observations , the log-likelihood is, up to an additive constant,
Its unique maximum is
Now and, using , . The central limit theorem followed by the delta method for gives
For , the delta method with gives
Thus the transformation stabilizes variance at large .
A Gaussian linear model for treats the transformed observations as having additive constant-variance errors and models as linear. A Poisson regression with square-root link retains the Poisson likelihood and models
The two procedures therefore have different likelihoods, fitted means, and constraints; transforming the response is not the same as transforming its conditional mean.
The central limit theorem gives
Applying the delta method to yields
For this denotes the degenerate distribution at zero, which is the correct limit at the scale.
Let
Its gradient, in the coordinate order , is
Apply the Delta method to the limit in part (c). Multiplying the displayed covariance matrix there on both sides by this gradient gives
Consequently the asymptotic distribution of the Gaussian sample correlation is
The transformation satisfies and . Applying the delta method to part c gives
Thus exponentiating half the robust location estimate produces an asymptotically normal scale estimator.
The delta method gives
If , this denotes the degenerate distribution at zero.
Since and , the central limit theorem gives
Applying the delta method to yields
Thus, with , an asymptotic interval centered at the MLE is
optionally intersected with . Slutsky's theorem gives coverage tending to .
Put . The geometric moments are
For , one has and
The central limit theorem for followed by the delta method therefore gives
For transition-intensity vector and generator , the expected infectious occupancy over three months from is the Markov reward model quantity
Insert the fitted intensities to obtain , evaluating the matrix exponential and integral numerically if necessary.
If has estimated covariance , calculate the numerical gradient at . The delta method gives estimated variance and the approximate confidence interval
Simulation from followed by transformation through gives a useful asymmetric alternative.
Since , the rate MLE has asymptotic variance . The delta method for the logarithm gives
Independence of the arms therefore yields
Under , the asymptotic variance of the log rate ratio is