For a balanced incomplete block design with , has eigenvalues on the zero-sum subspace and on the constant vector. Thus its rank is . As is , .
Let be the incidence matrix of a set system. The balanced incomplete block design overlap counts give . On the zero-sum subspace its eigenvalue is ; on the constant vector its eigenvalue is . Therefore has matrix rank . Since a product cannot have larger rank than either factor,
This proves Fisher's inequality for block designs directly. The incomplete, nontrivial hypotheses ensure the positive eigenvalue needed for this argument.
Count incidences of a treatment with a block containing it. Summing over blocks gives , while summing over treatments gives . Hence the balanced incomplete block design satisfies
A balanced incomplete block design has treatments and blocks, each containing distinct treatments with . Every treatment occurs in blocks and every distinct treatment pair occurs together in blocks. A symmetric balanced incomplete block design has , which implies by the first counting identity below.
For the separate even-order condition, let be the square treatment-by-block incidence matrix of a set system of a symmetric design. Its diagonal overlap counts are and its off-diagonal counts are , giving
The counting identities proved in the next slots show that and . Consequently the eigenvalues of are on the constant vector and on its -dimensional orthogonal complement. Taking determinants gives
If is even, the exponent is odd. The exponent of every prime in and in is even; therefore its exponent in the positive integer must be even. This proves the even-order symmetric design square obstruction:
It is a necessary condition, not a sufficiency claim.
A randomized complete block design groups similar experimental units into experimental blocks and places every treatment once in each experimental block, independently randomizing treatment allocation within experimental blocks. Blocking removes background between-block variation from within-block treatment comparisons. A balanced incomplete block design has treatments and experimental blocks of size , with each treatment occurring in experimental blocks and every distinct treatment pair together in experimental blocks. Counting treatment occurrences and co-occurrences gives and .
A Latin square of order places each of treatment symbols once in every row and every column of a array. This controls two nuisance directions. A Graeco-Latin square superposes two mutually orthogonal Latin squares: each square is Latin, and each ordered pair of symbols from the two squares occurs exactly once. The two symbol factors can then be fitted together with row and column effects. Random permutations of row labels, column labels and treatment symbols provide a suitable randomized allocation without destroying these balances.
For the operator experiment, fit the additive normal linear model
Here account for the two fabric directions, and is the operator effect. Additivity, independent equal-variance errors and the normal assumption supply the exact F-test; there is only one observation per cell, so an unrestricted interaction model is not available.
In the analysis of variance for a Latin square, the three centered factor spaces each have dimension . Together with the intercept the model has rank ten, leaving residual statistical degrees of freedom. The corrected total has . The completed calculation is
Under the null hypothesis , the operator F-statistic is
The upper five-percent critical value for is , so the operator effect is not significant at five percent in this model.
The decision is independent of the order of fitting rows, columns and operators. To prove the required orthogonality, let be any centered row coefficients and any centered operator coefficients. Their observation-space inner product is , since every operator occurs once in each row. The same calculation applies to columns and operators, and all row-column pairs occur once as well. Thus the three centered spaces are mutually orthogonal; their projections commute and sequential sums of squares equal adjusted sums of squares in every fitting order.
With machines included, the layout is a Graeco-Latin square: the machine symbols are Latin, and all sixteen operator-machine pairs occur exactly once. Fit
in addition to the earlier constraints. The machine space contributes three new statistical degrees of freedom and is orthogonal to the previous factor spaces, so the operator sum of squares remains . The residual sum of squares and statistical degrees of freedom become
Hence
This exceeds but not . Therefore the operator effect is now significant at five percent, but not at one percent. Adjustment for machines has removed substantial residual variation without altering the operator numerator. This illustrates why a balanced nuisance factor can change precision and significance even when it does not change the estimated treatment contrasts.
Write the binary incidence matrix of a set system as to avoid confusing it with a sample size. Then
counts the number of blocks containing treatment , and for , counts blocks containing both treatments. The displayed structure therefore makes every treatment replication equal to and every pair concurrence equal to .
Assume the usual nonempty block design, with and . If , any block containing one treatment contains every other treatment: otherwise that treatment's replication would exceed its concurrence with the omitted one. Every block is consequently complete, so and . With random allocation within blocks this is a randomized complete block design. If , the design is a balanced incomplete block design. Counting partners of a fixed treatment in its blocks gives
For , gives , hence . Moreover has eigenvalue on the -dimensional zero-sum subspace and eigenvalue on the constant direction. Thus , while . This proves , the Fisher's inequality for block designs. The case allows singleton blocks; their lack of connectedness prevents general treatment comparison after block adjustment, but does not invalidate the rank inequality.
First ignore the day information and use a one-way normal linear model with five group means and independent errors of common variance. There are observations. The analysis of variance has treatment, residual and corrected-total statistical degrees of freedom , respectively. The mean squares in ANOVA give
Under equal treatment means this has the F-distribution with degrees. It exceeds the given upper critical value , so the unblocked analysis rejects equality at .
The actual day allocation omits a different treatment on each of the five days. Hence it is a balanced incomplete block design with
Every pair occurs on exactly the three days omitting neither member. In the additive two-way analysis of variance, use an intercept, four independent treatment contrasts and four day contrasts. The design is connected, so the model rank is and its residual has statistical degrees of freedom. Day effects and treatment effects are not orthogonal, because a different treatment is absent each day; the unadjusted treatment sum of squares cannot be used for the adjusted test.
Enter days first and then treatments. The supplied sequential sums of squares give
The day, treatment-after-day, residual and corrected-total statistical degrees of freedom are . Under the additive normal linear model and the null of equal treatment effects, the adjusted numerator and residual squared lengths are independent chi-squared quantities, so the ratio has the F-distribution with degrees. It is below even the given critical value . After accounting for days there is no significant treatment difference at either or . This inference assumes independent homoscedastic normal errors and an additive day effect; with one observation per observed treatment–day cell, unrestricted interactions cannot also be fitted.
A balanced incomplete block design is symmetric when its number of blocks equals its number of treatments. The identity then gives . Its square incidence matrix of a set system obeys , with eigenvalues of multiplicity and of multiplicity one.