The Four Types of Estimable Functions

Examples

A One-Way Classification Model

For the model

upper Y equals mu plus upper A Subscript i Baseline plus epsilon i equals 1 comma 2 comma 3

the general form of estimable functions bold upper L bold-italic beta is (from the previous example)

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

Thus,

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

Tests involving only the parameters upper A 1, upper A 2, and upper A 3 must have an bold upper L of the form

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

Since this bold upper L for the A parameters involves only two symbols, hypotheses with at most two degrees of freedom can be constructed. For example, letting left-parenthesis upper L Baseline 2 comma upper L Baseline 3 right-parenthesis be left-parenthesis 1 comma 0 right-parenthesis and left-parenthesis 0 comma 1 right-parenthesis, respectively, yields

bold upper L equals Start 2 By 4 Matrix 1st Row 1st Column 0 2nd Column 1 3rd Column 0 4th Column negative 1 2nd Row 1st Column 0 2nd Column 0 3rd Column 1 4th Column negative 1 EndMatrix

The preceding bold upper L can be used to test the hypothesis that upper A 1 equals upper A 2 equals upper A 3. For this example, any bold upper L with two linearly independent rows with column 1 equal to zero produces the same sum of squares. For example, a joint test for linear and quadratic effects of A

bold upper L equals Start 2 By 4 Matrix 1st Row 1st Column 0 2nd Column 1 3rd Column 0 4th Column negative 1 2nd Row 1st Column 0 2nd Column 1 3rd Column negative 2 4th Column 1 EndMatrix

gives the same SS. In fact, for any bold upper L of full row rank and any nonsingular matrix bold upper K of conformable dimensions,

SS left-parenthesis upper H 0 colon bold upper L bold-italic beta equals 0 right-parenthesis equals SS left-parenthesis upper H 0 colon bold upper K bold upper L bold-italic beta equals 0 right-parenthesis

A Three-Factor Main-Effects Model

Consider a three-factor main-effects model involving the CLASS variables A, B, and C, as shown in Table 1.

Table 1: Three-Factor Main-Effects Model

Obs A B C
1 1 2 1
2 1 1 2
3 2 1 3
4 2 2 2
5 2 2 2


The general form of an estimable function is shown in Table 2.

Table 2: General Form of an Estimable Function for Three-Factor Main-Effects Model

Parameter Coefficient
mu (Intercept) L1
A1 L2
A2 L1L2
B1 L4
B2 L1L4
C1 L6
C2 L1 + L2L4 – 2 times L6
C3 L2 + L4 + L6


Since only four symbols (L1, L2, L4, and L6) are involved, any testable hypothesis will have at most four degrees of freedom. If you form an bold upper L matrix with four linearly independent rows according to the preceding rules, then testing bold upper L bold-italic beta equals bold 0 is equivalent to testing that normal upper E left-bracket bold upper Y right-bracket is uniformly 0. Symbolically,

SS left-parenthesis upper H 0 colon bold upper L bold-italic beta equals 0 right-parenthesis equals upper R left-parenthesis mu comma upper A comma upper B comma upper C right-parenthesis

In a main-effects model, the usual hypothesis of interest for a main effect is the equality of all the parameters. In this example, it is not possible to unambiguously test such a hypothesis because of confounding: any test for the equality of the parameters for any one of A, B, or C will necessarily involve the parameters for the other two effects. One way to proceed is to construct a maximum rank hypothesis (MRH) involving only the parameters of the main effect in question. This can be done using the general form of estimable functions. Note the following:

  • To get an MRH involving only the parameters of A, the coefficients of bold upper L associated with mu, B1, B2, C1, C2, and C3 must be equated to zero. Starting at the top of the general form, let L1 = 0, then L4 = 0, then L6 = 0. If C2 and C3 are not to be involved, then L2 must also be zero. Thus, A1A2 is not estimable; that is, the MRH involving only the A parameters has zero rank and upper R left-parenthesis upper A vertical-bar mu comma upper B comma upper C right-parenthesis equals 0.

  • To obtain the MRH involving only the B parameters, let L1 = L2 = L6 = 0. But then to remove C2 and C3 from the comparison, L4 must also be set to 0. Thus, B1B2 is not estimable and upper R left-parenthesis upper B vertical-bar mu comma upper A comma upper C right-parenthesis equals 0.

  • To obtain the MRH involving only the C parameters, let L1 = L2 = L4 =0. Thus, the MRH involving only C parameters is

    upper C Baseline italic 1 minus 2 times upper C Baseline italic 2 plus upper C Baseline italic 3 equals upper K left-parenthesis for any upper K right-parenthesis

    or any multiple of the left-hand side equal to K. Furthermore,

    SS left-parenthesis upper H 0 colon upper C Baseline italic 1 minus 2 times upper C Baseline italic 2 plus upper C Baseline italic 3 equals 0 right-parenthesis equals upper R left-parenthesis upper C vertical-bar mu comma upper A comma upper B right-parenthesis

A Multiple Regression Model

Suppose

normal upper E left-bracket upper Y right-bracket equals beta 0 plus beta 1 x 1 plus beta 2 x 2 plus beta 3 x 3

where the bold upper X prime bold upper X matrix has full rank. The general form of estimable functions is as shown in Table 3.

Table 3: General Form of Estimable Functions for a Multiple Regression Model When bold upper X prime bold upper X Matrix Is of Full Rank

Parameter Coefficient
beta 0 L1
beta 1 L2
beta 2 L3
beta 3 L4


For example, to test the hypothesis that beta 2 equals 0, let L1 = L2 = L4 = 0 and let L3 = 1. Then SSleft-parenthesis bold upper L bold-italic beta equals bold 0 right-parenthesis equals upper R left-parenthesis beta 2 vertical-bar beta 0 comma beta 1 comma beta 3 right-parenthesis. In this full-rank case, all parameters, as well as any linear combination of parameters, are estimable.

Suppose, however, that upper X Baseline 3 equals 2 x 1 plus 3 x 2. The general form of estimable functions is shown in Table 4.

Table 4: General Form of Estimable Functions for a Multiple Regression Model When bold upper X prime bold upper X Matrix Is Not of Full Rank

Parameter Coefficient
beta 0 L1
beta 1 L2
beta 2 L3
beta 3 2 times upper L Baseline italic 2 plus 3 times upper L Baseline italic 3


For this example, it is possible to test upper H 0 colon beta 0 equals 0. However, beta 1, beta 2, and beta 3 are not jointly estimable; that is,

StartLayout 1st Row 1st Column upper R left-parenthesis beta 1 vertical-bar beta 0 comma beta 2 comma beta 3 right-parenthesis 2nd Column equals 3rd Column 0 2nd Row 1st Column upper R left-parenthesis beta 2 vertical-bar beta 0 comma beta 1 comma beta 3 right-parenthesis 2nd Column equals 3rd Column 0 3rd Row 1st Column upper R left-parenthesis beta 3 vertical-bar beta 0 comma beta 1 comma beta 2 right-parenthesis 2nd Column equals 3rd Column 0 EndLayout
Last updated: December 09, 2022