Score test of homogeneity for survival data

If follow-up is made for subjects which are grouped into units, such as familial or spatial units then it may be interesting to test whether the groups are homogeneous (or independent for given explanatory variables). The effect of the groups is modelled as random and we consider a frailty proportional hazards model which allows to adjust for explanatory variables. We derive the score test of homogeneity from the marginal partial likelihood and it turns out to be the sum of a pairwise correlation term of martingale residuals and an overdispersion term. In the particular case where the sizes of the groups are equal to one, this statistic can be used for testing overdispersion. The asymptotic variance of this statistic is derived using counting process arguments. An extension to the case of several strata is given. The resulting test is computationally simple; its use is illustrated using both simulated and real data. In addition a decomposition of the score statistic is proposed as a sum of a pairwise correlation term and an overdispersion term. The pairwise correlation term can be used for constructing a statistic more robust to departure from the proportional hazard model, and the overdispesion term for constructing a test of fit of the proportional hazard model.

[1]  J. R. Batchelor,et al.  HL-A matching in treatment of burned patients with skin allografts. , 1970, Lancet.

[2]  N. Breslow.  Covariance analysis of censored survival data. , 1974, Biometrics.

[3]  D. Clayton.  A model for association in bivariate life tables and its application in epidemiological studies of familial tendency in chronic disease incidence , 1978 .

[4]  R. Sutton,et al.  β-BLOCKERS, AUTOIMMUNITY, AND RASHES , 1980, The Lancet.

[5]  David R. Cox,et al.  Some remarks on overdispersion , 1983 .

[6]  D. Clayton,et al.  Multivariate generalizations of the proportional hazards model , 1985 .

[7]  A. Yashin,et al.  Heterogeneity's ruses: some surprising effects of selection on population dynamics. , 1985, The American statistician.

[8]  P. Hougaard.  A class of multivanate failure time distributions , 1986 .

[9]  M. Gail,et al.  Radiation-induced thyroid neoplasms: evidence for familial susceptibility factors. , 1988, The Journal of clinical endocrinology and metabolism.

[10]  O. Aalen,et al.  Heterogeneity in survival analysis. , 1988, Statistics in medicine.

[11]  D. Harrington,et al.  Counting Processes and Survival Analysis , 1991 .

[12]  C. Dean.  Testing for Overdispersion in Poisson and Binomial Regression Models , 1992 .

[13]  Richard D. Gill,et al.  A counting process approach to maximum likelihood estimation in frailty models , 1992 .

[14]  Niels Keiding,et al.  Statistical Models Based on Counting Processes , 1993 .

[15]  D. Commenges,et al.  Test of homogeneity of binary data with explanatory variables. , 1994, Biometrics.

[16]  D. Commenges,et al.  Tests of Homogeneity for Generalized Linear Models , 1995 .

[17]  R. Gray.  Tests for Variation over Groups in Survival Data , 1995 .