Prints output for a proportional hazards model fitted by RPHFIT (R.W. Payne).
Option
PRINT = string tokens |
Controls printed output (model, deviance, summary, estimates, correlations, fittedvalues, accumulated, loglikelihood); default mode, summ, esti |
|---|
No parameters
Description
This procedure allows you to display further output for a proportional hazards models that has been fitted by procedure RPHFIT. The output is controlled by PRINT option with similar settings to those of the FIT directive, except that there is an extra setting loglikelihood to print -2 times the log-likelihood and the number of degrees of freedom in the model. The deviance produced for the terms in the regression model can be assessed using chi-square distributions as usual, but the residual deviance is not usable as the maximal model assumed by the generalized linear models method is inappropriate. So, the residual line is suppressed in the summary and accumulated analysis of deviance.
Option: PRINT.
Method
Further details of the method used by RPHFIT can be found in Aitkin et al. (1989).
Action with RESTRICT
None of the vectors must be restricted (and any restrictions will have been cancelled by RPHFIT).
Reference
Aitkin, M., Anderson, A., Francis, B. & Hinde, J. (1989). Statistical Modelling in GLIM. Oxford University Press.
See also
Procedures: KAPLANMEIER, RLIFETABLE, RPHFIT, RPHCHANGE, RPHKEEP, RPROPORTIONAL, RSTEST, RSURVIVAL.
Commands for: Survival analysis.
Example
CAPTION 'RPHDISPLAY example',\
'Data from Gehan (1965, Biometrika, 52, 203-223).';\
STYLE=meta,plain
VARIATE [VALUES=1,1,2,2,3,4,4,5,5,8,8,8,8,11,11,12,12,15,17,22,23,\
6,6,6,6,7,9,10,10,11,13,16,17,19,20,22,23,25,32,32,34,35] Time
& [VALUES=24(0),1,0,1,0,1,1,0,0,1,1,1,0,0,1,1,1,1,1] Censor
FACTOR [LABELS=!t(control,'6-mercaptopurine'); VALUES=21(1,2)] Treat
FACTOR [LEVELS=42; VALUES=1...42] Subject
RPHFIT [PRINT=*; TIMES=Time; SUBJECTS=Subject; CENSORED=Censor] Treat
RPHDISPLAY