The Four Types of Estimable Functions

General Form of an Estimable Function

This section demonstrates a shorthand technique for displaying the generating set for any estimable bold upper L. Suppose

bold upper X equals Start 6 By 4 Matrix 1st Row 1st Column 1 2nd Column 1 3rd Column 0 4th Column 0 2nd Row 1st Column 1 2nd Column 1 3rd Column 0 4th Column 0 3rd Row 1st Column 1 2nd Column 0 3rd Column 1 4th Column 0 4th Row 1st Column 1 2nd Column 0 3rd Column 1 4th Column 0 5th Row 1st Column 1 2nd Column 0 3rd Column 0 4th Column 1 6th Row 1st Column 1 2nd Column 0 3rd Column 0 4th Column 1 EndMatrix and bold-italic beta equals Start 4 By 1 Matrix 1st Row  mu 2nd Row  upper A 1 3rd Row  upper A 2 4th Row  upper A 3 EndMatrix

bold upper X is a generating set for bold upper L, but so is the smaller set

bold upper X Superscript bold asterisk Baseline equals Start 3 By 4 Matrix 1st Row 1st Column 1 2nd Column 1 3rd Column 0 4th Column 0 2nd Row 1st Column 1 2nd Column 0 3rd Column 1 4th Column 0 3rd Row 1st Column 1 2nd Column 0 3rd Column 0 4th Column 1 EndMatrix

bold upper X Superscript bold asterisk is formed from bold upper X by deleting duplicate rows.

Since all estimable bold upper Ls must be linear functions of the rows of bold upper X Superscript bold asterisk for bold upper L bold-italic beta to be estimable, an bold upper L for a single-degree-of-freedom estimate can be represented symbolically as

upper L Baseline italic 1 times left-parenthesis 1 1 0 0 right-parenthesis plus upper L Baseline italic 2 times left-parenthesis 1 0 1 0 right-parenthesis plus upper L Baseline italic 3 times left-parenthesis 1 0 0 1 right-parenthesis

or

bold upper L equals left-parenthesis upper L Baseline italic 1 plus upper L Baseline italic 2 plus upper L Baseline italic 3 comma upper L Baseline italic 1 comma upper L Baseline italic 2 comma upper L Baseline italic 3 right-parenthesis

For this example, bold upper L bold-italic beta is estimable if and only if the first element of bold upper L is equal to the sum of the other elements of bold upper L or if

bold upper L bold-italic beta equals left-parenthesis upper L Baseline italic 1 plus upper L Baseline italic 2 plus upper L Baseline italic 3 right-parenthesis times mu plus upper L Baseline italic 1 times upper A 1 plus upper L Baseline italic 2 times upper A 2 plus upper L Baseline italic 3 times upper A 3

is estimable for any values of L1, L2, and L3.

If other generating sets for bold upper L are represented symbolically, the symbolic notation looks different. However, the inherent nature of the rules is the same. For example, if row operations are performed on bold upper X Superscript bold asterisk to produce an identity matrix in the first 3 times 3 submatrix of the resulting matrix

bold upper X Superscript asterisk bold asterisk Baseline equals Start 3 By 4 Matrix 1st Row 1st Column 1 2nd Column 0 3rd Column 0 4th Column 1 2nd Row 1st Column 0 2nd Column 1 3rd Column 0 4th Column negative 1 3rd Row 1st Column 0 2nd Column 0 3rd Column 1 4th Column negative 1 EndMatrix

then bold upper X Superscript asterisk bold asterisk is also a generating set for bold upper L. An estimable bold upper L generated from bold upper X Superscript asterisk bold asterisk can be represented symbolically as

bold upper L equals left-parenthesis upper L Baseline italic 1 comma upper L Baseline italic 2 comma upper L Baseline italic 3 comma upper L Baseline italic 1 minus upper L Baseline italic 2 minus upper L Baseline italic 3 right-parenthesis

Note that, again, the first element of bold upper L is equal to the sum of the other elements.

With multiple generating sets available, the question arises as to which one is the best to represent bold upper L symbolically. Clearly, a generating set containing a minimum of rows (of full row rank) and a maximum of zero elements is desirable.

The generalized g 2-inverse left-parenthesis bold upper X prime bold upper X right-parenthesis Superscript minus of bold upper X prime bold upper X computed by the modified sweep operation (Goodnight 1979) has the property that left-parenthesis bold upper X prime bold upper X right-parenthesis Superscript minus Baseline bold upper X prime bold upper X usually contains numerous zeros. For this reason, in PROC GLM the nonzero rows of left-parenthesis bold upper X prime bold upper X right-parenthesis Superscript minus Baseline bold upper X prime bold upper X are used to represent bold upper L symbolically.

If the generating set represented symbolically is of full row rank, the number of symbols left-parenthesis upper L Baseline italic 1 comma upper L Baseline italic 2 comma ellipsis right-parenthesis represents the maximum rank of any testable hypothesis (in other words, the maximum number of linearly independent rows for any bold upper L matrix that can be constructed). By letting each symbol in turn take on the value of 1 while the others are set to 0, the original generating set can be reconstructed.

Last updated: December 09, 2022