Latent time‐varying factors in longitudinal analysis: a linear mixed hidden Markov model for heart rates

Longitudinal data are often segmented by unobserved time-varying factors, which introduce latent heterogeneity at the observation level, in addition to heterogeneity across subjects. We account for this latent structure by a linear mixed hidden Markov model. It integrates subject-specific random effects and Markovian sequences of time-varying effects in the linear predictor. We propose an expectationŰ-maximization algorithm for maximum likelihood estimation, based on data augmentation. It reduces to the iterative maximization of the expected value of a complete likelihood function, derived from an augmented dataset with case weights, alternated with weights updating. In a case study of the Survey on Stress Aging and Health in Russia, the model is exploited to estimate the influence of the observed covariates under unobserved time-varying factors, which affect the cardiovascular activity of each subject during the observation period.

[1]  Geoffrey J. McLachlan,et al.  Finite Mixture Models , 2019, Annual Review of Statistics and Its Application.

[2]  C. Cooper,et al.  Is grip strength a useful single marker of frailty? , 2003, Age and ageing.

[3]  A. Maruotti Mixed Hidden Markov Models for Longitudinal Data: An Overview , 2011 .

[4]  M. Malik,et al.  Sex differences in cardiac autonomic regulation and in repolarisation electrocardiography , 2013, Pflügers Archiv - European Journal of Physiology.

[5]  P. Stein,et al.  Insights from the study of heart rate variability. , 1999, Annual review of medicine.

[6]  J. Detilleux The analysis of disease biomarker data using a mixed hidden Markov model (Open Access publication) , 2008, Genetics, selection, evolution : GSE.

[7]  Maud Delattre Inference in Mixed Hidden Markov Models and Applications to Medical Studies , 2010 .

[8]  K. Sacco,et al.  Shared “Core” Areas between the Pain and Other Task-Related Networks , 2012, PloS one.

[9]  Jan Bulla,et al.  Computational issues in parameter estimation for stationary hidden Markov models , 2008, Comput. Stat..

[10]  Gary L Myers,et al.  Markers of inflammation and cardiovascular disease: application to clinical and public health practice: A statement for healthcare professionals from the Centers for Disease Control and Prevention and the American Heart Association. , 2003, Circulation.

[11]  G. Breithardt,et al.  Heart rate variability: standards of measurement, physiological interpretation and clinical use. Task Force of the European Society of Cardiology and the North American Society of Pacing and Electrophysiology. , 1996 .

[12]  X. Jouven,et al.  Gender-specific trends in heart rate in the general population from 1992–2007: a study of 226,288 French adults , 2013, European journal of preventive cardiology.

[13]  Clive W. J. Granger,et al.  Time Series Modelling and Interpretation , 1976 .

[14]  N. Unwin,et al.  Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol in Adults (Adult Treatment Panel III) Detection, Evaluation, and Treatment of High Blood Cholesterol Education Program (NCEP) Expert Panel on Executive Summary of the Third Report of the National , 2009 .

[15]  J. Nash Compact Numerical Methods for Computers , 2018 .

[16]  H. Tunstall-Pedoe,et al.  Estimation of ten-year risk of fatal cardiovascular disease in Europe: the SCORE project. , 2003, European heart journal.

[17]  Salvatore Ingrassia,et al.  Degeneracy of the EM algorithm for the MLE of multivariate Gaussian mixtures and dynamic constraints , 2011, Comput. Stat. Data Anal..

[18]  J. Mckenney,et al.  Executive Summary of The Third Report of The National Cholesterol Education Program (NCEP) Expert Panel on Detection, Evaluation, And Treatment of High Blood Cholesterol In Adults (Adult Treatment Panel III). , 2001, JAMA.

[19]  C L Feldman,et al.  Determinants of heart rate variability. , 1996, Journal of the American College of Cardiology.

[20]  D. Nunan,et al.  A Quantitative Systematic Review of Normal Values for Short‐Term Heart Rate Variability in Healthy Adults , 2010, Pacing and clinical electrophysiology : PACE.

[21]  A. Maruotti,et al.  A Multivariate Hidden Markov Model for the Identification of Sea Regimes from Incomplete Skewed and Circular Time Series , 2012 .

[22]  Christophe Biernacki,et al.  Choosing starting values for the EM algorithm for getting the highest likelihood in multivariate Gaussian mixture models , 2003, Comput. Stat. Data Anal..

[23]  A. Malliani,et al.  Heart rate variability. Standards of measurement, physiological interpretation, and clinical use , 1996 .

[24]  Giuseppe Boccignone,et al.  Gaussian Mixture Model of Heart Rate Variability , 2012, PloS one.

[25]  Eric Moulines,et al.  Inference in hidden Markov models , 2010, Springer series in statistics.

[26]  J. Vermunt,et al.  Latent class and finite mixture models for multilevel data sets , 2008, Statistical methods in medical research.

[27]  W. Zucchini,et al.  Hidden Markov Models for Time Series: An Introduction Using R , 2009 .

[28]  G. Celeux,et al.  Mixture of linear mixed models for clustering gene expression profiles from repeated microarray experiments , 2005 .

[29]  L. Køber,et al.  Resting, night-time, and 24 h heart rate as markers of cardiovascular risk in middle-aged and elderly men and women with no apparent heart disease. , 2013, European heart journal.

[30]  F. Lagona,et al.  Model-based clustering of multivariate skew data with circular components and missing values , 2012 .

[31]  Dylan S. Small,et al.  HIDDEN MARKOV MODELS FOR ALCOHOLISM TREATMENT TRIAL DATA , 2010, 1010.1410.

[32]  J. Cockcroft,et al.  Left ventricular mechanics in humans with high aerobic fitness: adaptation independent of structural remodelling, arterial haemodynamics and heart rate , 2012, The Journal of physiology.

[33]  Roland Langrock,et al.  Combining hidden Markov models for comparing the dynamics of multiple sleep electroencephalograms , 2013, Statistics in medicine.

[34]  Gérard Govaert,et al.  Assessing a Mixture Model for Clustering with the Integrated Completed Likelihood , 2000, IEEE Trans. Pattern Anal. Mach. Intell..

[35]  F. Lagona,et al.  Maximum likelihood estimation of bivariate circular hidden Markov models from incomplete data , 2013 .

[36]  Gilles Celeux,et al.  Combining Mixture Components for Clustering , 2010, Journal of computational and graphical statistics : a joint publication of American Statistical Association, Institute of Mathematical Statistics, Interface Foundation of North America.

[37]  R. Altman Mixed Hidden Markov Models , 2007 .

[38]  Dmitri A. Jdanov,et al.  Biological mechanisms of disease and death in Moscow: rationale and design of the survey on Stress Aging and Health in Russia (SAHR) , 2009, BMC public health.

[39]  New York Dover,et al.  ON THE CONVERGENCE PROPERTIES OF THE EM ALGORITHM , 1983 .

[40]  Peter E. Rossi,et al.  Case Studies in Bayesian Statistics , 1998 .