Random-effects regression modelling is proposed for analysis of correlated grouped-time survival data. Two analysis approaches are considered. The first treats survival time as an ordinal outcome, which is either right-censored or not. The second approach treats survival time as a set of dichotomous indicators of whether the event occurred for time periods up to the period of the event or censor. For either approach both proportional hazards and proportional odds versions of the random-effects model are developed, while partial proportional hazards and odds generalizations are described for the latter approach. For estimation, a full-information maximum marginal likelihood solution is implemented using numerical quadrature to integrate over the distribution of multiple random effects. The quadrature solution allows some flexibility in the choice of distributions for the random effects; both normal and rectangular distributions are considered in this article. An analysis of a dataset where students are clustered within schools is used to illustrate features of random-effects analysis of clustered grouped-time survival data.
Get full access to this article
View all access options for this article.
References
1.
Thompson WA.On the treatment of grouped observations in life studies . Biometrics1977; 33: 463-470 .
2.
Prentice RL , Gloeckler LA.Regression analysis of grouped survival data with application to breast cancer data . Biometrics1978; 34: 57-67 .
3.
Allison PD.Discrete-time methods for the analysis of event histories. In Leinhardt S ed. Sociological methodology. San Francisco, CA: Jossey-Bass , 1982: 61-98.
4.
Aitkin M , Longford N.Statistical modelling issues in school effectiveness studies (with discussion) . Journal of the Royal Statistical Society, Series A1986; 149: 1-43 .
5.
Bryk AS , Raudenbush SW.Hierarchical linear models: applications and data analysis methods. London: Sage , 1992.
6.
Goldstein H.Multilevel statistical models, 2nd edn.London: Edward Arnold , 1995.
7.
Hedeker D , Gibbons RD , Flay BR.Random-effects regression models for clustered data: with an example from smoking prevention research . Journal of Consulting and Clinical Psychology1994; 62: 757-765 .
8.
Vaupel JW , Manton KG , Stallard E.The impact of heterogeneity in individual frailty on the dynamics of mortality . Demography1979; 16: 439-454 .
9.
Lancaster T.Econometric methods for the duration of unemployment . Econometrica1979; 47: 939-956 .
10.
Clayton D , Cuzick J.Multivariate generalizations of the proportional hazards model (with discussion) . Journal of the Royal Statistical Society, Series A1985; 148: 82-117 .
11.
Self SG , Prentice RL.Incorporating random effects into multivariate relative risk regression models. In Moolgavkar SH , Prentice RL eds. Modern statistical methods in chronic disease epidemiology. New York: John Wiley , 1986: 167-178.
12.
Guo G , Rodriquez G.Estimating a multivariate proportional hazards model for clustered data using the EM algorithm. With an application to child survival in Guatemala . Journal of the American Statistical Association1992; 87: 969-976 .
13.
Paik MC , Tsai W-Y , Ottman R.Multivariate survival analysis using piecewise gamma frailty . Biometrics1994; 50: 975-988 .
14.
Shih JH , Louis TA.Assessing gamma frailty models for clustered failure time data . Lifetime Data Analysis1995; 1: 205-220 .
15.
Pickles A , Crouchley R.A comparison of frailty models for multivariate survival data . Statistics in Medicine1995; 14: 1447-1461 .
16.
Hougaard P.Frailty models for survival data . Lifetime data analysis1995; 1: 255-273 .
17.
Liang K-Y , Zeger SL.Longitudinal data analysis using generalized linear models . Biometrika1986; 73: 13-22 .
18.
Lee EW , Wei LJ , Amato DA.Cox-type regression analysis for large numbers of small groups of correlated failure time observations. In Klein JP , Goel PK eds. Survival analysis: state of the art. Dordrecht: Kluwer , 1992: 237-247.
19.
Wei LJ , Lin DY , Weissfeld L.Regression analysis of multivariate incomplete failure time data by modeling marginal distributions . Journal of the American Statistical Association1989; 84: 1065-1073 .
20.
Han A , Hausman JA.Flexible parametric estimation of duration and competing risk models . Journal of Applied Econometrics1990; 5: 1-28 .
21.
Scheike TH , Jensen TK.A discrete survival model with random effects: an application to time to pregnancy . Biometrics1997; 53: 318-329 .
22.
Ten Have TR.A mixed effects model for multivariate ordinal response data including correlated discrete failure times with ordinal responses . Biometrics1996; 52: 473-491 .
23.
Ten Have TR , Uttal DH.Subject-specific and population-averaged continuation ratio logit models for multiple discrete time survival profiles . Applied Statistics1994; 43: 371-384 .
24.
Guo SW , Lin DY.Regression analysis of multivariate grouped survival data . Biometrics1994; 50: 632-639 .
25.
McCullagh P.Regression models for ordinal data (with discussion) . Journal of the Royal Statistical Society, Series B1980; 42: 109-142 .
26.
Teachman JD , Call VRA , Carver KP.Marital status and the duration of joblessness among white men . Journal of Marriage and the Family1994; 56: 415-428 .
27.
D'Agostino RB , Lee M-L , Belanger AJ , Cupples LA , Anderson K , Kannel WB.Relation of pooled logistic regression to time dependent Cox regression analysis: the Framingham Heart Study . Statistics in Medicine1990; 9: 1501-1515 .
28.
Singer JD , Willett JB.It's about time: using discrete-time survival analysis to study duration and the timing of events . Journal of Educational and Behavioral Statistics1993; 18: 155-195 .
29.
Hedeker D , Gibbons RD.A random-effects ordinal regression model for multilevel analysis . Biometrics1994; 50: 933-944 .
30.
Preisler HK.Analysis of a toxicological experiment using a generalized linear model with nested random effects . International Statistical Review1989; 57: 145-159 .
31.
Lee Y , Nelder JA.Hierarchical generalized linear models (with discussion) . Journal of the Royal Statistical Society, Series B1996; 58: 619-678 .
32.
Gibbons RD , Hedeker D.Random-effects probit and logistic regression models for three-level data . Biometrics1997; 53: 1527-1537 .
33.
Bock RD , Aitkin M.Marginal maximum likelihood estimation of item parameters: an application of the EM algorithm . Psychometrika1981; 46: 443-459 .
34.
Bennett S.Analysis of survival data by the proportional odds model . Statistics in Medicine1983; 2: 273-277 .
35.
Collett D.Modelling survival data in medical research. New York: Chapman & Hall , 1994.
36.
Doksum KA , Gasko M.On a correspondence between models in binary regression analysis and in survival analysis . International Statistical Review1990; 58: 243-252 .
37.
Cupples LA , D'Agostino RB , Anderson K , Kannel WB.Comparison of baseline and repeated measure covariate techniques in the Framingham Heart Study . Statistics in Medicine1985; 7: 205-218 .
38.
Engel J.On the analysis of grouped extreme-value data with GLIM . Applied Statistics1993; 42: 633-640 .
39.
Läärä E , Matthews JNS.The equivalence of two models for ordinal data . Biometrika1985; 72: 206-207 .
40.
Terza JV.Ordinal probit: a generalization . Communications in Statistical Theory and Methods1985; 14: 1-11 .
41.
Peterson B , Harrell FE.Partial proportional odds models for ordinal response variables . Applied Statistics1990; 39: 205-217 .
42.
Cox C.Location-scale cumulative odds models for ordinal data: a generalized nonlinear model approach . Statistics in Medicine1995; 14: 1191-1203 .
43.
Hedeker D , Mermelstein RJ.A multilevel thresholds of change model for analysis of stages of change data . Multivariate Behavioral Research1998; 33: 427-455 .
44.
Hedeker D , Gibbons RD.MIXOR: a computer program for mixed-effects ordinal probit and logistic regression analysis . Computer Methods and Programs in Biomedicine1996; 49: 157-176 .
Neuhaus JM , Kalbfleisch JD , Hauck WW.A comparison of cluster-specific and population-averaged approaches for analyzing correlated binary data . International Statistical Review1991; 59: 25-35 .
47.
Zeger SL , Liang K-Y , Albert PS.Models for longitudinal data: a generalized estimating equation approach . Biometrics1988; 44: 1049-1060 .
48.
Flay B , Brannon B , Johnson Cet al.The Television School and Family Smoking Prevention and Cessation Project: I. Theoretical Basis and Program Development . Preventive Medicine1988; 17: 585-607 .
49.
Agresti A.Categorical data analysis. New York: John Wiley , 1990.
50.
Siddiqui O , Hedeker D , Flay BR , Hu FB.Intraclass correlation estimates in a school-based smoking prevention study: outcome and mediating variables, by sex and ethnicity . Americal Journal of Epidemiology1996; 144: 425-433 .
Tanner MA.Tools for statistical inference: methods for the exploration of posterior distributions and likelihood functions, 3rd edn.New York: Springer , 1996.
55.
Bock RD , Shilling S.High-dimensional full-information item factor analysis. In Berkane M ed. Latent variable modeling and applications to causality. New York: Springer , 1997: 163-176.
56.
Allison PD.Survival analysis using the SAS system: a practical guide. Cary, NC: SAS Institute , 1995.
57.
Magnus JR.Linear structures. London: Charles Griffin , 1988.
58.
Bock RD.Measurement of human variation: a two-stage model. In Bock RD ed. Multilevel analysis of educational data. San Diego, CA: Academic Press , 1989: 319-342.
