This example uses a data set from the National Health and Nutrition Examination Survey I (NHANES I) Epidemiologic Followup Study (NHEFS). The NHEFS is a national longitudinal survey that is conducted by the National Center for Health Statistics, the National Institute on Aging, and some other agencies of the Public Health Service in the United States. Some important objectives of this survey are to determine the relationships between clinical, nutritional, and behavioral factors; to determine mortality and hospital utilizations; and to monitor changes in risk factors for the initial cohort that represents the NHANES I population. A cohort of size 14,407, which includes all persons 25 to 74 years old who completed a medical examination at NHANES I in 1971–1975, was selected for the NHEFS. Personal interviews were conducted for every selected unit during the first wave of data collection from the year 1982 to 1984. Follow-up studies were conducted in 1986, 1987, and 1992. In the year 1986, only nondeceased persons 55 to 74 years old (as reported in the base year survey) were interviewed. The 1987 and 1992 NHEFS contain the entire nondeceased NHEFS cohort. Vital and tracing status data, interview data, health care facility stay data, and mortality data for all four waves are available for public use. See http://www.cdc.gov/nchs/nhanes/nhefs/nhefs.htm for more information about the survey and the data sets.
For illustration purposes, 1,018 observations from the 1987 NHEFS public use interview data are used to create the data set
cancer
. The observations are obtained from 10 strata that contain 596 PSUs. The sum of observation weights for these selected units
is over 19 million. Observation weights range from 359 to 129,359 with a mean of 18,747.69 and a median of 11,414. Several
observation weights have large values; therefore it is reasonable to rescale the observation weights to facilitate the optimization
routine. Different scaling techniques are proposed in the literature. For example, Binder (1992) uses scaled weights such that the sum of weights over the sampled units is one. Without loss of generality, the analysis
weights in this example are obtained by dividing each observation weight by a large number (130,000). Because of this rescaling,
you must be careful interpreting some results from PROC SURVEYPHREG.
The following variables are used in this example:
ObsNo
, unit identification
Strata
, stratum identification
PSU
, identification for primary sampling units
ObservationWt
, sampling weight associated with each unit
AnalysisWt
, obtained from the sampling weights by dividing each ObservationWt
by 130,000
Smoke
, smoking status (–1 = not applicable, 1 = never smoked, 2 = current or former smoker in 1982–1984 follow-up, and 3 = current
or former smoker in 1987 follow-up)
Age
, the event-time variable, defined as follows:
age of the subject when the first cancer was reported for subjects with reported cancer
age of the subject at death for deceased subjects without reported cancer
age of the subject as reported in 1987 follow-up (this value is used for nondeceased subjects who never reported cancer)
age of the subject for the entry year 1971–1975 survey if the subject has cancer (or is deceased) but the date of incident is not reported
Cancer
, cancer indicator (1 = cancer reported, 0 = cancer not reported)
BodyWeight
, body weight of the subject as reported in the 1987 follow-up, or an imputed body weight based on the subject’s age in the
entry year 1971–1975 survey
The following SAS statements create the data set cancer
. Note that BodyWeight
for a few observations (8%) is imputed based on Age
by using a deterministic regression imputation model (Särndal and Lundström (2005, chapter 12)). The imputed values are treated as observed values in this example. In other words, this example treats the
data set Cancer
as the observed data set.
data cancer; input ObsNo Strata PSU AnalysisWt ObservationWt Smoke Age Cancer BodyWeight; datalines; 1 3 002 0.02927 3805 2 53 1 175 2 3 002 0.04698 6107 2 77 0 175 3 3 039 0.02283 2968 2 50 0 160 4 3 084 0.23414 30438 2 52 0 145 5 3 007 0.03908 5081 1 80 0 127 6 3 009 0.02993 3891 1 62 0 180 7 3 009 0.02754 3580 2 50 0 157 8 3 022 0.02283 2968 2 56 0 142 9 3 050 0.18268 23748 2 60 0 140 ... more lines ... 1016 4 002 0.02068 2689 2 40 0 120 1017 4 092 0.35298 45888 2 52 0 166 1018 4 035 0.03344 4347 -1 58 0 156 ;
Suppose you want to study the occurrence of cancer for the base year survey population and its relation to smoking status
and body weight. The following statements request a proportional hazards regression of Age
on BodyWeight
and Smoke
with Cancer
as the censor indicator. The STRATA, CLUSTER, and WEIGHT statements identify the variance strata, PSUs, and analysis weights
respectively. The CLASS statement specifies that Smoke
is a categorical variable, and the MODEL statement provides information about the analysis model. The TIES= option in the
MODEL statement requests the Efron likelihood to handle tied events. If you do not specify the TIES= option in the MODEL statement,
then the procedure uses the Breslow likelihood. The PHISTORY option in the NLOPTIONS statement is used to display the iteration
history table. The ESTIMATE statement computes a contrast between subjects who are reported as current (or former) smokers
and the others. The EXP option in the ESTIMATE statement requests that the linear contrast be estimated in the exponential
scale, which is the hazard ratio. The TEST statement requests the Type 3 test for each effect that is specified in the MODEL
statement.
proc surveyphreg data = cancer; strata strata; cluster psu; weight analysiswt; class smoke; model age*cancer(0) = bodyweight smoke / ties = efron; nloptions phistory; estimate smoke 0.5 0.5 -0.5 -0.5 / exp; test ; run;
Some summary statistics are shown in Output 93.2.1. The “Model Information” table contains information about the model such as the names for the dependent and censoring variables, and the likelihood.
The “Number of Observations” table displays the number of observations and the sum of weights. A total of 1,018 observations are read from the Cancer
data set, but one observation is not used in the analysis because it has a zero sampling weight. The sum of weights is 146.81,
which gives an estimated population size of 19,085,105 (). Note that the estimated population size would be 19,085,151 if you use the sampling weights (ObservationWt
) instead of the analysis weights (AnalysisWt
). The difference is due to the rounding errors in AnalysisWt
. For simplicity, analysis weights are rounded at the fifth decimal place. The “Design Summary” table shows that there are 596 PSUs and 10 strata. From the censored summary tables, 11.7% subjects in the sample have reported
cancer and an estimated 11.6% subjects in the study population have cancer. The “Variance Estimation” table shows that the Taylor series linearization variance estimation method is used and the observation units with missing
values are excluded from the analysis. Note that the only missing unit in this data set has a zero sampling weight and hence
it is not included in the analysis.
Output 93.2.1: Model Information, Data Summary, Design Summary, and Information about Variance Estimation
Model Information | |
---|---|
Data Set | WORK.CANCER |
Dependent Variable | Age |
Censoring Variable | Cancer |
Censoring Value(s) | 0 |
Weight Variable | AnalysisWt |
Stratum Variable | Strata |
Cluster Variable | PSU |
Ties Handling | EFRON |
Number of Observations Read | 1018 |
---|---|
Number of Observations Used | 1017 |
Sum of Weights Read | 146.8085 |
Sum of Weights Used | 146.8085 |
Design Summary | |
---|---|
Number of Strata | 10 |
Number of Clusters | 596 |
Summary of the Number of Event and Censored Values |
|||
---|---|---|---|
Total | Event | Censored | Percent Censored |
1017 | 119 | 898 | 88.30 |
Summary of the Weighted Number of Event and Censored Values |
|||
---|---|---|---|
Total | Event | Censored | Percent Censored |
146.8085 | 17.01185 | 129.7966 | 88.41 |
Variance Estimation | |
---|---|
Method | Taylor Series |
Missing Values | Excluded |
The “Iteration History” table in Output 93.2.2 shows that the procedure converged after four iterations. The “Objective Function” column contains the value of the likelihood after every iteration. The “Objective Function Change” column measures the change in the objective function between iterations; however, this is not the monitored convergence criterion. The SURVEYPHREG procedure monitors several features simultaneously to determine whether to stop an optimization.
Output 93.2.2: Iteration History
Maximum Likelihood Iteration History | ||||||||||
---|---|---|---|---|---|---|---|---|---|---|
Iteration | Restarts | Function Calls |
Active Constraints |
Objective Function |
Objective Function Change |
Max Abs Gradient Element |
Ridge | Ratio Between Actual and Predicted Change |
||
1 | 0 | 4 | 0 | -63.34004 | 1.6501 | 21.9620 | 0 | 0.916 | ||
2 | 0 | 6 | 0 | -63.29819 | 0.0418 | 0.2005 | 0 | 1.052 | ||
3 | 0 | 8 | 0 | -63.29776 | 0.000430 | 0.00293 | 0 | 1.012 | ||
4 | 0 | 10 | 0 | -63.29776 | 1.528E-7 | 1.102E-6 | 0 | 1.000 |
Estimates for proportional hazards regression coefficients and their standard errors are shown in Output 93.2.3. The categorical variable Smoke
has four levels, and GLM parameterization is used by PROC SURVEYPHREG. You can use the PARAM= option in the CLASS statement
to specify other types of parameterizations. The estimated regression coefficient for BodyWeight
is 0.012 with a standard error of 0.003. The degrees of freedom for the t test are equal to the number of PSUs (596) minus the number of strata (10). The “Estimates” table displays the estimated contrast and the corresponding hypothesis test. The estimated value for the contrast is –0.75.
The estimated hazard for the nonsmokers is 0.47 times the estimated hazard for the current or former smokers. In this example
data set, the contrast of interest is not significant at 0.05 levels. The “Type III Tests of Model Effects” table displays the Type 3 analysis. The effect variable Smoke
has four levels. The F Value for Smoke
is 1.49 with three numerator degrees of freedom and 586 denominator degrees of freedom.
Output 93.2.3: Parameter Estimates
Analysis of Maximum Likelihood Estimates | ||||||
---|---|---|---|---|---|---|
Parameter | DF | Estimate | Standard Error | t Value | Pr > |t| | Hazard Ratio |
BodyWeight | 586 | 0.011920 | 0.003155 | 3.78 | 0.0002 | 1.012 |
Smoke -1 | 586 | -1.174048 | 0.739450 | -1.59 | 0.1129 | 0.309 |
Smoke 1 | 586 | -1.006515 | 0.578810 | -1.74 | 0.0826 | 0.365 |
Smoke 2 | 586 | -0.674183 | 0.558412 | -1.21 | 0.2278 | 0.510 |
Smoke 3 | 586 | 0 | . | . | . | 1.000 |
Type III Tests of Model Effects | ||||
---|---|---|---|---|
Effect | Num DF | Den DF | F Value | Pr > F |
BodyWeight | 1 | 586 | 14.27 | 0.0002 |
Smoke | 3 | 586 | 1.49 | 0.2160 |
Estimate | ||||||
---|---|---|---|---|---|---|
Label | Estimate | Standard Error | DF | t Value | Pr > |t| | Exponentiated |
Row 1 | -0.7532 | 0.3870 | 586 | -1.95 | 0.0521 | 0.4709 |