The SURVEYPHREG Procedure

Replicate Weights Method

The replicate weights variance estimation method is a general-purpose variance estimation method that uses the replicate weights and replicate coefficients that you provide by using the REPWEIGHTS statement and the REPCOEFS= option, respectively.

If you provide your own replicate weights in a REPWEIGHTS statement but do not specify replicate coefficients in a REPCOEFS= option, then the default replicate coefficient depends on the VARMETHOD= option in the PROC SURVEYPHREG statement as shown in the following table:

Value of VARMETHOD= Default Replicate Coefficient
None specified left-parenthesis upper R minus 1 right-parenthesis slash upper R
BOOTSTRAP 1 slash upper R
BRR 1 slash upper R
JACKKNIFE left-parenthesis upper R minus 1 right-parenthesis slash upper R

Let ModifyingAbove bold-italic beta With caret be the estimated proportional hazards regression coefficients from the full sample, and let ModifyingAbove bold-italic beta With caret Subscript r and alpha Subscript r be the estimated regression coefficients and the replicate coefficient for the rth replicate, respectively. PROC SURVEYPHREG estimates the covariance matrix of ModifyingAbove bold-italic beta With caret by

ModifyingAbove bold upper V With caret left-parenthesis ModifyingAbove bold-italic beta With caret right-parenthesis equals sigma-summation Underscript r equals 1 Overscript upper R Endscripts alpha Subscript r Baseline left-parenthesis ModifyingAbove bold-italic beta With caret Subscript r Baseline minus ModifyingAbove bold-italic beta With caret right-parenthesis left-parenthesis ModifyingAbove bold-italic beta With caret Subscript r Baseline minus ModifyingAbove bold-italic beta With caret right-parenthesis prime

with R degrees of freedom, where R is the number of replicates.

If you specify the CENTER=REPLICATES method-option, then PROC SURVEYPHREG computes the covariance matrix of ModifyingAbove bold-italic beta With caret by

ModifyingAbove bold upper V With caret left-parenthesis ModifyingAbove bold-italic beta With caret right-parenthesis equals sigma-summation Underscript r equals 1 Overscript upper R Endscripts alpha Subscript r Baseline left-parenthesis ModifyingAbove bold-italic beta With caret Subscript r Baseline minus ModifyingAbove bold-italic beta With caret Subscript r Baseline overbar right-parenthesis left-parenthesis ModifyingAbove bold-italic beta With caret Subscript r Baseline minus ModifyingAbove bold-italic beta With caret Subscript r Baseline overbar right-parenthesis prime

where ModifyingAbove bold-italic beta With caret Subscript r Baseline overbar is the average of the replicate estimates as follows:

ModifyingAbove bold-italic beta With caret Subscript r Baseline overbar equals StartFraction 1 Over upper R EndFraction sigma-summation Underscript r equals 1 Overscript upper R Endscripts ModifyingAbove bold-italic beta Subscript r Baseline With caret

If one or more components of ModifyingAbove bold-italic beta With caret Subscript r cannot be calculated for some replicates, then the variance estimator uses only the replicates for which the proportional hazards regression coefficients can be estimated. Estimability and nonconvergence are two common reasons why ModifyingAbove bold-italic beta With caret Subscript r might not be available for a replicate sample even if ModifyingAbove bold-italic beta With caret is defined for the full sample. Let upper R Subscript a be the number of replicates where ModifyingAbove bold-italic beta With caret Subscript r are available, and let upper R minus upper R Subscript a be the number of replicates where ModifyingAbove bold-italic beta With caret Subscript r are not available. Without loss of generality, assume that ModifyingAbove bold-italic beta With caret Subscript r are available only for the first upper R Subscript a replicates; then the jackknife variance estimator is

ModifyingAbove bold upper V With caret left-parenthesis ModifyingAbove bold-italic beta With caret right-parenthesis equals sigma-summation Underscript r equals 1 Overscript upper R Subscript a Baseline Endscripts alpha Subscript r Baseline left-parenthesis ModifyingAbove bold-italic beta With caret Subscript r Baseline minus ModifyingAbove bold-italic beta With caret right-parenthesis left-parenthesis ModifyingAbove bold-italic beta With caret Subscript r Baseline minus ModifyingAbove bold-italic beta With caret right-parenthesis prime

with upper R Subscript a degrees of freedom. Alternatively, you can use the VADJUST=AVGREPSS option in the MODEL statement to use the average sum of squares for the invalid replicate samples. For more information, see Variance Adjustment Factors.

Last updated: December 09, 2022