Why Work on the Log Scale?
A parameter such as a rate, variance, or mean is often known to be positive. A conventional Wald interval adds and subtracts a margin on the parameter’s original scale, so its lower endpoint can fall below zero. That is not merely inconvenient: a negative endpoint cannot represent a possible value of a positive parameter.
When an estimator is positive, taking its logarithm offers another route. The logarithm maps positive values to the whole real line, where a normal approximation can be used to form a familiar symmetric interval. Exponentiating the two endpoints returns an interval on the original scale. That interval has positive endpoints and is generally asymmetric around the estimate.
The earlier material on delta-method intervals provides the main approximation behind this construction: a smooth transformation of an asymptotically normal estimator is often asymptotically normal as well. Here the transformation is \(g(\theta)=\log\theta\), and the focus is the resulting interval and its interpretation.
Construction and Conditions
Fix a confidence level \(1-\alpha\), where \(0<\alpha<1\), and let \(z_{1-\alpha/2}\) be the \(1-\alpha/2\) quantile of the standard normal distribution. Suppose an estimated standard error \(\widehat{\operatorname{se}}_{\log}\) is available for \(\log\hat{\theta}\). The ordinary two-sided interval on the log scale is
Exponentiating its lower and upper endpoints gives the interval
The estimated standard error is on the log scale. If the standard error is instead first derived for \(\hat{\theta}\), the delta method often gives the approximation \(\widehat{\operatorname{se}}_{\log}\approx\widehat{\operatorname{se}}(\hat{\theta})/\hat{\theta}\). This follows because the derivative of \(\log\theta\) is \(1/\theta\), evaluated at the estimate for a plug-in standard error. A model may also provide a direct standard error for the log estimator.
For the approximation to justify a confidence interval, the log estimator must have an appropriate large-sample normal approximation, and the estimated standard error must consistently estimate its scale. The following result states one convenient version of these requirements.
Proof. On the event \(\hat{\theta}>0\), the logarithm is strictly increasing. Taking logarithms of the two endpoint inequalities shows that the coverage event is equivalent to
Since \(\widehat{\operatorname{se}}_{\log}>0\) with probability tending to one under the stated consistency condition, divide by it on that event. Equivalently, the resulting event is
By the assumed convergence in distribution and Slutsky’s theorem, the fraction converges in distribution to \(N(0,1)\). The standard normal distribution assigns probability \(1-\alpha\) to the interval \([-z_{1-\alpha/2},z_{1-\alpha/2}]\), and has probability zero at either endpoint. Therefore the probability of the displayed event converges to \(1-\alpha\). The event \(\hat{\theta}\leq0\) has probability tending to zero, so it does not change the limit. \(\square\)
What the Endpoints Mean
Write \(c=z_{1-\alpha/2}\widehat{\operatorname{se}}_{\log}\), so \(c\geq0\). The endpoints are \(\hat{\theta}e^{-c}\) and \(\hat{\theta}e^c\). Both are positive whenever \(\hat{\theta}>0\). Their geometric mean is exactly the estimate:
Thus the estimate is the center on the log scale, or equivalently the geometric center on the original scale. The endpoints are not generally equally far from \(\hat{\theta}\) in ordinary units. When \(c>0\), the upper distance is \(\hat{\theta}(e^c-1)\), while the lower distance is \(\hat{\theta}(1-e^{-c})\). Their difference is
The inequality follows because \(e^c+e^{-c}>2\) for \(c>0\). Hence the interval is strictly asymmetric around \(\hat{\theta}\) when \(c>0\). When \(c=0\), the endpoints coincide with the estimate, so there is no strict asymmetry. The multiplicative interpretation is direct: the upper endpoint is \(\hat{\theta}\) multiplied by \(e^c\), and the lower endpoint is \(\hat{\theta}\) multiplied by \(e^{-c}\).
Worked Example: An Exponential Rate
Suppose \(X_1,\ldots,X_n\) are independent exponential observations with rate \(\lambda>0\), so \(E[X_i]=1/\lambda\). For \(n=100\), suppose the observed sum is \(\sum_{i=1}^{100}X_i=50\). Construct an approximate 95% log-scale interval for \(\lambda\).
The likelihood is proportional to \(\lambda^n\exp(-\lambda\sum_i X_i)\), so its maximizer is \(\hat{\lambda}=n/\sum_iX_i=100/50=2\). For this model the Fisher information in one observation is \(1/\lambda^2\), so the asymptotic standard error of \(\hat{\lambda}\) is approximately \(\lambda/\sqrt n\). Applying the delta method to \(\log\hat{\lambda}\), whose derivative is \(1/\lambda\), gives an asymptotic standard error \(1/\sqrt n\). Here use \(\widehat{\operatorname{se}}_{\log}=1/\sqrt{100}=0.1\).
With \(z_{0.975}=1.96\), the log-scale margin is \(c=1.96(0.1)=0.196\). The endpoints are
The approximate 95% interval is \([1.644,2.433]\). Both endpoints are positive; the estimate \(2\) is the geometric, not arithmetic, center of the interval.
Worked Example: A Poisson Rate
Suppose a count \(K\) observed over exposure \(T\) has a Poisson distribution with mean \(\lambda T\), where the rate \(\lambda>0\) is unknown. If \(K=25\) events occur over \(T=10\) exposure units, construct a 95% log-scale interval using the large-count approximation.
The rate estimate is \(\hat{\lambda}=K/T=25/10=2.5\). Since the variance of a Poisson count is its mean, the approximate variance of \(\hat{\lambda}\) is \(\lambda/T\). The delta method for the logarithm gives an approximate variance of \((1/\lambda^2)(\lambda/T)=1/(\lambda T)\). Substituting \(\hat{\lambda}\) gives \(\widehat{\operatorname{se}}_{\log}=1/\sqrt{\hat{\lambda}T}=1/\sqrt K=1/5=0.2\).
The log-scale margin is \(c=1.96(0.2)=0.392\). Therefore the interval is
This is a large-sample approximation, not an exact Poisson interval. The calculation illustrates how the log standard error depends on the observed count and why the count must be positive to use this particular plug-in formula.
Worked Example: A Positive Mean
Five positive measurements are \(2,2,4,6,6\). Treat them as independent observations from a population with positive mean \(\mu\) and finite variance. Construct a 95% log-scale interval using the usual sample standard error as a plug-in estimate.
The sample mean is \(\bar X=(2+2+4+6+6)/5=4\). The deviations from the mean are \(-2,-2,0,2,2\), whose squared values sum to \(16\). Thus the sample variance is \(s^2=16/(5-1)=4\), and \(s=2\). The central limit theorem gives an approximate standard error \(s/\sqrt{5}\) for \(\bar X\). Applying the delta method to \(\log\bar X\) gives
The margin is \(c=1.96(0.2236)\approx0.4383\), so the endpoints are
The resulting approximate interval is \([2.581,6.200]\). Because this illustration has only five observations, the large-sample approximation may be inaccurate; the calculation demonstrates the construction, not a guarantee of good finite-sample coverage.
Practical Checks and a Common Mistake
A log-scale interval is useful when positive values span a range where relative differences are more meaningful than absolute differences. For example, multiplying a rate estimate by 1.2 represents the same proportional change whether the original rate is small or large. The interval’s endpoints reflect that multiplicative scale.
Before reporting one, check that the estimated parameter is positive for the observed data, that the standard error corresponds to the log scale, and that the normal approximation is credible for the sample size and model. If a chosen estimator can be zero in finite samples, the logarithm is not defined at those samples; this construction cannot be computed there without changing the estimator or method.
- Common mistake: adding and subtracting the log-scale margin from the original estimate. That produces \(\hat{\theta}\pm c\), which is not the interval derived from normal inference for \(\log\theta\). The correction is to exponentiate the log endpoints, yielding \(\hat{\theta}e^{-c}\) and \(\hat{\theta}e^c\).
- Common mistake: using the standard error of \(\hat{\theta}\) as if it were the standard error of \(\log\hat{\theta}\). The delta-method plug-in approximation is \(\widehat{\operatorname{se}}(\hat{\theta})/\hat{\theta}\), not simply \(\widehat{\operatorname{se}}(\hat{\theta})\).
- Common mistake: claiming strict asymmetry whenever the interval is exponentiated. The margin satisfies \(c\geq0\). The interval is strictly asymmetric around the estimate when \(c>0\); if \(c=0\), both endpoints equal the estimate.
Check Your Understanding
Use these questions to check both the construction and its interpretation.
- Write the endpoints of a 95% log-scale interval when \(\hat{\theta}=3\) and \(\widehat{\operatorname{se}}_{\log}=0.15\). Are they symmetric around 3 in ordinary units?
- Why does exponentiating the endpoints preserve the coverage event when the parameter is positive?
- If an estimate has standard error \(0.4\) on its original scale and its estimated value is \(2\), what is the delta-method plug-in standard error for its logarithm?
- What happens to the interval endpoints when the estimated log-scale standard error equals zero?
- Why is a log-scale interval not directly computable from an observed estimate equal to zero?