Next: Design orderings based on Up: Statistical Properties Previous: Robustness properties   Contents

### Computational details

As has already been explained, the elements of statistical_properties are quantities which can be calculated starting from the information matrix . There are three fundamental calculations: the canonical variances, the pairwise variances, and the canonical efficiency factors.

The canonical variances are the inverses of the eigenvalues of , eigenvalues of zero corresponding to canonical variances of . Thus we need the roots of the polynomial . As is a rational matrix, this polynomial admits a factorization into irreducible factors over the rational field. Thus, in theory, the multiplicities of the canonical variances can be determined exactly, even if some of the values themselves are irrational. If the eigenvalues of are numerically extracted directly without factoring the characteristic polynomial, then the problem of inexact counts of those eigenvalues can arise.

Pairwise variances are defined above in terms of the Moore-Penrose inverse of : . In fact, any generalized inverse of can be used, from which Let be an all-ones matrix. If is connected, then is invertible for any and is a generalized inverse of (the same operation can be carried out for the connected components of if is disconnected). Thus pairwise variances can be calculated by inversion of a rational, nonsingular matrix.

Efficiency factors are defined as eigenvalues of the matrix , which can certainly be irrational. Extracting the roots of with respect to , that is, solving the equation , produces values for satisfying and otherwise . Thus efficiency factors can be found by extracting roots of a symmetric, rational matrix, involving the same computational issues as for the canonical variances.

The number of infinite canonical variances equals the number of connected components of less 1 (this being zero for any connected designs). Numerical extraction of eigenvalues of can potentially produce, at a given level of precision, values indistinguishable from zero that are in actuality positive, consequently producing an erroneous number of infinite canonical variances. This approximation error is prohibited by cross-checking against the connected indicator.

Next: Design orderings based on Up: Statistical Properties Previous: Robustness properties   Contents
Peter Dobcsanyi 2003-12-15