This Article 
   
 Share 
   
 Bibliographic References 
   
 Add to: 
 
Digg
Furl
Spurl
Blink
Simpy
Google
Del.icio.us
Y!MyWeb
 
 Search 
   
A Unified Approach for Modeling Longitudinal and Failure Time Data, With Application in Medical Monitoring
February 1996 (vol. 18 no. 2)
pp. 109-123

Abstract—This paper considers biomedical problems in which a sample of subjects, for example clinical patients, is monitored through time for purposes of individual prediction. Emphasis is on situations in which the monitoring generates data both in the form of a time series and in the form of events (development of a disease, death, etc.) observed on each subject over specified intervals of time. A Bayesian approach to the combined modeling of both types of data for purposes of prediction is presented. The proposed method merges ideas of Bayesian hierarchical modeling, nonparametric smoothing of time series data, survival analysis, and forecasting into a unified framework. Emphasis is on flexible modeling of the time series data based on stochastic process theory. The use of Markov Chain Monte Carlo simulation to calculate the predictions of interest is discussed. Conditional independence graphs are used throughout for a clear presentation of the models. An application in the monitoring of transplant patients is presented.

[1] E. Arjas and D. Gasbarra, "Nonparametric Bayesian inference from right censored survival data, using the Gibbs sampler," Statistica Sinica (to appear).
[2] C. Berzuini, R. Bellazzi, S. Quaglini, and D.J. Spiegelhalter, "Bayesian networks for patient monitoring," Artificial Intelligence in Medicine, vol. 4, pp. 243-260, 1992.
[3] C. Berzuini, D. Clayton, and L. Bernardinelli, "Bayesian inference on the Lexis diagram," Bull. Int'l Statistical Institute, vol. LV, pp. 149-165, 1993.
[4] C. Berzuini and D. Clayton, "Bayesian analysis of survival on multiple time scales," Statistics in Medicine, vol. 13, pp. 823-838, 1994.
[5] C. Berzuini, N. Best, W. Gilks, and C. Larizza, "Dynamic graphical models and Markov chain Monte Carlo methods," technical report, Dipartimento di Informatica e Sistemistica, Universitádi Pavia, Italy, 1994.
[6] N.E. Breslow and N. Day, Statistical Methods in Cancer Research, Vol. II—The Design and Analysis of Cohort Studies.Lyon, France: I.A.R.C. Scientific Publications, 1987.
[7] N.E. Breslow and D. Clayton, "Approximate inference in generalized linear mixed models," J. Am. Statistical Soc., vol. 88, pp. 9-25, 1993.
[8] E.M. Chi and G.C. Reinsel, "Models for longitudinal data with random effects and AR(1) errors," J. Am. Statistical Assoc., vol. 84, pp. 452-459, 1989.
[9] D. Clayton, "The analysis of event history data: A review of progress and outstanding problems," Statistics in Medicine, vol. 7, pp. 819-841, 1988.
[10] D.G. Clayton, "A Monte Carlo method for Bayesian inference in frailty models," Biometrics, vol. 47, pp. 467-485, 1991.
[11] D.G. Clayton and M. Hills, Statistical Models in Epidemiology.Oxford, England: Oxford Science Publications, Oxford Univ. Press, 1993.
[12] B.R. Cullis and C.A. McGilchrist, "A model for the analysis of growth data from designed experiments," Biometrics, vol. 46, pp. 131-142, 1990.
[13] A.P. Dawid, "Present position and potential developments: Some personal views. Statistical theory: The prequential approach (with discussion)," J. Royal Statistical Soc., Series A, vol. 147, pp. 178-292, 1984.
[14] P.J. Diggle, "An approach to the analysis of repeated measurements," Biometrics, vol. 44, pp. 959-971, 1988.
[15] T. Fearn, "A Bayesian approach to growth curves," Biometrika, 62, pp. 89-100, 1975.
[16] D. Gamerman and M. West, "An application of dynamic survival models in unemployment studies," The Statistician, vol. 36, pp. 269-274, 1987.
[17] D. Gamerman and H.S. Migon, "Dynamic hierarchical models," J. Royal Statistical Soc., Series B, vol. 55, no. 3, pp. 629-642, 1993.
[18] S. Geisser, "The predictive sample reuse method with application," J. Am. Statistical Assoc., vol. 70, pp. 320-328, 1975.
[19] A.E. Gelfand and A.F.M. Smith, "Sampling based approaches to calculating marginal densities," J. Am. Statistical Assoc., vol. 85, pp. 398-409, 1990.
[20] A.E. Gelfand, D.K. Dey, and H. Chang, "Model determination using predictive distributions with implementation via sampling-based methods," Bayesian Statistics 4, J.M. Bernardo, J.O. Berger, A.P. Dawid, and A.F.M. Smith, eds., pp. 147-167, Oxford Univ. Press, 1992.
[21] A.E. Gelfand, S.K. Sahu, and B.P. Carlin, "Efficient parametrizations for generalized linear mixed models," Proc. Fifth Int'l Meeting Bayesian Statistics, Alicante, pp. 48-74, 1994.
[22] A. Gelman and D.B. Rubin, "Inference from iterative simulation using multiple sequences," Statistical Science, vol. 7, pp. 457-472, 1992.
[23] S. Geman and D. Geman, "Stochastic relaxation, Gibbs distributions, and the Bayesian restoration of images," IEEE Trans. Pattern Analysis and Machine Intelligence, vol. 6, pp. 721-741, 1984.
[24] W.R. Gilks and P. Wild, "Adaptive rejection sampling for Gibbs sampling," Applied Statistics, vol. 41, pp. 337-348, 1992.
[25] W.R. Gilks, D.G. Clayton, D.J. Spiegelhalter, N.G. Best, A.J. McNeil, L.D. Sharples, and A.J. Kirby, "Modeling complexity: Applications of Gibbs sampling in medicine," J. Royal Statistical Soc., Series B, vol. 55, no. 1, pp. 39-52, 1993.
[26] W.R. Gilks, N.G. Best, and K.K.C. Tan, "Adaptive rejection Metropolis sampling for Gibbs sampling," technical report, Biostatistics Unit, Medical Research Council, Cambridge, U.K., 1994.
[27] K. Gordon and A.F.M. Smith, "Modeling and monitoring biomedical time series," J. Am. Statistical Assoc., vol. 85, no. 410, pp. 328-337, 1990.
[28] P.J. Green, "Penalized likelihood for general semiparametric regression models," Int'l Statistical Rev., vol. 55, pp. 245-259, 1987.
[29] P.J. Harrison and C.F. Stevens, "Bayesian forecasting (with discussion)," J. Royal Statistical Soc., Series B, vol. 38, pp. 205-247, 1976.
[30] W.K. Hastings, "Monte Carlo sampling methods using Markov chains and their applications," Biometrika, vol. 57, pp. 97-109, 1970.
[31] J.D. Kalbfleisch, "Nonparametric Bayesian analysis of survival time data," J. Royal Statistical Soc., Series B, vol. 40, pp. 214-221, 1978.
[32] A. Kong, J.S. Liu, and W.H. Wong, "Sequential imputations and Bayesian missing data problems," J. Am. Statistical Assoc., vol. 89, pp. 278-288, 1994.
[33] N. Lange, B.P. Carlin, and A.E. Gelfand, "Hierarchical Bayes models for the progression of HIV infection using longitudinal CD4 T-cell numbers," J. Am. Statistical Assoc., vol. 87, no. 419, pp. 615-632, 1992.
[34] S.L. Lauritzen, A.P. Dawid, B.N. Larsen, and H.G. Leimer, "Independence properties of directed Markov fields," Networks, vol. 20, pp. 491-505, 1990.
[35] S.L. Lauritzen and D.J. Spiegelhalter, "Local computations with probabilities on graphical structures and their application to expert systems," J. Royal Statistical Soc., Series B, vol. 50, no. 2, pp. 202-203, 1988.
[36] T. Leonard, "Density estimation, stochastic processes, and prior information," with discussion, J. Royal Statistical Soc., Series B, vol. 40, pp. 113-146, 1978.
[37] D.V. Lindley and A.F.M. Smith, "Bayes estimates for the linear model," J. Royal Statistical Soc., Series B, vol. 34, pp. 1-42, 1972.
[38] J.S. Liu and R. Chen, "Blind deconvolution via sequential imputations," J. Am. Statistical Assoc., vol. 90, to appear.
[39] N. Metropolis, A.W. Rosenbluth, M.N., Rosenbluth, A.H. Teller, and E. Teller, "Equations of state calculations by fast computing machine," J. Chemical Physics, vol. 21, pp. 1,087-1,091, 1953.
[40] H.S. Migon and D. Gamerman, "Generalized exponential growth models—A Bayesian approach," J. Forecasting, vol. 12, pp. 573-584, 1993.
[41] J. Pearl, "Evidential reasoning using stochastic simulation of causal models," Artificial Intelligence, vol. 32, pp. 245-257, 1987.
[42] S. Richardson and W.R. Gilks, "Conditional independence models for epidemiological studies with covariate measurement error," Statistics in Medicine, vol. 12, pp. 1,703-1,722, 1993.
[43] D.B. Rubin, "A noniterative sampling/importance resampling alternative to the data augmentation algorithm for creating a few imputations when fractions of missing information are modest: The SIR algorithm," comment on Tanner and Wong, "The calculation of posterior distributions by data augmentation, with discussion," J. Am. Statistical Assoc., vol. 52, pp. 543-546, 1987.
[44] A.F.M. Smith and M. West, "Monitoring renal transplants: An application of the multiprocess Kalman filter," Biometrics, vol. 39, pp. 867-878, 1983.
[45] A.F.M. Smith and G.O. Roberts, "Bayesian computation via the Gibbs sampler and related Markov chain Monte Carlo methods," J. Royal Statistical Soc., Series B, vol. 55, no. 1, pp. 3-23, 1993.
[46] J.M.G. Taylor, W.G. Cumberland, and J.P. Sy, "A Stochastic model for analysis of longitudinal AIDS data," J. Am. Statistical Soc., vol. 89, pp. 727-736, 1994.
[47] A. Thomas, D.J. Spiegelhalter, and W.R. Gilks, "BUGS: A program to perform Bayesian inference using Gibbs sampling," Bayesian Statistics 4, J.M. Bernardo, J.O. Berger, A.P. Dawid, and A.F.M. Smith, eds. Oxford, U.K.: Clarendon Press, pp. 837-842.
[48] L. Tierney, "Exploring posterior distributions using Markov chains," Computing Science and Statistics: Proc. 23rd Symp. Interface, pp. 563-570, 1991.
[49] J.W. Vaupel, K.G. Manton, and E. Stallard, "The impact of heterogeneity in individual frailty on the dynamics of mortality," Demography, vol. 16, pp. 439-454, 1979.
[50] J.C. Wakefield, A.F.M. Smith, A. Racine-Poon, and A.E. Gelfand, "Bayesian analysis of linear and nonlinear population models using the Gibbs sampler," Applied Statistics, to be published.
[51] M. West and P.J. Harrison, "Bayesian forecasting and dynamic models," Springer-Verlag, 1989.
[52] J. Whittaker, Graphical Models in Applied Multivariate Statistics.New York: John Wiley, 1990.
[53] M. West, P.J. Harrison, and H.S. Migon, "Dynamic generalized linear models and Bayesian forecasting," J. Am. Statistical Assoc., vol. 80, pp. 73-97, 1985.

Index Terms:
Bayesian inference, statistical forecasting, analysis of time series data, analysis of failure data, Markov chain Monte Carlo methods, conditional independence graphs, model determination, medical monitoring.
Citation:
Carlo Berzuini, Cristiana Larizza, "A Unified Approach for Modeling Longitudinal and Failure Time Data, With Application in Medical Monitoring," IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 18, no. 2, pp. 109-123, Feb. 1996, doi:10.1109/34.481537
Usage of this product signifies your acceptance of the Terms of Use.