The GLMPOWER Procedure

CONTRAST Statement

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

The CONTRAST statement enables you to define custom Type III hypothesis tests by specifying an bold upper L vector or matrix for testing either the hypothesis bold upper L bold-italic beta equals 0 (for univariate models) or the hypothesis bold upper L bold upper B bold upper M equals 0 (for multivariate models). The bold upper L matrix consists of one or more between-subject contrasts.

To use this feature, you must be familiar with the details of the model parameterization that PROC GLM uses. For more information, see the section Parameterization of PROC GLM Models in Chapter 53, The GLM Procedure. All the elements of the bold upper L matrix can be given, or if only certain portions of the bold upper L matrix are given, PROC GLMPOWER constructs the remaining elements from the context (in a manner similar to that in rule 4 in the section Construction of Least Squares Means in Chapter 53, The GLM Procedure).

There is no limit to the number of CONTRAST statements that you can specify, but they must appear after the MODEL statement. Each power analysis includes tests for all CONTRAST statements.

You can specify the following arguments:

label

identifies the contrast in the output. A label is required for every contrast that is specified. Labels must be enclosed in single or double quotation marks.

effect

identifies an effect that appears in the MODEL statement, or the INTERCEPT effect. You do not need to include all effects that appear in the MODEL statement.

values

are constants that are elements of the bold upper L matrix associated with the effect.

You can specify the following option in the CONTRAST statement after a slash (/):

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 times number for any row in the contrast, then bold upper L 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.

As stated previously, the CONTRAST statement enables you to define custom hypothesis tests. If the hypothesis is testable in a univariate model, then the hypothesis sum of squares, SS(upper H 0 colon bold upper L bold-italic beta equals 0), is computed as

left-parenthesis bold upper L bold b right-parenthesis prime left-parenthesis bold upper L left-parenthesis bold upper X prime bold upper X right-parenthesis Superscript minus Baseline bold upper L prime right-parenthesis Superscript negative 1 Baseline left-parenthesis bold upper L bold b right-parenthesis

where bold b equals left-parenthesis bold upper X prime bold upper X right-parenthesis Superscript minus Baseline bold upper X prime bold y.

For testable hypotheses in a multivariate model, the usual multivariate tests are defined by using

bold upper H equals bold upper M prime left-parenthesis bold upper L bold upper B right-parenthesis prime left-parenthesis bold upper L left-parenthesis bold upper X prime bold upper X right-parenthesis Superscript minus Baseline bold upper L prime right-parenthesis Superscript negative 1 Baseline left-parenthesis bold upper L bold upper B right-parenthesis bold upper M

where bold upper B equals left-parenthesis bold upper X prime bold upper X right-parenthesis Superscript minus Baseline bold upper X prime bold upper Y and bold upper Y is the matrix of multivariate responses or dependent variables.

The degrees of freedom associated with the hypothesis are equal to the row rank of bold upper L. The sum of squares computed in this situation is equivalent to the sum of squares computed using an bold upper L matrix with any row deleted that is a linear combination of previous rows.

Multiple-degrees-of-freedom hypotheses can be specified by separating the rows of the bold upper L matrix with commas.

Last updated: December 09, 2022