Codex @codex  0

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.

A homogeneity study with replicated test portions samples several containers of a reference material and measures separately prepared portions from each. The replication separates variation between containers from variation remaining within a container and the analytical procedure. Repeated detector readings of one prepared portion do not supply the same information.

The 2006 IUPAC protocol recommends randomly selecting at least ten distribution containers, homogenizing each selected container separately, and analyzing two portions from each in randomized order under a specified repeatability condition of measurement. It also requires adequate analytical precision and checks for trends, excessive rounding and discordant duplicate pairs.^[ref]

An original balanced random effect model isAssume independent zero-mean container effects with variance ${\sigma }_b^2$ and independent equal-variance residuals with variance ${\sigma }_e^2$. Residuals include analytical scatter and whatever within-container variation remains at the chosen portion size. One-way analysis of variance then givesso a method-of-moments estimate is ${\hat{\sigma }}_b^2=(M{S}_B-M{S}_W)/r$. A negative estimate can occur through sampling variation; it is not negative physical variance or proof of perfect reference material homogeneity.

For a small illustrative calculation, take duplicate pairs $(9,11)$, $(11,13)$ and $(13,15)$ in one measurement unit. The container means are 10, 12 and 14. The within-container sum of squares is 6 on three degrees of freedom, so $M{S}_W=2$; the between-container sum is 16 on two degrees of freedom, so $M{S}_B=8$. Thus ${\hat{\sigma }}_b^2=(8-2)/2=3$. These three invented containers explain the arithmetic; they do not meet the protocol's recommended ten-container study.

Averaging $r$ independent portions from one container gives variance ${\sigma }_b^2+{\sigma }_e^2/r$. More replication reduces the residual term without removing the common container effect. Randomizing the measurement order helps prevent instrumental drift from masquerading as between-container variation.

Fitness for purpose requires comparing material variation with an appropriate tolerance, not merely failing to reject zero variation. The protocol explicitly distinguishes sufficient homogeneity from the stricter hypothesis that sampling variance is exactly zero.^[ref] Keep the portion size, sampling design and uncertainty with the conclusion.

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.

A nested wringing-repeatability study separates variation between independently prepared gauge block wringings from readout scatter within one unchanged wring. Repeating a detector reading many times does not sample the contact variability introduced by separating and rejoining the surfaces.

For $m$ wrings with $r$ readings each, an original balanced random effect model iswhere independent wring effects have variance ${\tau }^2$ and independent readout errors have variance ${\sigma }^2$. Assume stable conditions and independence between those levels. One-way analysis of variance givesThus $M{S}_W$ estimates readout variance and $(M{S}_B-M{S}_W)/r$ estimates the between-wring variance in this model. A negative variance estimate reflects sampling variability or model limitations, not a physically negative contact variance.

For a small invented data set, duplicate readings in nm relative to an arbitrary reference are $(-3,-1)$, $(1,3)$ and $(5,7)$. Wring means are $-2,2,6$; $M{S}_W=2$ and $M{S}_B=32$, so the method-of-moments between-wring variance is $15\mathrm{n}{\mathrm{m}}^2$. These are arithmetic examples, not a recommended sample size or a NIST calibration result.

For the overall mean of the balanced study,With stipulated $\tau =6\mathrm{nm}$, $\sigma =2\mathrm{nm}$, five wrings and four readings per wring, its standard deviation is $\sqrt{7.4}\approx 2.72\mathrm{nm}$. Increasing to forty readings on the same five wrings lowers it only to about 2.69 nm. Fifty wrings with four readings each would instead give about 0.86 nm, if the same model and stable preparation remain valid.

A shared reference or persistent preparation contribution does not disappear under either averaging rule. An independent common uncertainty $u_0$ adds $u_0^2$ to the displayed variance. This connects the nested study with the wringing-film covariance model and the gauge-block reference uncertainty floor.

NIST's handbook documents contact-layer dependence on preparation and time, including evolution after thermal equilibration.^[ref] Record rewringing events, timestamps and surface preparation; use appropriate controls for drift and damage. A nested contact-study allocation should balance new wrings against repeat readings, rather than treating every reading as an independent contact experiment.

Use the account of one-way analysis of variance in part a of the Experimental Psychology paper. The procedure, its assumptions and the scope of the omnibus conclusion are the same in this variant.

Use one-way analysis of variance to compare a quantitative response variable across levels of one categorical factor; for example, compare task accuracy across three independently assigned teaching methods. The ordinary independent-groups procedure tests the null hypothesis that all group population means are equal. Its F-test compares variation associated with the group factor against within-group error variation; a significant result indicates that not all population means are equal, without identifying which groups differ.^[ref]

The classical one-way normal linear model assumes independent observations and within-group errors approximately following a normal distribution with a common variance. Repeated observations on the same people require a different dependence structure; pronounced unequal variances may call for Welch's one-way analysis of variance. Planned contrasts or appropriately controlled follow-up comparisons identify scientifically relevant differences. Report the means, differences and confidence intervals, not just significance. Causal interpretation depends on the design of experiments, particularly randomization, rather than on the F-test alone.

A heteroscedastic extension of independent-groups one-way analysis of variance, allowing unequal group variances. It does not make paired or clustered observations independent.