A multistate approach for estimating the incidence of human immunodeficiency virus by using HIV and AIDS French surveillance data

Before the introduction of highly active antiretroviral therapy (HAART) in 1996, monitoring the epidemics of HIV infection was extensively based on back-calculating the number of new HIV infections from AIDS surveillance data and knowledge of the AIDS incubation period distribution. However, the increase in the AIDS incubation period induced by HAART has complicated the use of AIDS diagnosis data only, and made it necessary to supplement these data with HIV diagnosis data. We explore the advantage of combining HIV and AIDS surveillance data in the context where the HIV diagnosis data are available only for the most recent years. Extending the earlier work of Aalen et al. (Statist. Med. 1997; 16(19):2191-2210) based on a discrete-time Markov model that describes simultaneously disease progression, HIV diagnosis and treatment intake, we propose a penalized likelihood approach to estimate smooth HIV incidence, together with HIV diagnosis rates. The smoothing parameter is chosen using an approximated cross-validation criterion. In a simulation study, we show that incorporation of HIV test information improves the precision of the estimation of the incidence of infections in recent periods. The method is illustrated using HIV and AIDS surveillance data collected by the Institut de Veille Sanitaire, France.

[1]  I. Marschner,et al.  Using time of first positive HIV test and other auxiliary data in back-projection of AIDS incidence. , 1994, Statistics in medicine.

[2]  R. Salamon,et al.  [Human immunodeficiency virus infection and AIDS in Aquitaine. 10 years' experience of a hospital information system, 1985-1995. Le Groupe d'Epidemiologie Clinique du SIDA en Aquitaine (GECSA)]. , 1997, Presse medicale.

[3]  Rodolphe Thiébaut,et al.  A multistate approach for estimating the incidence of human immunodeficiency virus by using data from a prevalent cohort study , 2005 .

[4]  G. Wahba.  Bayesian "Confidence Intervals" for the Cross-validated Smoothing Spline , 1983 .

[5]  B. Yandell,et al.  Automatic Smoothing of Regression Functions in Generalized Linear Models , 1986 .

[6]  D. Costagliola,et al.  Including pre-AIDS mortality in back-calculation model to estimate HIV prevalence in France, 2000 , 2006, European Journal of Epidemiology.

[7]  D Commenges,et al.  A penalized likelihood approach for arbitrarily censored and truncated data: application to age-specific incidence of dementia. , 1998, Biometrics.

[8]  V. De Gruttola,et al.  Nonparametric analysis of truncated survival data, with application to AIDS , 1988 .

[9]  D. Marquardt.  An Algorithm for Least-Squares Estimation of Nonlinear Parameters , 1963 .

[10]  A Alioum,et al.  A direct approach for correcting AIDS incidence: variance formula and comparison with other methods. , 1995, Statistics in medicine.

[11]  John W. Ward,et al.  1993 revised classification system for HIV infection and expanded surveillance case definition for AIDS among adolescents and adults. , 1993, MMWR. Recommendations and reports : Morbidity and mortality weekly report. Recommendations and reports.

[12]  A. Senthilselvan.  Penalized Likelihood Estimation of Hazard and Intensity Functions , 1987 .

[13]  Odd O Aalen,et al.  Bayesian back‐calculation using a multi‐state model with application to HIV , 2005, Statistics in medicine.

[14]  N E Day,et al.  A Markov model for HIV disease progression including the effect of HIV diagnosis and treatment: application to AIDS prediction in England and Wales. , 1997, Statistics in medicine.

[15]  Mitchell H. Gail,et al.  A Method for Obtaining Short-Term Projections and Lower Bounds on the Size of the AIDS Epidemic , 1988 .

[16]  R Bellocco,et al.  Joint analysis of HIV and AIDS surveillance data in back-calculation. , 2000, Statistics in medicine.

[17]  Benoit Liquet,et al.  Choice between Semi‐parametric Estimators of Markov and Non‐Markov Multi‐state Models from Coarsened Observations , 2007 .