The PHREG Procedure

Type 3 Tests and Joint Tests

For models that use less-than-full-rank parameterization (as specified by the PARAM=GLM option in the CLASS statement), a Type 3 test of an effect of interest (main effect or interaction) is a test of the Type III estimable functions that are defined for that effect. When the model contains no missing cells, performing the Type 3 test of a main effect corresponds to testing the hypothesis of equal marginal means. For more information about Type III estimable functions, see Chapter 53, The GLM Procedure, and Chapter 16, The Four Types of Estimable Functions. Also see Littell, Freund, and Spector (1991).

For models that use full-rank parameterization, all parameters are estimable when there are no missing cells, so it is unnecessary to define estimable functions. The standard test of an effect of interest in this case is the joint test that the values of the parameters associated with that effect are zero. For a model that uses effects parameterization (as specified by the PARAM=EFFECT option in the CLASS statement), performing the joint test for a main effect is equivalent to testing the equality of marginal means. For a model that uses reference parameterization (as specified by the PARAM=REF option in the CLASS statement), performing the joint test is equivalent to testing the equality of cell means at the reference level of the other model effects. For more information about the coding scheme and the associated interpretation of results, see Muller and Fetterman (2002, Chapter 14).

If there is no interaction term, the Type 3 test of an effect for a model that uses GLM parameterization is the same as the joint test of the effect for the model that uses full-rank parameterization. In this situation, the joint test is also called the Type 3 test. For a model that contains an interaction term and no missing cells, the Type 3 test of a component main effect under GLM parameterization is the same as the joint test of the component main effect under effect parameterization. Both test the equality of cell means. But this Type 3 test differs from the joint test under reference parameterization, which tests the equality of cell means at the reference level of the other component main effect. If some cells are missing, you can obtain meaningful tests only by testing a Type III estimable function, so in this case you should use GLM parameterization.

The results of a Type 3 test or a joint test do not depend on the order in which you specify the terms in the MODEL statement.

The following statistics can be used to test the null hypothesis upper H Subscript 0 upper L Baseline colon bold upper L bold-italic beta equals bold 0, where bold upper L is a matrix of known coefficients. Under large sample likelihood theory (Kalbfleisch and Prentice 1980, Section 3.4), each of the following statistics has each of the following statistics has an asymptotic chi-square distribution with p degrees of freedom, where p is the rank of bold upper L. Let bold-italic beta overTilde Subscript bold upper L be the maximum likelihood of bold-italic beta under the null hypothesis upper H Subscript 0 bold upper L; that is,

l left-parenthesis bold-italic beta overTilde Subscript bold upper L Baseline right-parenthesis equals max Underscript bold upper L bold-italic beta equals 0 Endscripts l left-parenthesis bold-italic beta right-parenthesis

Likelihood Ratio Statistic

chi Subscript normal upper L normal upper R Superscript 2 Baseline equals 2 left-bracket l left-parenthesis ModifyingAbove bold-italic beta With caret right-parenthesis minus l left-parenthesis bold-italic beta overTilde Subscript bold upper L Baseline right-parenthesis right-bracket

Score Statistic

chi Subscript upper S Superscript 2 Baseline equals left-bracket StartFraction partial-differential l left-parenthesis bold-italic beta overTilde Subscript bold upper L Baseline right-parenthesis Over partial-differential bold-italic beta EndFraction right-bracket prime left-bracket minus StartFraction partial-differential squared l left-parenthesis bold-italic beta overTilde Subscript bold upper L Baseline right-parenthesis Over partial-differential bold-italic beta squared EndFraction right-bracket Superscript negative 1 Baseline left-bracket StartFraction partial-differential l left-parenthesis bold-italic beta overTilde Subscript bold upper L Baseline right-parenthesis Over partial-differential bold-italic beta EndFraction right-bracket

Wald’s Statistic

chi Subscript upper W Superscript 2 Baseline equals left-parenthesis bold upper L ModifyingAbove bold-italic beta With caret right-parenthesis prime left-bracket bold upper L ModifyingAbove bold upper V With caret left-parenthesis ModifyingAbove bold-italic beta With caret right-parenthesis bold upper L prime right-bracket Superscript negative 1 Baseline left-parenthesis bold upper L ModifyingAbove bold-italic beta With caret right-parenthesis

where ModifyingAbove bold upper V With caret left-parenthesis ModifyingAbove bold-italic beta With caret right-parenthesis is the estimated covariance matrix, which can be the model-based covariance matrix left-bracket minus StartFraction partial-differential squared l left-parenthesis ModifyingAbove bold-italic beta With caret right-parenthesis Over partial-differential bold-italic beta squared EndFraction right-bracket Superscript negative 1 or the sandwich covariance matrix upper V Subscript upper S Baseline left-parenthesis ModifyingAbove bold-italic beta With caret right-parenthesis. For more information, see the section Robust Sandwich Variance Estimate.

Last updated: March 08, 2022