Introduction to Statistical Modeling with SAS/STAT Software


Estimable Functions

A function $\bL \bbeta $ is said to be estimable if there exists a linear combination of the expected value of $\bY $, such as $\bK \mr{E}[\bY ]$, that equals $\bL \bbeta $. Since $\mr{E}[\bY ] = \bX \bbeta $, the definition of estimability implies that $\bL \bbeta $ is estimable if there is a matrix $\bK $ such that $\bL = \bK \bX $. Another way of looking at this result is that the rows of $\bX $ form a generating set from which all estimable functions can be constructed.

The concept of estimability of functions is important in the theory and application of linear models because hypotheses of interest are often expressed as linear combinations of the parameter estimates (for example, hypotheses of equality between parameters, $\beta _1=\beta _2 \Leftrightarrow \beta _1-\beta _2 = 0$). Since estimability is not related to the particular value of the parameter estimate, but to the row space of $\bX $, you can test only hypotheses that consist of estimable functions. Further, because estimability is not related to the value of $\bbeta $ (Searle, 1971, p. 181), the choice of the generalized inverse in a situation with rank-deficient $\bX ’\bX $ matrix is immaterial, since

\[  \bL \widehat{\bbeta } = \bK \bX \widehat{\bbeta } = \bK \bX \left(\bX ’\bX \right)^{-}\bX ’\bY  \]

where $\bX \left(\bX ’\bX \right)^{-}\bX $ is invariant to the choice of generalized inverse.

$\bL \bbeta $ is estimable if and only if $\bL (\bX ’\bX )^{-}(\bX ’\bX ) = \bL $ (see, for example, Searle 1971, p. 185). If $\bX $ is of full rank, then the Hermite matrix $(\bX ’\bX )^{-}(\bX ’\bX )$ is the identity, which implies that all linear functions are estimable in the full-rank case.

See Chapter 15: The Four Types of Estimable Functions, for many details about the various forms of estimable functions in SAS/STAT.