The GLM Procedure

ESTIMATE Statement

  • ESTIMATE ’label’ effect values <…effect values> </ options>;

The ESTIMATE statement enables you to estimate linear functions of the parameters by multiplying the vector bold upper L by the parameter estimate vector bold b, resulting in bold upper L bold b. All of the elements of the bold upper L vector might be given, or, if only certain portions of the bold upper L vector are given, the remaining elements are constructed by PROC GLM from the context (in a manner similar to rule 4 discussed in the section Construction of Least Squares Means).

The linear function is checked for estimability. The estimate bold upper L bold b, where bold b equals left-parenthesis bold upper X prime bold upper X right-parenthesis Superscript minus Baseline bold upper X prime bold y, is displayed along with its associated standard error, StartRoot bold upper L left-parenthesis bold upper X prime bold upper X right-parenthesis Superscript minus Baseline bold upper L prime s squared EndRoot, and t test. If you specify the CLPARM option in the MODEL statement, confidence limits for the true value are also displayed.

There is no limit to the number of ESTIMATE statements that you can specify, but they must appear after the MODEL statement. In the ESTIMATE statement,

label

identifies the estimate on the output. A label is required for every contrast specified. Labels must be enclosed in quotes.

effect

identifies an effect that appears in the MODEL statement, or the INTERCEPT effect. The INTERCEPT effect can be used as an effect when an intercept is fitted in the model. You do not need to include all effects that are in the MODEL statement.

values

are constants that are the elements of the bold upper L vector associated with the preceding effect. For example, 

estimate 'A1 VS A2' A  1 -1;

forms an estimate that is the difference between the parameters estimated for the first and second levels of the CLASS variable A.

You can specify the following options in the ESTIMATE statement after a slash (/):

DIVISOR=number

specifies a value by which to divide all coefficients so that fractional coefficients can be entered as integer numerators. For example, you can use 

estimate '1/3(A1+A2) - 2/3A3' a 1 1 -2 / divisor=3;

instead of 

estimate '1/3(A1+A2) - 2/3A3' a 0.33333 0.33333 -0.66667;
E

displays the entire bold upper L vector. This option is useful in confirming the ordering of parameters for specifying bold upper L.

SINGULAR=number

tunes the estimability checking. If ABSleft-parenthesis bold upper L minus bold upper L bold upper H right-parenthesis greater-than upper C timesnumber, then the bold upper L vector is declared nonestimable. bold upper H is the left-parenthesis bold upper X prime bold upper X right-parenthesis Superscript minus Baseline bold upper X prime bold upper X matrix, and C is ABSleft-parenthesis bold upper L right-parenthesis except for rows where bold upper L is zero, and then it is 1. The default value for the SINGULAR= option is 10 Superscript negative 4. Values for the SINGULAR= option must be between 0 and 1.

See also the section Specification of ESTIMATE Expressions.

Last updated: December 09, 2022