For the latest version of JMP Help, visit JMP.com/help.


Reliability and Survival Methods > Fit Proportional Hazards > Example Using Multiple Effects and Multiple Levels
Publication date: 11/29/2021

Example Using Multiple Effects and Multiple Levels

This example uses a proportional hazards model for the sample data, VA Lung Cancer.jmp. The data were collected from a randomized clinical trial, where males with inoperable lung cancer were placed on either a standard or a novel (test) chemotherapy (Treatment). The primary interest of this trial was to assess if the treatment type has an effect on survival time; special interest is given to the type of tumor (Cell Type). See Prentice (1973) and Kalbfleisch and Prentice (2002) for more information about this data set.

For the proportional hazards model, covariates include the following:

Whether the patient had undergone previous therapy (Prior)

Age of the patient (Age)

Time from lung cancer diagnosis to beginning the study (Diag Time)

A general medical status measure (KPS)

Age, Diag Time, and KPS are continuous measures and Cell Type, Treatment, and Prior are categorical (nominal) variables. The four nominal levels of Cell Type include Adeno, Large, Small, and Squamous.

This example illustrates the results for a model with more than one effect and a nominal effect with more than two levels. Risk ratios are also demonstrated, with example calculations for risk ratios for a continuous effect and risk ratios for an effect that has more than two levels.

1. Select Help > Sample Data Library and open VA Lung Cancer.jmp.

2. Select Analyze > Reliability and Survival > Fit Proportional Hazards.

3. Select Time as Time to Event.

4. Select censor as Censor.

5. Select Cell Type, Treatment, Prior, Age, Diag Time, and KPS and click Add.

6. Click Run.

7. (Optional) Click the disclosure icon on the Baseline Survival at mean title bar to close the plot, and click the disclosure icon on Effect Summary to close the report.

Figure 15.6 Report Window for Proportional Hazards Model with Multiple Effects and Levels 

Report Window for Proportional Hazards Model with Multiple Effects and Levels

Note the following about the results:

In the Whole Model report, the low Prob>Chisq value (<.0001) indicates that there is a difference in survival time when at least one of the effects is included in the model.

In the Effect Likelihood Ratio Tests report, the Prob>ChiSq values indicate that KPS and at least one of the levels of Cell Type are significant; however, the Treatment, Prior, Age, and Diag Time effects are not significant.

8. Click the Proportional Hazards Fit red triangle and select Risk Ratios.

Figure 15.7 Risk Ratios Report 

Risk Ratios Report

Risk Ratios for Multiple Effects and Multiple Levels

Figure 15.7 shows the Risk Ratios for the continuous effects (Age, Diag Time, KPS) and the nominal effects (Cell Type, Treatment, Prior). For illustration, focus on the continuous effect, Age, and the nominal effect with four levels (Cell Type) for the VA Lung Cancer.jmp sample data.

For the continuous effect, Age, in the VA Lung Cancer.jmp sample data, the risk ratios are calculated as follows:

Unit Risk Ratios

exp(β) = exp(-0.0085494) = 0.991487

Range Risk Ratios

exp[β(xmax - xmin)] = exp(-0.0085494 * 47) = 0.669099

For the nominal effect, Cell Type, all pairs of levels are calculated and are shown in the Risk Ratios for Cell Type table. Note that for a categorical variable with k levels, only k -1 design variables, or levels, are used. In the Parameter Estimates table, parameter estimates are shown for only three of the four levels for Cell Type (Adeno, Large, and Small). The Squamous level is not shown, but it is calculated as the negative sum of the other estimates. Here are two example Risk Ratios for Cell Type calculations:

Large/Adeno = exp(βLarge)/exp(βAdeno) = exp(-0.2114757)/exp(0.57719588) = 0.4544481

Squamous/Adeno = exp[-(βAdeno + βLarge + βSmall)]/exp(βAdeno)

= exp[-(0.57719588 + (-0.2114757) + 0.24538322)]/exp(0.57719588) = 0.3047391

Reciprocal shows the value for 1/Risk Ratio.

Want more information? Have questions? Get answers in the JMP User Community (community.jmp.com).