 Research
 Open Access
 Published:
Markov modeling in hepatitis B screening and linkage to care
Theoretical Biology and Medical Modelling volumeÂ 14, ArticleÂ number:Â 11 (2017)
Abstract
Background
With up to 240 million people chronically infected with hepatitis B worldwide, including an estimated 2 million in the United States, widespread screening is needed to link the infected to care and decrease the possible consequences of untreated infection, including liver cancer, cirrhosis and death. Screening is currently fraught with challenges in both the developed and developing world. New pointofcare tests may have advantages over standardofcare tests in terms of costeffectiveness and linkage to care. Stochastic modeling is applied here for relative utility assessment of pointofcare tests and standardofcare tests for screening.
Methods
We analyzed effects of pointofcare versus standardofcare testing using Markov models for disease progression in individual patients. Simulations of large cohorts with distinctly quantified models permitted the assessment of particular screening schemes. The validity of the trends observed is supported by sensitivity analyses for the simulation parameters.
Results
Increased utilization of pointofcare screening was shown to decrease hepatitis Brelated mortalities and increase life expectancy at low projected expense.
Conclusions
The results suggest that standardofcare screening should be substituted by pointofcare tests resulting in improved linkage to care and decrease in longterm complications.
Background
With up to 240 million people chronically infected with hepatitis B virus (HBV) worldwide [1], including an estimated 2 million people in the United States [2, 3], widespread testing to identify the infected is needed in order to link them to care and decrease the possible consequences of untreated HBV infection, which include approximately 500,000 to 1.2 million deaths yearly from liver cirrhosis and its complications, including primary liver cancer [1]. Limitations related to funding and access to commercially available tools for chronic HBV testing are particularly important in developing countries where the burden of chronic HBV is heaviest. Success of traditional standardofcare (SOC) testing for HBV infection hinges on the existence of a systematic process of following up test results that return several days after testing, notifying patients of results, and arranging for followups to discuss antiviral therapy, a system that requires resources that are limited in developing regions.
The development of rapid pointofcare (POC) tests for HBV has the potential to address many of these limiting factors and establish a more effective medical care model for chronic HBV. In a recent study of patients undergoing HBV screening, the performance characteristics of NanoSignÂ® HBs POC chromatographic immunoassay was compared with standard commercial laboratory HBsAg testing (Quest Diagnostics EIA). The POC tests yielded a sensitivity of 73.7% and a specificity of 97.8% [4]. In a metaanalysis evaluating the accuracy of POC testing, Shivkumar et al. reported POC testing sensitivity of 9398% and specificity of 93â€“99% [5]. Furthermore, the low cost ($0.50) and rapid turnaround (20 min from phlebotomy to test results) of POC tests give them the potential to significantly improve the widespread implementation of HBV screening, especially in resourcelimited regions.
Modeling in HBV analysis and treatment is an active research topic and multiple approaches have been considered recently [6,7,8]. A variety of mathematical modeling strategies have been used to address in particular the costeffectiveness of HBV screening, using predominantly combinations of decision trees and/or Markov chain models [9]. In this paper, we propose timevarying Markov chain models of detailed structure, reflecting disease propagation in individuals to quantify the effects of largescale utilization of POC tests to succeed the SOC screening model.
Methods
In comparing effectiveness of POC and SOC screening strategies for HBV, we made use of two Markov models with identical structure but different transition probabilities. Each of these models has six states capturing the medical progression of individual patients and is formed by aggregating states from a more detailed Markov model describing chronic HBV disease progression of individuals. The two aggregated models were used to simulate consequences of POC or SOC utilization in HBV screening strategies on large populations of individuals. Before they can be iterated numerically, the two Markov models rely on the specification of certain numerical values dealing with the rate of uptake of POC, the rate of infected patients seeking medical care, death rates, and so on. Some of these numbers can be determined (at least within a range) from the medical literature, which we did. Others are more hypothetical or might be the outcome of policy initiatives. The utility of the models is in their low computational cost and attendant capacity for iteration with many possible candidate values and the determination of the sensitivity of the observed behavior to the specific parameter values. Where available, the transition parameters in our models were selected from the literature. The remaining model parameters were estimated and their effects on the overall results analyzed in terms of sensitivity.
For the aggregated models, we considered the six patient states depicted in Fig. 1, where arrows symbolize state transitions admissible in a single time step. At year t, a patient is in state i with probability Ï€ _{ i,t }. Arranging these probabilities into a state rowvector, we have \( {\varPi}_t=\left[\begin{array}{cccc}\hfill {\pi}_{1, t}\hfill & \hfill {\pi}_{2, t}\hfill & \hfill \dots \hfill & \hfill {\pi}_{6, t}\hfill \end{array}\right] \). This vector is then propagated over time via Î _{ tâ€‰+â€‰1}â€‰=â€‰Î _{ t } P _{ t }, where P _{ t } denotes the Markov state transition matrix at time t with element p _{ ij,t } denoting the conditional probability of a patient in state i at year t transitioning to state j at year tâ€‰+â€‰1. Evidently, each row of the state transition matrix sums to one at all times. Special cases are p _{ ij,t }â€‰=â€‰0 for inadmissible transitions and p _{ ij,t }â€‰=â€‰1 for certain transitions. For instance, a transition from being immune to having an undetected HBV infection is inadmissible, whereas a deceased patient is going to remain so.
As illustrated by the connection between States 3 and 4 in Fig. 1, we assume that no patient starting treatment ever abandons this treatment. This assumption is justified in that it alters the transition probabilities p _{36,t } and p _{46,t } of dying from HBV with or without medical treatment by a relatively small degree, which is covered by the sensitivity analysis described below. Model States 3 and 4 in this model are aggregated states from more detailed Markov models described below. This common model structure for both SOC and POC screening policies takes P _{ t } of form
where p _{12,t }â€‰=â€‰p _{12} and p _{34,t }â€‰=â€‰p _{34} are constant transition probabilities chosen depending on the screening strategy at hand. The constant transition probability p _{13} is presumed independent of the screening policies while p _{36} and p _{46} are results of the aggregation procedure outlined below. Timevariation of the state transition matrix P _{ t } is caused solely by the varying propensity for death, p _{5,t } as age advances, which is modeled through linear inter and extrapolation of annual mortality rates for individuals in the USA [10]. However, even though all time variations are induced by variations in p _{5,t }, notice that the changing mortality rates affect the first four transition probabilities on the diagonal of the state transition matrix, being the probabilities to remain in the respective nonabsorbing states. For illustration, whenever p _{5,t } increases at some time, p _{22,t } has to decrease by just as much to ensure that the sum of all transition probabilities from State 2 is one at all times t. That is, we require p _{22,t }â€‰=â€‰1â€‰âˆ’â€‰p _{5,t } for all times t. These timeadjustments have to be made for all rows of the state transition matrix corresponding to nonabsorbing states (i.e., the first four rows of P _{ t }).
Special cases comprise the third and fourth rows of P _{ t }, which capture transitions emerging from the two aggregated states corresponding to HBV disease progression in untreated and treated forms, respectively. We next discuss aggregation of a more complex model capturing the natural history of chronic HBV to estimate the transition probabilities corresponding to State 3 and 4 of the aggregated model in Fig. 1.
To form the aggregation resulting in State 3, untreated infection, we captured the disease progression of HBV without treatment using a different Markov model with its own states and admissible transitions as depicted in Fig. 2. Transition probabilities in this model are based on a literature review and subsequent weighting of the annual probabilities reported in the references using the GRADE criteria [11] for assessing the quality of each study. The resulting transition probabilities are summarized in Table 1, where HCC connotes hepatocellular carcinoma (liver cancer). We refer to the auxiliary model corresponding to Fig. 2 equipped with annual transition probabilities summarized in Table 1 as the disease model, while we refer to the more compact sixstate model described above as the aggregated model.
Absorption in the disease model means transition to the union of the States 4 (Infection under treatment), 5 (death, HBVunrelated) and 6 (death, HBVrelated) of the aggregated model, which can be reached from any state in the disease model.

1)
Initiation of treatment occurs via p _{34,t }â€‰=â€‰p _{34} defined for the aggregated model above.

2)
Death unrelated to HBV in the disease model is assumed timeinvariant with annual probability of 0.1%, corresponding to an individual of age 25â€“34 years [10].

3)
HBVrelated death follows the probabilities listed in Table 1.
The total absorption probability at any state in the disease model is then the sum of the three aforementioned component probabilities. Having the purpose of aggregation in mind, we do not need to distinguish the absorbing states in the disease model any further. The technical tool used to aggregate the disease model into State 3 of the aggregated model is the fundamental matrix Nâ€‰=â€‰(Iâ€‰âˆ’â€‰Q)^{âˆ’â€‰1} of the disease model, where Q is the matrix obtained by extracting all rows and columns of the state transition matrix corresponding to transient states. The fundamental matrix N allows a number of useful deductions about the Markov chain, such as expected numbers of occupancies in transient states until absorption and corresponding variances.
A particularly useful property of the fundamental matrix is that element n _{ ij } equals the expected number of years spent at state j of the disease model when starting in state i [12]. That is, row i of N accumulates the expected numbers of years in each disease state given the process is initiated in state i. Given a sufficiently large number of patients, the normalized version of this rowvector can be interpreted as the average fraction of time spent at each state until absorption, where normalization refers to scaling the vector such that its components have sum one. We can now estimate the probability of death caused by untreated HBV by forming this normalized vector from the first row of the fundamental matrix and using its components to obtain a weighted sum of the probabilities for HBVrelated death in Table 1. This weighted sum is used for p _{36}, concluding the aggregation of the disease model into State 3 of the aggregated model. Assuming probabilities p _{34}â€‰=â€‰15% of initiating medical treatment and 0.1% for HBVunrelated death in the disease model, this procedure yields the estimate p _{36}â€‰=â€‰1.35%. Notice that this probability depends implicitly on the screening policy employed via variation of p _{34}.
State 4, Infection under treatment of the aggregated model can be viewed as an aggregation of the same states used to form the disease model (Fig. 2), although with annual transition probabilities differing from those summarized in Table 1 to reflect effects of treatment. Moreover, the probability of absorption from this disease model under treatment would be decreased by the amount of p _{34}. The effect of these transition probabilities for the disease model under treatment is a value for the HBVrelated mortality rate under treatment in the aggregated model, p _{46}. To obtain this transition probability, we correct the probability for absorption in the disease model used above by p _{34}, but keep using the values in Table 1. To adjust for the favorable effects of medical intervention, we introduce a scaling parameter Î±â€‰âˆˆâ€‰(0,â€‰1) and estimate p _{46}â€‰=â€‰Î±p ^{*}_{46} , where p ^{*}_{46} denotes the probability obtained after aggregation with the values in Table 1. For instance, scaling factor Î±â€‰=â€‰0.25 and fixed annual probability 0.1% for HBVunrelated death result in the estimate p _{46}â€‰=â€‰0.54%. This scaling approach is chosen as we focus on effects of screening policies rather than treatment options. Summarizing the modeling and aggregation procedure, the state transition matrix for the aggregated models takes on the structure
with HBVunrelated mortality rates p _{5,t } from the literature [10] and the transition probabilities in the third and fourth rows depending on the aggregation procedure outlined above. In the following, we use this transition matrix structure for simulation and corresponding sensitivity analyses based on a number of constants in the state transition matrix, namely the transition probabilities p _{12}, p_{13}, p _{34} and the scaling parameter Î±. As mentioned above, the effects of SOC and POC screening strategies are compared using altered transition probabilities p _{12} and p _{34} in the aggregated model, which in terms also affects the state aggregation yielding p _{36} as well as the respective transition probabilities on the diagonal of P _{ t }. Higher utilization of POC screening with subsequent immunization in uninfected cases and initiation of medical treatment in infected cases, respectively, is anticipated to increase both p _{12} and p _{34}, albeit to different degrees. To model these changes, we take p _{12}â€‰â†’â€‰p ^{SOC}_{12} and p _{34}â€‰â†’â€‰p ^{SOC}_{34} in the SOC case. In the POC case, we take p _{12}â€‰â†’â€‰p ^{POC}_{12} â€‰=â€‰Î²p ^{SOC}_{12} and p _{34}â€‰â†’â€‰p ^{POC}_{34} â€‰=â€‰Î³p ^{SOC}_{34} , employing additional scaling parameters Î² and Î³, each greater than one.
The approach taken to analyze POC/SOC utilization effectiveness using the quantified aggregated model is to model a population of a large number of individuals starting with an initial probability distribution over the six states and then to propagate the Markov chain until the collective probability of the death states is nearly one. We presume the population comprises 100,000 initially uninoculated and uninfected 10yearolds and the evolution of the Markov chain over time yields the anticipated proportions of the aging population in each state. The assessment of various performance measures such as mortality rates, years under treatment or life expectancy under SOC and POC screening policies is then tracked via evolution of the probability vector Î _{ t }, which now admits the interpretation as the proportions of the population occupying each disease state, since the population is presumed large.
As mentioned above, we use an annual probability of 0.1% for HBVunrelated death in the disease model. The remaining transition probabilities in the disease model are according to Table 1. The timevarying probabilities for HBVunrelated death in the aggregated model, p _{5,t }, are obtained via inter and extrapolation of mortality data for individuals in the U.S. [10]. The remaining simulation parameters to be chosen are the scaling constants Î±, Î² and Î³ as well as the transition probabilities p_{13}, p ^{SOC}_{12} and p ^{SOC}_{34} . In the following, we use the nominal parameter values p ^{SOC}_{12} â€‰=â€‰0.2%, p ^{SOC}_{34} â€‰=â€‰15%, p _{13}â€‰=â€‰0.15%, Î±â€‰=â€‰0.25, Î²â€‰=â€‰5 and Î³â€‰=â€‰2 unless otherwise specified. Sensitivity analyses around the nominal parameter values specified above are displayed in Table 2, one for the inoculationrelated parameters p ^{SOC}_{12} and Î², one for the treatmentrelated parameters p ^{SOC}_{34} and Î³ and one for the screeningunrelated parameters p _{13} and Î±. Each sensitivity analysis uses nominal values for the remaining parameters.
Results
Simulation results based on our cohort of 100,000 initially uninfected and uninoculated 10yearolds and the nominal simulation parameters specified above are displayed in Fig. 3. As anticipated, utilization of POC screenings reduces the numbers of untreated infections and HBVrelated mortalities at all times. This holds for all setups of the simulation parameters, as long as the scaling parameters are confined to their respective boundaries, that is 0â€‰<â€‰Î±â€‰<â€‰1, Î²â€‰>â€‰1 and Î³â€‰>â€‰1. The number of infections under treatment using POC screening is initially higher but lower in average than with SOC screening. This trend is a result of both the increased inoculation rate p ^{POC}_{12} â€‰>â€‰p ^{SOC}_{12} and the increased rate for initiation of treatment p ^{POC}_{34} â€‰>â€‰p ^{SOC}_{34} . Initially, about the same number of people gets infected under each screening strategy, while more of those infected individuals are linked to medical treatment in the POC case. As the population ages, a larger fraction of the cohort has been inoculated in the POC case, which results in a decreased number of new infections. This in turns leads to a smaller number of patients with untreated infections that can potentially be linked to care, resulting in lower average numbers of patients receiving medical treatment under the POC screening setup. However, the ratio of people linked to care over those infected without treatment is significantly higher in the POC case. The results from our modeling serve to quantify and bound these results, which are a logical consequence of the modelâ€™s structure.
The data displayed in Fig. 3 are based on the particular set of nominal parameter values specified above, with the resultant findings extended to different parameter combinations. Sensitivity analyses around the nominal parameter values specified above are displayed in Table 2 for the following varied parameter combinations: the inoculationrelated parameters p ^{SOC}_{12} and Î²; the treatmentrelated parameters p ^{SOC}_{34} and Î³; and the screeningunrelated parameters p _{13} and Î±. For presentation purposes, each sensitivity analysis presumes nominal values for the remaining simulation parameters. However, the trends summarized in the following paragraphs extend to combined sensitivity analysis. The indicators listed to evaluate the simulations are the total numbers of HBVrelated mortalities and life expectancies under POC and SOC screening policies as well as relative improvements gained by implementing POC screenings. For instance, the nominal parameter values result in improvements of 27.8% in HBVrelated mortality numbers 0.18% in life expectancy, respectively. Sensitivity is interpreted as variation in these two relative measures for the benefit in POC screening utilization. The reason for the seemingly low changes in life expectancy is that only a fraction of the population ever gets infected with HBV, while the change in life expectancy for infected individuals is larger.
In the first sensitivity analysis, only the treatment effectiveness factor Î± and the infection rate p _{13} vary from their nominal values. These are the simulation parameters presumed independent of the screening policies employed. As we can see, each of the tested combinations of these two parameters yields improvements of at least 25.1% in total HBVrelated death numbers and 0.12% in life expectancy in the POC screening case. In general, we observe trends for larger improvements in HBVrelated mortality numbers towards higher treatment effectiveness (i.e., lower value for Î±). The infection rate p _{13} has only minor influence on the improvement in HBVrelated death totals, while increasing the gains in life expectancy at higher infection rates. The comparatively small influence of p _{13} on the relative improvements in mortality numbers is not surprising as p _{13} only changes the proportions of the population ever to become infected, but not the change of course for patients after being infected. The treatment effectiveness factor Î±, however, is strongly linked to potential gains in POC screening by the improved linkage to care and thus affects relative improvements in mortality numbers to a greater extent. The reason for the strong sensitivity of the gains in life expectancy to the infection rate is that if a larger fraction of the population becomes infected, the relative weight of the improvements for this fraction on the entire population grows.
The second sensitivity study focuses on variations of the inoculationrelated simulation parameters p ^{SOC}_{12} and Î². Using the parameter values in Table 2, we gain improvements of at least 21.98% in mortality numbers and at least 0.17% in life expectancies when implementing POC screenings. As expected, both scaling factor Î² and base inoculation rate p ^{SOC}_{12} have significant influence on the two measures of improvement obtainable using POC tests. However, even for low scaling factors Î² and high base inoculation rates, notable benefits of POC test utilization are apparent. In the third sensitivity study, the treatmentrelated simulation parameters p ^{SOC}_{34} and Î³ are varied. Improvements are at least 22.51% in HBVrelated mortality numbers and 0.13% in life expectancy, while both parameters appear to have similar influence on the two measures.
Discussion and conclusions
Chronic HBV is a worldwide problem, with millions of new people infected each year and a large population of chronically infected patients facing health care consequences both short and longterm. However, many chronic HBV patients remain asymptomatic and millions worldwide are unaware of their infections. The importance of early detection via HBV screening of highrisk individuals hinges on the ability to implement effective antiviral therapy to prevent progression of liver disease leading to complications such as cirrhosis and hepatocellular carcinoma. While commercially available serologic immunoassays are widely used for HBV screening, the availability and access to these testing tools for resourcelimited regions or marginalized populations such as the homeless and immigrants are suboptimal. Furthermore, the effort associated with following up on SOC test results, patient callback and counseling can be considerable and create further hurdles for implementing effective screening programs. Recent development of POC tests for HBV holds promise, and previous studies have reported satisfactory sensitivity and specificity of POC testing when compared with SOC testing. However, few studies have used a modeling approach that not only takes into account the performance characteristics of POC testing, but also the natural history of untreated HBV infection to evaluate accurately the added benefit of POC testing over SOC testing. Given the significantly lower cost and more rapid turnaround time associated with POC testing for HBV, the replacement of SOC testing by POC testing has the potential to improve HBV screening programs by promoting greater access and improving linkage to care.
Using Markov modeling based on a comprehensive literature review, our current study demonstrates that POC testing is associated with significantly lower HBVrelated mortality and greater life expectancy when compared with SOC testing. In conclusion, the simulation results under various parameter selections indicate that a significant improvement is obtainable via replacement of SOC screening by new POC tests. The clinical impact of POC testing may be even greater in resourcelimited regions and among marginalized populations where health care access and followup after testing are obstacles to the effective implementation of HBV screening programs. In a future study, additional measures such as morbidity and expected cost of treatment will be analyzed based on additional data regarding cost and effectiveness of medical treatment as well as costs of POC and SOC screening implementation.
Abbreviations
 HBV:

Hepatitis B virus
 POC:

Pointofcare
 SOC:

Standardofcare
References
Lavanchy D. Hepatitis B, virus epidemiology, disease burden, treatment, and current and emerging prevention and control measures. J Viral Hepat. 2004;11:97â€“107.
Gish RG, Gadano AC. Chronic hepatitis B: current epidemiology in the Americas and implications for management. J Viral Hepat. 2006;13:787â€“98.
Kowdley KV, Wang CC, Welch S, Roberts H, Brosgart CL. Prevalence of chronic hepatitis B among foreignborn persons living in the United States by country of origin. Hepatology. 2012;56:422â€“33.
Gish RG, Gutierrez JA, NavarroCazarez N, Giang K, Adler D, Tran B, et al. A simple and inexpensive pointofcare test for hepatitis B surface antigen detection: serological and molecular evaluation. J Viral Hepat. 2014;21:905â€“8.
Shivkumar S, Peeling R, Jafari Y, Joseph L, Pai NP. Rapid pointofcare firstline screening tests for hepatitis B infection: a metaanalysis of diagnostic accuracy (1980â€“2010). Am J Gastroenterol. 2012;107:1306â€“13.
Owolabi KM. Numerical solution of diffusive HBV model in a fractional medium. Springerplus. 2016;5:1643.
Shlomai A, Schwartz RE, Ramanan V, Bhatta A, De Jong YP, Bhatia SN, et al. Modeling host interactions with hepatitis B virus using primary and induced pluripotent stem cellderived hepatocellular systems. Proc Natl Acad Sci U S A. 2014;111:12193â€“8.
Cheng L, Li F, Bility MT, Murphy CM, Su L. Modeling hepatitis B virus infection, immunopathology and therapy in mice. Antiviral Res. 2015;121:1â€“8.
Geue C, Wu O, Xin Y, Heggie R, Hutchinson S, Martin NK, et al. CostEffectiveness of HBV and HCV Screening StrategiesA Systematic Review of Existing Modelling Techniques. PLoS One. 2015;10:e0145022.
Centers for Disease Control and Prevention NCFHS 2013. Compressed Mortality File 1999â€“2010 on CDC WONDER Online Database. http://wonder.cdc.gov/cmficd10.html.
Guyatt GH, Oxman AD, Vist GE, Kunz R, FalckYtter Y, AlonsoCoello P, et al. GRADE: an emerging consensus on rating quality of evidence and strength of recommendations. Br Med J. 2008;336:924â€“6.
Kemeny JG, Snell JL. Finite Markov chains. New York: Springer; 1976.
Eckman MH, Kaiser TE, Sherman KE. The Costeffectiveness of Screening for Chronic Hepatitis B Infection in the United States. Clin Infect Dis. 2011;52:1294â€“306.
Toy M, Veldhuijzen IK, De Man RA, Richardus JH, Schalm SW. Potential impact of longterm nucleoside therapy on the mortality and morbidity of active chronic hepatitis B. Hepatology. 2009;50:743â€“51.
Toy M, Salomon JA, Jiang H, Gui HL, Wang H, Wang JS, et al. Population health impact and costeffectiveness of monitoring inactive chronic hepatitis B and treating eligible patients in Shanghai, China. Hepatology. 2014;60:46â€“55.
Veldhuijzen IK, Toy M, Hahne SJM, De Wit GA, Schalm SW, De Man RA, et al. Screening and early treatment of migrants for chronic hepatitis B virus infection is costeffective. Gastroenterology. 2010;138:522â€“30.
Wong WWL, Woo G, Heathcote EJ, Krahn M. Cost effectiveness of screening immigrants for hepatitis B. Liver Int. 2011;31:1179â€“90.
Lacey LF, Gane E. The costeffectiveness of longterm antiviral therapy in the management of HBeAgpositive and HBeAgnegative chronic hepatitis B in Singapore. J Viral Hepat. 2007;14:751â€“66.
Chen JD, Yang HI, Iloeje UH, You SL, Lu SN, Wang LY, et al. Carriers of Inactive Hepatitis B Virus Are Still at Risk for Hepatocellular Carcinoma and LiverRelated Death. Gastroenterology. 2010;138:1738â€“47.
Lok A, Lai CL, Wu PC, Leung E. Spontaneous hepatitis B e antigen to antibody seroconversion and reversion in Chinese patients with chronic hepatitis B virus infection. Gastroenterology. 1987;92:1839â€“43.
Yuen MF, Wong DKH, Fung J, Ip P, But D, Hung I, et al. HBsAg seroclearance in chronic hepatitis B in Asian patients: Replicative level and risk of hepatocellular carcinoma. Gastroenterology. 2008;135:1192â€“9.
Arase Y, Ikeda K, Suzuki F, Suzuki Y, Saitoh S, Kobayashi M, et al. Longterm outcome after hepatitis B surface antigen seroclearance in patients with chronic hepatitis B. Am J Med. 2006;119:71e9â€“16.
Chen YC, Sheen IS, Chu CM, Liaw YF. Prognosis following spontaneous HBsAg seroclearance in chronic hepatitis B patients with or without concurrent infection. Gastroenterology. 2002;123:1084â€“9.
Liu J, Yang HI, Lee MH, Lu SN, Jen CL, BatrlaUtermann R, et al. Spontaneous seroclearance of hepatitis B seromarkers and subsequent risk of hepatocellular carcinoma. Gut. 2014;63:1648â€“57.
Aggarwal R, Ghoshal UC, Naik SR. Assessment of costeffectiveness of universal hepatitis B immunization in a lowincome country with intermediate endemicity using a Markov model. J Hepatol. 2003;38:215â€“22.
Hutton DW, So SK, Brandeau ML. Costeffectiveness of nationwide hepatitis B catchup vaccination among children and adolescents in China. Hepatology. 2010;51:405â€“14.
Hutton DW, Tan D, So SK, Brandeau ML. Costeffectiveness of screening and vaccinating Asian and Pacific Islander adults for hepatitis B. Ann Intern Med. 2007;147:460â€“9.
Acknowledgements
Writing and editing support was provided by independent medical editor Lark Lands, Ph.D., and funded by Robert Gish, M.D. The authors would like to express their gratitude for her invaluable assistance in preparing the manuscript for publication.
Funding
Funds were provided through a special projects fund at the University of California, San Diego.
Availability of data and materials
All data not explicitly included in the article is simulation data, which the authors will make available upon request.
Authorsâ€™ contributions
MAS contributed substantially to the conception and design of the study, the interpretation of data, both the drafting and the critical revision of the manuscript for important intellectual content, and the final approval of the version to be published; he was the lead researcher who was responsible for modeling procedure, simulation, and sensitivity analysis; he agrees to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the article are appropriately investigated and resolved. KDJ contributed substantially to the conception and design of the study, the analysis and interpretation of data, both the drafting and the critical revision of the manuscript for important intellectual content, and the final approval of the version to be published; he contributed substantially to the comprehensive literature reviews used to quantify the models, and contributed to model procedure; he agrees to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the article are appropriately investigated and resolved. JMF contributed substantially to the conception and design of the study, both the drafting and the critical revision of the manuscript for important intellectual content, and the final approval of the version to be published; he contributed substantially to modeling procedure, simulation, and sensitivity analysis; he agrees to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the article are appropriately investigated and resolved. RJW contributed substantially to the conception and design of the study; analysis and interpretation of data, both the drafting and the critical revision of the manuscript for important intellectual content, and the final approval of the version to be published; he contributed substantially to the comprehensive literature reviews used to quantify the models and to all the medical discussion elements in the manuscript; he agrees to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the article are appropriately investigated and resolved. RRB contributed substantially to the conception and design of the study, the analysis and interpretation of data, both the drafting and the critical revision of the manuscript for important intellectual content, and the final approval of the version to be published; he contributed substantially to modeling procedure, simulation, and sensitivity analysis; he agrees to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the article are appropriately investigated and resolved. RGG contributed substantially to the concept and design of the study, both the drafting and the critical revision of the manuscript for important intellectual content, and the final approval of the version to be published; he contributed substantially to the manuscriptâ€™s discussion of hepatitis B testing, including standardofcare and pointofcare tests, hepatitis B epidemiology and disease progression, and approaches to hepatitis B screening and linkage to care; he agrees to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the article are appropriately investigated and resolved.
Competing interests
The authors list as possible competing interests the following: Robert G. Gish has had Grants/Research Support from Gilead Sciences, and Merck & Co.; Dr. Gish has performed as Consultant and/or Advisor to Akshaya Pharmaceuticals, Arbutus Biopharma Corporation, Arrowhead Research Corporation, BristolMyers Squibb, ContraVir Pharmaceuticals, Enyo Pharma, Gilead Sciences, HumAbs BioMed, Ionis Pharmaceuticals, Merck & Co., Nanogen Biopharmaceutical, and Novira Therapeutics; Dr. Gish has current activity with the scientific or clinical advisory boards of Arrowhead Research Corporation, Merck & Co., ContraVir Pharmaceuticals, Gilead Sciences, Isis Pharmaceuticals, Enyo Pharma, HumAbs BioMed, and Nanogen Biopharmaceutical; Dr. Gish is a member of the Speakers Bureau for BristolMyers Squibb, Gilead Sciences, and Merck & Co.; Dr. Gish has stock options with Arrowhead Research Corporation. All other authors have no competing interests.
Consent for publication
This study does not contain any individual personâ€™s data in any form.
Ethics approval and consent to participate
This study does not involve human participants, human data or human tissue.
Publisherâ€™s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Sehr, M.A., Joshi, K.D., Fontanesi, J.M. et al. Markov modeling in hepatitis B screening and linkage to care. Theor Biol Med Model 14, 11 (2017). https://doi.org/10.1186/s1297601700576
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1297601700576
Keywords
 Hepatitis B virus
 Screening
 Markov modeling
 Pointofcare
 Standardofcare
 Testing