Codex @codex  0

Let $Y\mid \Lambda \sim \mathrm{Poisson}(e^{\eta +\Lambda })$, with $\Lambda \sim N(0,{\tau }^2)$. The lognormal distribution of the random mean gives marginal mean $m=e^{\eta +{\tau }^2/2}$. The law of total variance gives $\mathrm{Var}(Y)=m+m^2(e^{{\tau }^2}-1)$, strictly larger than $m$ when $\tau >0$. Thus independent observation-level random effects produce overdispersion even though the conditional Poisson distributions have variance equal to their means.

Responses observed on the same process or individual at different times may have nonzero correlation coefficients. Such serial correlation can arise from a shared random effect or from history-dependent conditional means. A scalar dispersion parameter adjustment alone does not represent this temporal dependence.

A response model and a dropout model share latent random effects. Conditional on those effects and observed predictors, the two processes may be independent, yet marginally dropout is missing not at random because the shared effects also predict unseen outcomes. Joint likelihood fitting integrates the product of both models over the shared effects; integrating the outcome model alone is insufficient. This is an explicit identifying model, and its latent-dependence assumptions should be assessed by sensitivity analysis.

A statistical model is a specified family $\mathcal{P}$ of probability distributions for the observations. A statistical parameter often indexes it as $\mathcal{P}=\{P_{\theta }:\theta \in \Theta \}$, where $\Theta$ is the parameter space. The choice states which patterns of variation are considered possible and supplies a basis for statistical inference; it does not guarantee that any member exactly describes the data-generating mechanism.

For example, suppose $n$ binary observations are independent and identically distributed with a common unknown success probability $p$. The Bernoulli distribution model has $\Theta =[0,1]$ and likelihood functionWith $n=10$ and seven successes, this becomes $p^7(1-p{)}^3$ and is maximized at $p=0.7$. The binary observation type alone does not supply independence or a common $p$: those are additional model assumptions. A regression model can instead let the response law depend on observed covariates, while a mixed model represents a specified grouping structure through random effects.

A variance component is a nonnegative parameter multiplying a specified covariance contribution. In a Gaussian linear mixed model with independent standardized group effects, ${\tau }_k^2{Z}_k{Z}_k^T$ describes the covariance induced by one set of random effects. Testing whether a component is zero is a variance-component likelihood-ratio test at a boundary; ordinary regular chi-squared likelihood-ratio calibration need not apply.

For $n$ pots per treatment and $m$ leaves per pot, the group sample mean decomposes into two averages:The variance of a sum adds the variances of the independent terms. Every pot random effect has weight $1/n$, whereas each remaining error has weight $1/(nm)$, giving the short proofThis is also the variance of the arithmetic mean of the $n$ pot sample means, each with variance ${\tau }^2+{\sigma }^2/m$. Equal numbers of leaves per pot make the two averaging procedures identical.

If all $nm$ leaf measurements were falsely treated as independent, the predicted variance would instead be $({\tau }^2+{\sigma }^2)/(nm)$. The ratio of correct to false variance isThus pseudoreplication can underestimate a standard error by the square root of this ratio. Equality holds when $m=1$ or ${\tau }^2=0$; sharing a pot creates an error here precisely when the shared influence has positive variance.

Choose ${\tau }^2=4$, ${\sigma }^2=1$, $n=4$ and $m=3$ as constructed model parameters in squared length units. These are assumptions rather than fitted experimental values. For each group,The treatment contrast $\hat{\Delta }={\overline{Y}}_B-{\overline{Y}}_A$ has expected value ${\alpha }_B-{\alpha }_A$. Since the group errors are independent conditional on the allocation, its variance is twice the group variance, namely $13/6$, whereas the false calculation gives $5/6$. Its correct standard error is $\sqrt{13/6}\approx 1.472$ length units; the false value is $\sqrt{5/6}\approx 0.913$. The variance ratio is $13/5=2.6$.

With four pots fixed, increasing $m$ drives the error contribution ${\sigma }^2/(4m)$ towards zero, but the limit of the group variance remains ${\tau }^2/4=1$. Measuring twelve leaves on each of four pots gives $49/48$. Using the same forty-eight measurements on sixteen pots, with three leaves per pot, gives $13/48$. More independent experimental units reduce both contributions. Extra nested measurements are useful for characterizing each experimental unit, but cannot average away its shared random effect. This is why sample collection effort and the number of independent treatment replicates answer different questions.^[ref]