Let be iid with density , set , and put and . The finite-second-moment assumption gives, for ,
The sample variance consistently estimates . The central limit theorem and Slutsky theorem therefore give the asymptotic confidence interval
Its confidence probability is approximately , conventionally written ; the printed omits the factor . For , almost surely and the estimator has zero error. This interval concerns iid Monte Carlo estimators; correlated simulation output requires a long-run variance of a stationary process estimate instead of the iid sample variance.
Under fit.2, observations in each period are independent Poisson random variables with a common mean, so each sample variance should fluctuate around its sample mean. The plot instead has variances roughly 130–300 for means roughly 80–120, well above the line variance equals mean. This is overdispersion.
Day-to-day traffic variation can produce a Poisson mixture: a random daily intensity increases the marginal variance and can correlate observations from the same day. A useful replacement is a negative binomial regression, with variance , or a Poisson generalized linear mixed model with a day random intercept. A day-factor effect can also be included if inference is restricted to the observed days. A Quasi-Poisson regression estimates a dispersion multiplier for uncertainty, but does not itself explain the shared-day dependence and has no ordinary likelihood-based AIC. The means can retain the period-factor or selected jump structure.
An unbiased estimator of a parameter satisfies for every permissible value of . Since the sample mean is linear,
so is unbiased.
Use the identity
The two expectations on the right are and , respectively. Hence
Thus
The unbiased sample variance instead divides the same sum of squares by .
Because every column of has sample mean zero, every column of also has sample mean zero. The spectral theorem for real symmetric matrices gives
Hence for the sample covariance is , while the sample variance of is
Thus the are pairwise uncorrelated sample principal components, ordered by decreasing sample variance.
The given identities imply
Thus, writing
we have the exact identity
The sample variance satisfies , so
The Slutsky theorem and part (c) now give
Let the columns of a centered design matrix be variables and write by the spectral theorem for real symmetric matrices. The columns of are the sample principal-component score vectors. They satisfy , so distinct score vectors have zero sample covariance and has sample variance under the divisor- convention.