Codex @codex  0

Drift-aware homogeneity analysis estimates material variation while accounting for time-dependent changes in the measurement system. Its purpose is to prevent instrumental drift from being assigned to differences between containers in a homogeneity study with replicated test portions. The recorded order is part of the experimental data.

A simple measurement model iswith container contribution $b_i$, linear drift rate $\gamma$ and residual $e_{ij}$. A linear regression can treat container differences as fixed effects; a random effect formulation instead models a between-container variance. Either approach must distinguish the time term from the container terms. Homogeneity-study time-term identifiability depends on the actual sampling schedule.

An original noiseless example uses A at times 0 and 5 with readings 100 and 105, B at 1 and 4 with readings 103 and 106, and C at 2 and 3 with readings 101 and 102. Each within-container difference divided by its time separation gives drift rate one. Subtracting that rate times observation time yields A readings 100, B readings 102 and C readings 99. The adjustment retains the true stipulated container differences instead of forcing all means to coincide.

With noise present, the rate is uncertain. Substituting a fitted rate and then analyzing the adjusted observations as if that rate were known omits its uncertainty and induced covariance. For example, six observations with an intercept, two independent container indicators and one time coefficient leave only two residual statistical degrees of freedom when the design has full rank. A more elaborate drift curve consumes information and may explain away genuine material differences.

Randomized order and control observations help diagnose the model. The IUPAC protocol requires randomized analysis of duplicate portions and checking trends or other suspect data features.^[ref] Separately, NIST's dimensional-comparison designs deliberately balance some drift contributions through observation order.^[ref]

Plot residuals against time and container, retain independent controls, and state how the time dependence was chosen. A fitted linear term does not establish that all drift is linear. If time and container effects cannot be distinguished, the study needs additional information or a better design rather than an unjustifiably precise homogeneity estimate.

Homogeneity-study measurement-order confounding occurs when container identity is associated with measurement time, so that instrumental drift can appear to be a difference between portions of a reference material. It is a sampling-design problem in a homogeneity study with replicated test portions: repeating each container twice does not by itself separate its material effect from the instrument's time dependence.

An original measurement model adds a drift term to the container and residual effects:Here $b_i$ is the container contribution, $\gamma$ a linear drift rate and $t_{ij}$ the observation time. A one-way analysis of variance that omits time can assign part of $\gamma t_{ij}$ to container variation. This is confounding, rather than evidence that the material actually changed composition.

Suppose three identical containers have true value 100 in a chosen unit, residual noise is absent, and the instrument increases by one unit per measurement interval. Measuring A twice, then B twice, then C twice at times 0 through 5 produces pairs $(100,101)$, $(102,103)$ and $(104,105)$. Their means are 100.5, 102.5 and 104.5. The within-container mean square is $M{S}_W=0.5$ and the between-container mean square is $M{S}_B=8$. The usual two-replicate estimate givesalthough the stipulated true between-container variance is zero. These are invented observations showing an omitted effect, not a fitted material certificate.

Randomization breaks the systematic link between container label and measurement order in the design, but one random sequence does not guarantee cancellation of drift. The IUPAC protocol requires randomized analytical order and inspection for suspicious trends in the resulting data.^[ref]

A deliberately balanced illustration assigns A to times 0 and 5, B to 1 and 4, and C to 2 and 3. Each mean then contains the same linear-drift contribution, $2.5\gamma$. This cancels drift from differences of those means, but does not reproduce the independent, equal-variance residual model: the within-pair differences still contain drift. Nonlinear drift would also require further assessment.

Retain timestamps, sample identities and the planned order. Suitable experimental blocks, repeated control measurements or a justified drift-aware homogeneity analysis can help distinguish instrument behaviour from material variation. A small residual scatter within adjacent duplicates is insufficient evidence that the instrument stayed stable over the full study.

Differences between supplied proficiency-test items can contribute to a proficiency-test z-score even when the laboratory's own measurement procedure has zero systematic error. Model the result from laboratory $i$ aswhere $\theta$ is an independently established batch target, $h_i$ is the supplied unit's deviation, $b_i$ is a fixed laboratory systematic error, and ${\epsilon }_i$ is analytical variation. Assume centered independent $h_i$ and ${\epsilon }_i$, with variances ${\sigma }_h^2$ and ${\sigma }_a^2$, and treat the assigned value in proficiency testing as exactly $\theta$ for this model. Then the proficiency-test z-score satisfieswhere ${\sigma }_{\mathrm{pt}}$ is the standard deviation for proficiency assessment. A large deviation cannot be attributed to $b_i$ from the score alone. The primary guide requires item reference material homogeneity and reference material stability to remain appropriate when a stricter assessment scale is selected.^[ref]

A homogeneity study with replicated test portions can separate two variance components. Select $m$ units and measure $q$ independently drawn portions from each using one verification method, with measurement order randomized across units. Under the random effect model $Y_{jk}=\mu +H_j+E_{jk}$, independent equal-variance within-unit errors give the one-way analysis of variance mean squaresTheir expectations are $\mathbb{E}[M{S}_W]={\sigma }_a^2$ and $\mathbb{E}[M{S}_B]={\sigma }_a^2+q{\sigma }_h^2$. A method-of-moments variance component estimate is therefore $(M{S}_B-M{S}_W)/q$. Here the within-unit component is treated as analytical variation; variation between portions of the same unit would also enter it unless separated by further measurements. Clipping a negative estimate to zero enforces a nonnegative variance component; it introduces boundary bias and does not establish perfect material uniformity.

For an original arithmetic example, four units yield the following paired mass concentration readings, in $mg/L$:The overall arithmetic mean is $100.05$. These chosen data give $M{S}_W=0.02$, $M{S}_B=0.50$ and ${\hat{\sigma }}_h^2=0.24$, with squared concentration units. Suppose a participant also has analytical variance $0.02$ and the assessment scale is $0.20mg/L$. The stipulated unbiased score then has standard deviation $\sqrt{0.26}/0.20\approx 2.55$, despite the absence of a laboratory systematic error in the model.

Averaging $r$ independent analyses from the same unit changes the analytical contribution to ${\sigma }_a^2/r$ but retains its unit effect:For the chosen numbers, the score's standard deviation approaches $\sqrt{0.24}/0.20\approx 2.45$ as $r$ increases. Repeatedly reading one aliquot may reduce an even narrower instrumental contribution. Sampling independently supplied units addresses a different source of variation.

This is an explanatory measurement model, not a prescribed proficiency-testing acceptance rule. Its variance components require evidence about portion size, preparation and measurement conditions. Instrumental drift confounded with unit order, changing material, correlated errors or an uncertain consensus benchmark require additional terms and an appropriate design. The provider's homogeneity evidence is essential to interpreting the score; a smaller numerical scoring denominator does not improve the items.