- Research article
- Open Access
- Open Peer Review
Progression of liver cirrhosis to HCC: an application of hidden Markov model
BMC Medical Research Methodology volume 11, Article number: 38 (2011)
Health service databases of administrative type can be a useful tool for the study of progression of a disease, but the data reported in such sources could be affected by misclassifications of some patients' real disease states at the time. Aim of this work was to estimate the transition probabilities through the different degenerative phases of liver cirrhosis using health service databases.
We employed a hidden Markov model to determine the transition probabilities between two states, and of misclassification. The covariates inserted in the model were sex, age, the presence of comorbidities correlated with alcohol abuse, the presence of diagnosis codes indicating hepatitis C virus infection, and the Charlson Index. The analysis was conducted in patients presumed to have suffered the onset of cirrhosis in 2000, observing the disease evolution and, if applicable, death up to the end of the year 2006.
The incidence of hepatocellular carcinoma (HCC) in cirrhotic patients was 1.5% per year. The probability of developing HCC is higher in males (OR = 2.217) and patients over 65 (OR = 1.547); over 65-year-olds have a greater probability of death both while still suffering from cirrhosis (OR = 2.379) and if they have developed HCC (OR = 1.410). A more severe casemix affects the transition from HCC to death (OR = 1.714). The probability of misclassifying subjects with HCC as exclusively affected by liver cirrhosis is 14.08%.
The hidden Markov model allowing for misclassification is well suited to analyses of health service databases, since it is able to capture bias due to the fact that the quality and accuracy of the available information are not always optimal. The probability of evolution of a cirrhotic subject to HCC depends on sex and age class, while hepatitis C virus infection and comorbidities correlated with alcohol abuse do not seem to have an influence.
The evolution of chronic degenerative disease is characterized by progression through intermediate states to advanced disease and death. For these diseases, survival analysis must take into account the various transitions from one state to the next, as well as a series of prognostic variables that can have an influence on each event including death. For example, liver cirrhosis can evolve to hepatocellular carcinoma, and the presence of comorbidities, exposure to hepatitis B or C virus, as well as alcohol consumption and age, can influence the terminal event. In fact, liver cirrhosis is well known to consist of a diffuse alteration of the liver structure resulting from protracted processes of liver inflammation and necrosis of different natures. The main causes of liver cirrhosis are chronic viral hepatitis B or C and the consumption of alcohol. In particular, alcohol abuse can halve the time of onset of cirrhosis in a patient already affected by chronic viral hepatitis (from about 20-30 years to 10-15 years). Hepatocellular carcinoma occurs at a rate of 1% to 4% per year after cirrhosis is established  and cirrhosis underlies HCC in approximately 80%-90% of cases worldwide . Stochastic multistate or competing models, like Markov chains, are those best suited to the analysis of such phenomena [3–6]. It would, of course, be too long and costly to program clinical studies, or indeed prospective and follow-up trials, to study the natural history of chronic degenerative diseases in order to be able to apply multistate models correctly. In fact, in patients affected by these diseases passage from one state to the next often occurs after fairly long intervals.
Retrospective studies, based on the use of health service databases of administrative type (HSDBA), can be a valid alternative, despite the limits posed by the fact that the quality and accuracy of the available information are not always optimal. For example, using the hospital discharge sheets (HDS) database, patients can be followed from the probable diagnosis of onset of the disease through the subsequent worsening states. Moreover, by means of linkage with the death certificates database (DCDB), it is ultimately possible to trace the cause of death, if applicable.
A common problem when applying Markov models to HSDBA is that assessment of the disease state of an individual can be subject to classification errors. A similar classification problem frequently occurs in research in the field of social sciences, where wide use is made of a series of models that take into account latent states [7–9].
Multistate models that exploit the properties of Markov chains offer a useful methodological structure for describing complex time-dependent outcomes . The procedure for estimation of the probability of transition most widely adopted in Markov models is the Cox proportional model. This model can describe survival time in function of a multitude of prognostic factors , under the fundamental assumption of the proportionality of hazards, in other words that the examined factors will have a constant impact over time on the risk of death.
Jackson et al. (2003)  described a procedure for simultaneously estimating the transition rates and the probabilities of misclassification in a hidden Markov model, supplying software for implementing the multistate hidden Markov model in the R Project programming environment. Aim of the present work is to study the pathway leading subjects affected by liver cirrhosis to develop hepatocellular carcinoma and to death, determining the transition probabilities through the degenerative disease states and verifying whether, apart from being risk factors for the onset of cirrhosis, chronic viral hepatitis B and C infection and alcohol abuse also have a role in the process leading cirrhotic subjects to develop hepatocellular carcinoma and/or to death.
Retrospective observational studies can be conducted using the HSDBA, to assess the natural history of the disease in a group of subjects. In fact, the HDS database can be employed to individuate all subjects admitted to hospital one or more times, as well as the duration of the interval between one hospitalization and the next. For each hospitalization, the disease state observed and the covariates of interest can be traced. However. the disease state may be affected by an unknown degree of error, due both to an imperfect diagnosis and to incorrect classification. Moreover, due to the irregular nature of follow-up, observations of the actual time of entry into a disease state are frequently interval censored. The proposed Markov model allows us to take into account the above characteristics.
In a sample of n subjects, for the ith subject we assume that the following variables are observed at the jth visit:
Tij chronological time of clinical visit
Yij binary disease outcome measurement
Zij vector of covariates.
The observed disease outcome measures Yij are subject to error. The actual underlying disease status is assumed to be a process evolving in continuous time and is denoted (X(t), t > 0). This process is unobserved or 'hidden', and will be modeled as a continuous time two-state Markov process, where the states are interpreted as the presence or absence of the disease manifestation. Let Y1 j and T1 j denote the sequence from 1 to j of observed disease states and observation times for an individual i. The Markov assumption for the hidden disease process is given by
where the quantity Pxj-1, xj denotes the probability of transition to occupy state x at time Tj = tj given that the process was in state xj-1 at tj-1 and that the transition probabilities are assumed to be stationary. We also assume that, conditional on the state of the hidden process at time tj, an observation Yj is independent of all previous observations and the hidden process prior to time tj:
When Y is binary, f(y|x) can be interpreted as the probability of correctly or incorrectly classifying the disease state given the true state of the subject. The conditional independence of misclassification probabilities at successive time points is a strong assumption. However, dependence between successive measurements is confounded with dependence between the true disease states. It is not possible to disentangle these two features of the model with misclassified data .
If the disease status is observed accurately, then X(Tj) and Yj coincide, and the model is reduced to a pure continuous time Markov process.
At any time t and for each pair of states r and s, the transition from one state to the next and the time when this transition occurs are regulated by the transition intensity qrs. The transition intensity qrs represents the instantaneous risk of moving from state r to state s:
The intensities form a matrix Q whose rows sum to zero, so that the diagonal entries are defined by .
On the basis of the available data it is necessary to estimate the matrix Q of transition intensity according to the method described by Kalbfleisch and Lawless , and Kay . In addition, at each level of the model explicative variables can be included and it is possible to use a proportional hazards model to relate transition intensities qrs(t) to time t with the covariates z(t),
The new matrix of the transition probabilities P(t), that takes into account any necessary covariates, can be calculated by taking the matrix exponential of the scaled transition intensity matrix (see, for example, Cox and Miller ):
The hazards proportionality can be verified using the Schoenfeld residuals of the model [17, 18], defined as the value x ik of the covariate K for individual i who actually died at time ti minus the expected value, where the expected value is given by and p i is the probability of death of individual i at time ti. The graph, that represents the trend of the Schoenfeld residuals calculated for each individual and each covariate, can be used to directly visualize the hazards ratio . Assuming proportionality of the hazards, the Schoenfeld residuals are independent of time. The presence of a linear relationship with time, an index of non proportionality, can be tested by performing a simple linear regression and a trend test: a slope significantly different from zero would be evidence against proportionality and an increasing (decreasing) trend would indicate an increasing (decreasing) hazards ratio over time.
In our model the estimation of the parameters was obtained using the maximum likelihood method (MLE), and assuming that the chance variables are independent and normally distributed. We considered a two-state disease model: cirrhosis and hepatocellular carcinoma and the absorbent state, death, that is irreversible (Figure 1). State 1 is that of patients with a diagnosis of liver cirrhosis (codes ICD9cm 571.2 and 571.5); State 2 is that of patients with a diagnosis of hepatocellular carcinoma (codes ICD9cm 155.0 and 155.2), regardless of the concomitant presence of a diagnosis of cirrhosis; State 3, absorbent, is that of deceased subjects. The few cases of patients who underwent liver surgery, including transplant, were excluded (codes ICD9cm for procedures 50.21÷50.99). Supposing, therefore, that progression of cirrhosis to HCC is an irreversible procedure, but in any case taking into account the possibility of misclassification between the first and second state and vice versa, we estimated the following transition intensities (Q) and misclassification (E) matrices:
The covariates inserted in the model were age, sex, the presence of disease correlated to alcohol abuse, the presence of diagnostic codes correlated to hepatitis C virus infection and the Charlson Index assessing the subject's clinical severity.
The Charlson index, developed in 1987  and adapted to health data banks by Deyo et al. , is based on ICD9-CM diagnosis codes and contains 17 categories of comorbidity, each with an associated weight ranging from 1 to 6; the sum of all the weights gives the value of the index that, being determined in this way, takes into account both the number of comorbidities and their severity. The Charlson Index was divided into 2 classes: less than or equal to 3, greater than 3, as also was age: less than or equal to 65 years, older than 65 years. We selected 65 as the age cut-off because this is considered in the literature to be the mean age for development of liver cancer [21, 22] and also because it corresponded to the median value of the observed distribution of cases of HCC, guaranteeing the stability of the model.
The analysis was conducting using the electronic Hospital Discharge Sheet coming from all the Apulia hospitals for the years 2000-2006, and Death Certificate DataBase related to all deaths events in Apulia for the years 2000-2006. From both databases were selected cases with the following ICD-9-CM (International Statistical Classification of Diseases and Related Health Problems 9threvision Clinical Modification) as principal or secondary diagnosis:
Liver cirrhosis (LC): 571.2; 571.5
Hepatocellular carcinoma (HCC): 155.0; 155.2.
To establish the starting point of the first state, defined as of the first hospitalization with LC diagnosis, were selected patients in whom the onset of cirrhosis presumably occurred between 01/01/2000 and 31/12/2000. For this purpose, we eliminated patients who had been hospitalized at least once for cirrhosis or hepatocellular carcinoma in the years 1998 and 1999. After identification the clinical course of the selected cirrhotic patients was reconstructed searching others hospitalizations with diagnosis of LC or HCC in HDS databases, or death correlated with LC or HCC in HDS and DCDB databases.
In total, 1925 patients were identified who had a presumed onset of cirrhosis in 2000 and had undergone at least one transition of state by 31/12/2006 (Table 1). In 33 patients, after a hospital admission for hepatocellular carcinoma, a subsequent hospitalization occurred with a diagnosis of only cirrhosis despite not having undergone surgery. These subjects were considered to have been misclassified. Table 2 shows the frequency of the covariates considered in the model. Estimates of the mean time of persistence in each state, calculated according to the procedure established by Jackson et al. [12, 23], are reported in Table 3. We calculated the times of persistence also for the different levels of the covariates. The shortest time of persistence in state 1 (months = 24.99; CI =18.35-34.04) was observed in male subjects aged > 65 years, with a Charlson Index > 3 and free from diseases correlated to alcohol abuse or hepatitis C virus. The longest estimated time (months = 178.92; CI = 151.13-211.82) was observed in female subjects aged ≤ 65 years, with a Charlson Index ≤ 3, the presence of disease correlated to alcohol abuse and absence of hepatitis C virus infection. The shortest time of persistence in state 2 (months = 15.63; CI = 1.95-125.38) was demonstrated in male subjects aged > 65 years, with a Charlson Index > 3 and the presence of diseases correlated to alcohol abuse and hepatitis C virus infection, whereas the longest estimated time (months = 48.91; CI = 37.99-62.97) was found in male patients aged ≤ 65 years, with a Charlson Index ≤ 3, and absence of diseases correlated either to alcohol abuse or hepatitis C virus infection. For each covariate, the proportionality of hazards was verified using the Schoenfeld residuals method (Table 4). The overall test shows strong evidence of the proportionality of hazards and all the variables contribute to this proportionality.
The parameters estimated for the hidden Markov model are reported in Table 5. The estimated intensity matrix demonstrates that cirrhotic patients have twice the probability (0.0151/0.0071) of developing a liver cancer than of dying without developing a tumour. Moreover, the probability of death is four-fold higher in a subject with a liver cancer than in a subject with only cirrhosis (0.0284/0.0071).
The estimated odds ratios are reported in Table 6.
Figure 2 shows the trend over time of the estimated probabilities of transition in cirrhotic subjects free from comorbidities correlated to alcohol abuse or HCV infection. The probabilities of transition in Figures 2a and 2c are referred to cirrhotic subjects, male and female, aged 65 years or younger and with a Charlson Index of 3 or less.
The probabilities of transition from one state to the next are generally higher in males; in particular, the probability of progressing to HCC is never higher than 15.7% in females, whereas it reaches 30% in males after 4 years. In the first 40 months after the onset of liver cirrhosis, the probability of developing hepatocellular carcinoma is greater than the probability of death in males: 15.7% versus 6.8% after one year, 24.5% versus 15.5% after two years. In females, instead, the probability of death is little lower than that of developing HCC in the first few months, while already after two years the probability of death is higher (13.4% versus 11.8%).
The male-female differences in the probabilities of transition from cirrhosis to death grow over the first 3 years (being 10% at 36 months) but remain constant thereafter. The probability of developing HCC in the short term is higher in males than females (after 2 years, it is 25.8% in males versus 13.2% in females), but the values tend to converge in the longer term (6.3% in males versus 5.5% in females after 7 years).
The probabilities of death, both for subjects with only cirrhosis and with HCC, are generally higher than in subjects with a less severe casemix. In females the casemix does not seem to affect the probability of transition cirrhosis-HCC, whereas in males the more severe casemix has a long term effect when the probability of developing HCC declines in favor of a greater probability of death.
By exploiting the properties of Markov chains applied to a stochastic multistate model, we have calculated the temporal intensities of transition during the degenerative course of chronic liver cirrhosis. It was also possible to determine the time hazards of degeneration of the liver disease until death. Various studies have employed Cox models with the principal aim of determining the risk factors for the progression of cirrhosis to HCC [24, 25]. In 2000, Degos et al.  studied progression to HCC and death using an "illness-disease-death" in a small cohort of subjects with a diagnosis of HCV-related cirrhosis, estimating the time of the events with the Kaplan Meier method. Then, in 2007 Ioannou et al.  applied the Cox proportional hazards model to a large administrative database of cirrhotic patients to determine the incidence of hepatocellular carcinoma, but without calculating the probabilities of transition from the first to the second state. Moreover, they acknowledged that a limit of their study was the risk of error when filling out the diagnosis codes ICD-9. In our study, using the data contained in the administrative health services databases, the multistate hidden Markov model we applied enabled us to identify the possible misclassification errors that can occur.
Previous studies based on surveillance programs for hepatocellular carcinoma in cirrhotic subjects reported an incidence of hepatocellular carcinoma of 1.5, 2.5, 6.7 per 100 subjects a year. In our study we found an incidence of HCC of 1.5% per year. Ioannou et al. (2007) , who conducted their analysis on an administrative database as we have done, found an incidence of 2.4% per year. The difference between these two values (1.5% vs 2.4%) is probably due to the different reference population, but also to insertion in our model of the misclassification matrix. In fact, the probability of misclassifying subjects with HCC as subjects with cirrhosis alone was revealed to be 14.08%, while the reverse error, i.e. misclassifying cirrhotic subjects as affected by HCC, was 2.73%.
The estimated parameters for the covariate "age class" were statistically significant for all the transitions. "Sex" was significant only for the transition of state from cirrhosis to HCC, while the Charlson Index had an effect only on the transition from HCC to death.
In various longitudinal studies it has been shown that advanced age and the male sex are associated with an increased risk of HCC in cirrhotic subjects [24, 31–33], as was also shown in the present study. In fact, male cirrhotic subjects have approximately twice the probability of developing hepatocellular carcinoma as compared to female cirrhotic patients, while elderly patients (aged ≥ 65 years) have a higher risk of degeneration of the liver disease, and especially of dying while affected by cirrhosis.
The ample presence of concomitant diseases (Charlson Index ≥ 3) increases the risk of death in subjects with HCC. The Charlson Index was not found to have an incidence on the transition from cirrhosis to HCC. The insertion of specific comorbidities among the model covariates, such as diabetes mellitus [27, 34, 35], could help to identify the co-morbid conditions that may become risk factors for progression to HCC in cirrhotic subjects.
Alcohol is proposed to cause HCC mainly because it causes cirrhosis, whereas its association with HCC without cirrhosis is controversial and it probably has no direct carcinogenic role . In a case-control study conducted in 2007, Kumar et al.  showed that although HCV RNA positivity and alcohol abuse significantly increased the risk of hepatocellular carcinoma among cirrhotic patients, no significant risk increase was evident in the absence of cirrhosis. In our study the estimates of the parameters related to hepatitis C virus and to the presence of morbidity correlated to alcohol abuse did not result significant. This result likely confirms that hepatitis C virus and alcohol abuse are risk factors for the onset of cirrhosis , but once the cirrhosis has become established they do not have an influence on the development of HCC.
Despite the advantages of the structured hidden Markov model approach, there are some limitations. As is typical of population dynamic models, collinearity between parameter estimates can lead to identifiability problems . As a consequence, maximum likelihood estimates can sometimes yield implausible parameter values, and maximization algorithms may fail to converge . The problem becomes particularly severe when time series are short and data scarce. Moreover, addition to the model of further covariates could exacerbate this problem. However, it seems likely that such problems with the proposed structured hidden Markov models might be overcome by adopting a Bayesian formulation.
Markov model proved to be a useful tool for analysis of the course of a chronic degenerative disease like liver cirrhosis. In particular, the hidden Markov model that takes into account the risk of misclassification is well suited to the analysis of administrative health data because it can capture bias due to the problem that the data quality is not always optimal, as well as enabling the study of the effect of different covariates on the transitions of state.
A further improvement of the model could be that of predicting the state of subjects who have undergone surgery and appear cured or temporarily HCC-free. In our case, due to the paucity of such observations, it was not possible to consider a model like this.
Di Bisceglie AM: Hepatitis C and hepatocellular carcinoma. Hepatology. 1997, 26 (3 suppl 1): 34S-38S. 10.1002/hep.510260706.
Fattovich G, Stroffolini T, Zagni I, Donato F: Hepatocellular carcinoma in cirrhosis: incidence and risk factors. Gastroenterology. 2004, 127: S35-S50. 10.1053/j.gastro.2004.09.014.
Putter H, Fiocco M, Geskus RB: Tutorial in biostatistics: competing risks and multi-state models. Statistics in Medicine. 2007, 26: 2389-2430. 10.1002/sim.2712.
Serio G, Morabito A: Considerations on a staging process for chronic diseases. Rivista di Statistica Applicata. 1988, 21 (23): 335-348.
Serio G, Morabito A: Stochastic survival model with covariates in cancer. Modelling of Biomedical Systems. 1986, 91-96.
Chen HH, Duffy SW, Tabar L: A Markov chain method to estimate the tumour progression rate from preclinical to clinical phase, sensitivity and positive predictive value for mammography in breast cancer screening. The Statistician. 1996, 45 (3): 307-317. 10.2307/2988469.
Jackson CH, Sharples LD: Hidden Markov models for the onset and progression of bronchiolitis obliterans syndrome in lung transplant recipients. Statistics in Medicine. 2002, 21: 113-128. 10.1002/sim.886.
Clogg CC: Latent class models: recent developments and prospects for the future. Handbook of Statistical Modeling in the Social Sciences. Edited by: Arminger G, Clogg CC, Sobel ME. 1995, Plenum, New York, 311-359.
MacDonald IL, Zucchini W: Hidden Markov Models and Other Models for Discrete-valued Time Series. 1997, Chapman and Hall: New York
Scheike TH, Zhang M: Direct Modelling of Regression Effects for Transition Probabilities in Multistate Models. Scandinavian Journal of Statistics. 2007, 34: 17-32. 10.1111/j.1467-9469.2006.00544.x.
Cox D: Regression Models and Life-Tables. Journal of the Royal Statistical Society, Series B. 1972, 34: 187-220.
Jackson CH, Sharples LD, Thompson SG, Duffy SW, Couto E: Multistate Markov models for disease progression with classification error. Journal of the Royal Statistical Society, Series D - The Statistician. 2003, 52 (2): 193-209. 10.1111/1467-9884.00351.
Bureau A, Shiboski S, Hughes JP: Applications of continuous time hidden Markov models to the study of misclassified disease outcomes. Statistics in Medicine. 2003, 22: 441-462. 10.1002/sim.1270.
Kalbfleisch JD, Lawless JF: The analysis of panel data under a Markov assumption. Journal of the American Statistical Association. 1985, 80 (392): 863-871. 10.2307/2288545.
Kay R: A Markov model for analysing cancer markers and disease states in survival studies. Biometrics. 1986, 42: 855-865. 10.2307/2530699.
Cox DR, Miller HD: The Theory of Stochastic Processes. 1965, Chapman and Hall, London
Grambsch P, Therneau T: Proportional Hazards Tests and Diagnostics Based on Weighted Residuals. Biometrika. 1994, 81: 515-26. 10.1093/biomet/81.3.515.
Schoenfeld D: Chi-squared goodness of fit test for the proportional hazards regression model. Biometrika. 1981, 67: 147-53.
Charlson ME, Pompei P, Ales KL, MacKenzie CR: A new method of classifying prognostic comorbidity in longitudinal studies: development and validation. J Chron Dis. 1987, 40: 373-383. 10.1016/0021-9681(87)90171-8.
Deyo RA, Cherkin DC, Ciol MA: Adapting a clinical comorbidity index for use with ICD-9- CM administrative databases. J Clin Epidemiol. 1992, 45: 613-619. 10.1016/0895-4356(92)90133-8.
Ryder SD: Guidelines for the diagnosis and treatment of hepatocellular carcinoma (HCC) in adults. Gut. 2003, 52 (Suppl. III): iii1-8.
Hepatic Carcinoma, Primary Author: Keith E Stuart, MD, Chairman, Department of Hematology and Oncology, Lahey Clinic Coauthor(s): Zsofia K Stadler, MD. [http://emedicine.medscape.com/article/282814-overview]
Duffy SW, Chen HH: Estimation of mean sojourn time in breast cancer screening using a Markov chain model of entry to and exit from preclinical detectable phase. Statistics in Medicine. 1995, 14: 1531-1543. 10.1002/sim.4780141404.
Zaman SN, Melia WM, Johnson RD, Portmann BC, Johnson PJ, Williams R: Risk factors in development of hepatocellular carcinoma in cirrhosis: prospective study of 613 patients. Lancet. 1985, 1: 1357-60. 10.1016/S0140-6736(85)91785-4.
Sarbah SA, Gramlich T, Younoszai A, et al: Risk factors for hepatocellular carcinoma in patients with cirrhosis. Dig Dis Sci. 2004, 49: 850-853. 10.1023/B:DDAS.0000030098.75759.32.
Degos F, Christidis C, Ganne-Carrie N, Farmachidi JP, Degott C, Guettier C, Trinchet JC, Beaugrand M, Chevret S: Hepatitis C virus related cirrhosis: time to occurrence of hepatocellular carcinoma and death. Gut. 2000, 47: 131-136. 10.1136/gut.47.1.131.
Ioannou GN, Splan MF, Weiss NS, McDonald GB, Beretta L, Lee SP: Incidence and predictors of hepatocellular carcinoma in patients with cirrhosis. Clin Gastroenterol Hepatol. 2007, 5: 938-945. 10.1016/j.cgh.2007.02.039. 945.e1-945.e4
Fattovich G, Giustina G, Degos F, et al: Morbidity and mortality in compensated cirrhosis type C: a retrospective follow-up study of 384 patients. Gastroenterology. 1997, 112: 463-472. 10.1053/gast.1997.v112.pm9024300.
Colombo M, de Franchis R, Del Ninno E, et al: Hepatocellular carcinoma in Italian patients with cirrhosis. N Engl J Med. 1991, 325: 675-680. 10.1056/NEJM199109053251002.
N'Kontchou G, Paries J, Htar MT, et al: Risk factors for hepatocellular carcinoma in patients with alcoholic or viral C cirrhosis. Clin Gastroenterol Hepatol. 2006, 4: 1062-1068.
Sangiovanni A, Del Ninno E, Fasani P, De Fazio C, Ronchi G, Romeo R, Morabito A, de Franchis R, Colombo M: Increased survival of cirrhotic patients with a hepatocellular carcinoma detected during surveillance. Gastroenterology. 2004, 126: 1005-1014. 10.1053/j.gastro.2003.12.049.
Benvegnù L, Chemello L, Noventa F, Fattovich G, Pontisso P, Alberti A: Retrospective analysis of the effect of interferon therapy on the clinical outcome of patients with viral cirrhosis. Cancer. 1998, 83: 901-909.
Bolondi L, Sofia S, Siringo S, Gaiani S, Casali A, Zironi G, Piscaglia F, Gramantieri L, Zanetti M, Sherman M: Surveillance programme of cirrhotic patients for early diagnosis and treatment of hepatocellular carcinoma: a cost effectiveness analysis. Gut. 2001, 48: 251-259. 10.1136/gut.48.2.251.
La Vecchia C, Negri E, Decarli A, et al: Diabetes mellitus and the risk of primary liver cancer. Int J Cancer. 1997, 73: 204-207. 10.1002/(SICI)1097-0215(19971009)73:2<204::AID-IJC7>3.0.CO;2-#.
El-Serag HB, Tran T, Everhart JE: Diabetes increases the risk of chronic liver disease and hepatocellular carcinoma. Gastroenterology. 2004, 126: 460-468. 10.1053/j.gastro.2003.10.065.
London WT, McGlynn KA: Liver cancer. Cancer Epidemiology and Prevention. Edited by: Schottenfeld D, Fraumeni JF Jr. 1996, New York: Oxford University Press, 772-93.
Kumar M, Kumar R, Hissar SS, Saraswat MK, Sharma BC, Sakhuja P, Sarin SK: Risk factors analysis for hepatocellular carcinoma in patients with and without cirrhosis: a case -control study of 213 hepatocellular carcinoma patients from India. J Gastroenterol Hepatol. 2007, 22 (7): 1104-11. 10.1111/j.1440-1746.2007.04908.x.
Pradat P, Voirin N, Tillmann HL, Chevallier M, Trepo C: Progression to cirrhosis in hepatitis C patients: an age-dependent process. Liver Int. 2007, 27 (3): 335-339. 10.1111/j.1478-3231.2006.01430.x.
Brookhart MA, Hubbard AE, Van Der Laan MJ, Colford JM, Eisenberg JN: Statistical estimation of parameters in a disease transmission model: analysis of a Cryptosporidium outbreak. Statistics in Medicine. 2002, 21: 3627-38. 10.1002/sim.1258.
Cooper B, Lipsitch M: The analysis of hospital infection data using hidden Markov models. Biostatistics. 2004, 5 (2): 223-237. 10.1093/biostatistics/5.2.223.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2288/11/38/prepub
The authors declare that they have no competing interests.
NB conceived the study, conducted the analysis, wrote the manuscript. PT conceived the study, collaborated in the analysis, drafted the manuscript. GS supervised the analysis and reviewed the manuscript. All authors read and approved the final manuscript.
Nicola Bartolomeo, Paolo Trerotoli and Gabriella Serio contributed equally to this work.