The HPLMIXED Procedure

Statistical Properties

If bold upper G and bold upper R are known, ModifyingAbove bold-italic beta With caret is the best linear unbiased estimator (BLUE) of bold-italic beta, and ModifyingAbove bold-italic gamma With caret is the best linear unbiased predictor (BLUP) of bold-italic gamma (Searle 1971; Harville 1988, 1990; Robinson 1991; McLean, Sanders, and Stroup 1991). Here, "best" means minimum mean squared error. The covariance matrix of left-parenthesis ModifyingAbove bold-italic beta With caret minus bold-italic beta comma ModifyingAbove bold-italic gamma With caret minus bold-italic gamma right-parenthesis is

bold upper C equals Start 2 By 2 Matrix 1st Row 1st Column bold upper X prime bold upper R Superscript negative 1 Baseline bold upper X 2nd Column bold upper X prime bold upper R Superscript negative 1 Baseline bold upper Z 2nd Row 1st Column bold upper Z prime bold upper R Superscript negative 1 Baseline bold upper X 2nd Column bold upper Z prime bold upper R Superscript negative 1 Baseline bold upper Z plus bold upper G Superscript negative 1 EndMatrix Superscript minus

where Superscript minus denotes a generalized inverse ( Searle 1971).

However, bold upper G and bold upper R are usually unknown and are estimated by using one of the aforementioned methods. These estimates, ModifyingAbove bold upper G With caret and ModifyingAbove bold upper R With caret, are therefore simply substituted into the preceding expression to obtain

ModifyingAbove bold upper C With caret equals Start 2 By 2 Matrix 1st Row 1st Column bold upper X prime ModifyingAbove bold upper R With caret Superscript negative 1 Baseline bold upper X 2nd Column bold upper X prime ModifyingAbove bold upper R With caret Superscript negative 1 Baseline bold upper Z 2nd Row 1st Column bold upper Z prime ModifyingAbove bold upper R With caret Superscript negative 1 Baseline bold upper X 2nd Column bold upper Z prime ModifyingAbove bold upper R With caret Superscript negative 1 Baseline bold upper Z plus ModifyingAbove bold upper G With caret Superscript negative 1 EndMatrix Superscript minus

as the approximate variance-covariance matrix of left-parenthesis ModifyingAbove bold-italic beta With caret minus bold-italic beta comma ModifyingAbove bold-italic gamma With caret minus bold-italic gamma). In this case, the BLUE and BLUP acronyms no longer apply, but the word empirical is often added to indicate such an approximation. The appropriate acronyms thus become EBLUE and EBLUP.

McLean and Sanders (1988) show that ModifyingAbove bold upper C With caret can also be written as

ModifyingAbove bold upper C With caret equals Start 2 By 2 Matrix 1st Row 1st Column ModifyingAbove bold upper C With caret Subscript 11 Baseline 2nd Column ModifyingAbove bold upper C With caret Subscript 21 Superscript prime Baseline 2nd Row 1st Column ModifyingAbove bold upper C With caret Subscript 21 Baseline 2nd Column ModifyingAbove bold upper C With caret Subscript 22 EndMatrix

where

StartLayout 1st Row 1st Column ModifyingAbove bold upper C With caret Subscript 11 2nd Column equals left-parenthesis bold upper X prime ModifyingAbove bold upper V With caret Superscript negative 1 Baseline bold upper X right-parenthesis Superscript minus Baseline 2nd Row 1st Column ModifyingAbove bold upper C With caret Subscript 21 2nd Column equals minus ModifyingAbove bold upper G With caret bold upper Z prime ModifyingAbove bold upper V With caret Superscript negative 1 Baseline bold upper X ModifyingAbove bold upper C With caret Subscript 11 Baseline 3rd Row 1st Column ModifyingAbove bold upper C With caret Subscript 22 2nd Column equals left-parenthesis bold upper Z prime ModifyingAbove bold upper R With caret Superscript negative 1 Baseline bold upper Z plus ModifyingAbove bold upper G With caret Superscript negative 1 Baseline right-parenthesis Superscript negative 1 Baseline minus ModifyingAbove bold upper C With caret Subscript 21 Baseline bold upper X prime ModifyingAbove bold upper V With caret Superscript negative 1 Baseline bold upper Z ModifyingAbove bold upper G With caret EndLayout

Note that ModifyingAbove bold upper C With caret Subscript 11 is the familiar estimated generalized least squares formula for the variance-covariance matrix of ModifyingAbove bold-italic beta With caret.

Last updated: December 09, 2022