Loading [MathJax]/extensions/TeX/boldsymbol.js
Research article Special Issues

DeepMRMP: A new predictor for multiple types of RNA modification sites using deep learning

  • RNA modification plays an indispensable role in the regulation of organisms. RNA modification site prediction offers an insight into diverse cellular processing. Regarding different types of RNA modification site prediction, it is difficult to tell the most relevant feature combinations from a variant of RNA properties. Thereby, the performance of traditional machine learning based predictors relied on the skill of feature engineering. As a data-driven approach, deep learning can detect optimal feature patterns to represent input data. In this study, we developed a predictor for multiple types of RNA modifications method called DeepMRMP (Multiple Types RNA Modification Sites Predictor), which is based on the bidirectional Gated Recurrent Unit (BGRU) and transfer learning. DeepMRMP makes full use of multiple RNA site modification data and correlation among them to build predictor for different types of RNA modification sites. Through 10-fold cross-validation of the RNA sequences of H. sapiens, M. musculus and S. cerevisiae, DeepMRMP acted as a reliable computational tool for identifying N1-methyladenosine (m1A), pseudouridine (Ψ), 5-methylcytosine (m5C) modification sites.

    Citation: Pingping Sun, Yongbing Chen, Bo Liu, Yanxin Gao, Ye Han, Fei He, Jinchao Ji. DeepMRMP: A new predictor for multiple types of RNA modification sites using deep learning[J]. Mathematical Biosciences and Engineering, 2019, 16(6): 6231-6241. doi: 10.3934/mbe.2019310

    Related Papers:

    [1] Abdul Saboor Karzad, Moussa Leblouba, Zaid A. Al-Sadoon, Mohamed Maalej, Salah Altoubat . Modeling the flexural strength of steel fibre reinforced concrete. AIMS Materials Science, 2023, 10(1): 86-111. doi: 10.3934/matersci.2023006
    [2] Albert Uchenna Ude, Che Husna Azhari . Lateral crashworthiness response of bombyx mori fibre/glass–fibre/epoxy hybrid composite cylindrical tubes-experimental. AIMS Materials Science, 2019, 6(6): 1227-1239. doi: 10.3934/matersci.2019.6.1227
    [3] Zulzamri Salleh, Md Mainul Islam, Jayantha Ananda Epaarachchi, Haibin Su . Mechanical properties of sandwich composite made of syntactic foam core and GFRP skins. AIMS Materials Science, 2016, 3(4): 1704-1727. doi: 10.3934/matersci.2016.4.1704
    [4] Sanna Virtanen, Lisa Wikström, Kirsi Immonen, Upi Anttila, Elias Retulainen . Cellulose kraft pulp reinforced polylactic acid (PLA) composites: effect of fibre moisture content. AIMS Materials Science, 2016, 3(3): 756-769. doi: 10.3934/matersci.2016.3.756
    [5] Md Shamsuddoha, Md Mainul Islam, Thiru Aravinthan, Allan Manalo, Luke P. Djukic . Effect of hygrothermal conditioning on the mechanical and thermal properties of epoxy grouts for offshore pipeline rehabilitation. AIMS Materials Science, 2016, 3(3): 832-850. doi: 10.3934/matersci.2016.3.832
    [6] Albert Uchenna Ude, Che Husna Azhari . Crashworthiness response of natural silk-fibre glass hybrid reinforced epoxy cylindrical composite tubes under quasi-static load. AIMS Materials Science, 2019, 6(5): 852-863. doi: 10.3934/matersci.2019.5.852
    [7] M. Rajanish, N. V. Nanjundaradhya, Ramesh S. Sharma, H. K. Shivananda, Alok Hegde . Directional Interlaminar Shear Strength (ILSS) of nano-modified epoxy/unidirectional glass fibre composite. AIMS Materials Science, 2018, 5(4): 603-613. doi: 10.3934/matersci.2018.4.603
    [8] Timothy K. Mulenga, Albert U. Ude, Chinnasamy Vivekanandhan . Concise review on the mechanical characteristics of hybrid natural fibres with filler content. AIMS Materials Science, 2020, 7(5): 650-664. doi: 10.3934/matersci.2020.5.650
    [9] M. A. Mendizabal, Maitane Garcia, Luis Palenzuela, Enrique Hernández . Obtaining preforms by additive fused deposition modelling (FDM) extrusion technology for the manufacture of high-performance composites. AIMS Materials Science, 2022, 9(3): 481-497. doi: 10.3934/matersci.2022028
    [10] Sohel Rana, Subramani P, Raul Fangueiro, Antonio Gomes Correia . A review on smart self-sensing composite materials for civil engineering applications. AIMS Materials Science, 2016, 3(2): 357-379. doi: 10.3934/matersci.2016.2.357
  • RNA modification plays an indispensable role in the regulation of organisms. RNA modification site prediction offers an insight into diverse cellular processing. Regarding different types of RNA modification site prediction, it is difficult to tell the most relevant feature combinations from a variant of RNA properties. Thereby, the performance of traditional machine learning based predictors relied on the skill of feature engineering. As a data-driven approach, deep learning can detect optimal feature patterns to represent input data. In this study, we developed a predictor for multiple types of RNA modifications method called DeepMRMP (Multiple Types RNA Modification Sites Predictor), which is based on the bidirectional Gated Recurrent Unit (BGRU) and transfer learning. DeepMRMP makes full use of multiple RNA site modification data and correlation among them to build predictor for different types of RNA modification sites. Through 10-fold cross-validation of the RNA sequences of H. sapiens, M. musculus and S. cerevisiae, DeepMRMP acted as a reliable computational tool for identifying N1-methyladenosine (m1A), pseudouridine (Ψ), 5-methylcytosine (m5C) modification sites.


    In recent years, large-scale epidemic outbreaks involving asymptomatic patients have attracted widespread attention, especially COVID-19 since the end of 2019, which poses a huge threat to global human health, life safety and economic development. Due to the significant hazards caused by infectious diseases, countries have to adopt appropriate prevention and control strategies to contain virus transmission once an outbreak occurs. For example, some countries have implemented measures such as closing schools and businesses to restrict gatherings, and locking down cities to limit mobility (UK Government. Coronavirus (COVID-19): guidance and support. https://www.gov.uk/coronavirus. Accessed March 17, 2023). Evidently, rapid control of the outbreak became an important issue.

    Asymptomatic carriers of infectious diseases pose significant challenges for epidemic control, especially the COVID-19 pandemic. They are individuals who tested positive for the disease but do not have any clinical symptoms that they can perceive [1]. Based on this meaning of asymptomatic patients, we refer to pre-symptomatic patients and those who are asymptomatic throughout the duration of their illness collectively as asymptomatic patients in this work. While asymptomatic carriers may not be screened or diagnosed, they can still carry the virus and become a potential source of infection. Those patients may have longer exposure times and contact with more people than symptomatic patients, as they are unaware of their infection. Therefore, asymptomatic carriers are difficult to detect and isolate in a timely manner so that they directly cause the high transmission rate due to their hidden transmission characteristics. Neglecting the infection caused by asymptomatic carriers has led to continued transmission and repeated outbreaks in some countries and regions [2]. Some studies have shown that the existence of asymptomatic carriers is one of the important and key reasons for the difficulty in controlling the COVID-19 pandemic [3,4,5]. Therefore, it is an important topic to study the mechanisms of transmission with asymptomatic carriers and to find the key factors during the long-term transmission process of infectious diseases, so as to predict epidemic trends and develop effective prevention and control strategies.

    At present, research on the mechanisms of transmission with asymptomatic carriers mainly focuses on two classes: qualitative analysis and quantitative analysis. We review recent studies on infectious diseases with asymptomatic patients through taking COVID-19 as an example. For the first class, several scholars reviewed the characteristics of asymptomatic COVID-19 patients, mainly concerning contagiousness [6,7,8], infectious characteristics [8,9,10] and proportion of asymptomatic patients among all COVID-19 patients [8]. In terms of contagiousness, Nelson [11] and Yu et al. [12] found that asymptomatic patients are also infectious and there are cases of interpersonal infection. For the infectious characteristics study, Han et al. [13] and Hu et al. [14] found that asymptomatic patients had longer durations of virus presence, which posed a greater risk of infection. He et al. [15] found that asymptomatic patients had lower infection rates than symptomatic patients, and that patients infected by asymptomatic patients were more likely to be asymptomatic patients. Chen Yi et al. [16] found that the viral infection rates were about 6.30% and 4.11% for contact with symptomatic and asymptomatic patients, respectively. Zhou et al. [17] found that some asymptomatic patients transformed into symptomatic patients, and some patients remained asymptomatic patients. Kimball et al. [18] studied the proportion of asymptomatic patients who are transformed to symptomatic patients. Lauer et al. [19] observed that asymptomatic patients transformed to symptomatic patients after approximately 5 days of latency. For the problem of the proportion of asymptomatic patients, several authors investigated the proportion of asymptomatic patients in various sizes of the population of COVID-19 patients. Wu et al. [20] reviewed 72,314 patients with COVID-19 and found that the proportion of asymptomatic patients was about 1%. Using an example of 565 COVID-19 patients, Nishiura et al. [21] employed a Bayes theorem to estimate that the proportion of asymptomatic patients was 30.8%. Day [22] illustrated that the proportion of asymptomatic COVID-19 patients may reach four-fifths based on a published sample of 166 COVID-19 patients in China.

    For the second class of study, the quantitative study, considering the transmission process of COVID-19 with asymptomatic patients, most research used dynamical system theory, simulation or the complex network approach. Using dynamical systems theory, Sang Maosheng et al. [23] and Agrawal et al. [24] developed a modified susceptible–infectious–recovered (SIR) model to fit the transmission trend of COVID-19 considering asymptomatic patients. Masaki et al. [25] constructed a new mathematical model called SIIR model and found that the conditions for population immunity became more restrictive when considering the limited duration of immunity in asymptomatic patients. Some scholars [26,27,28,29,30] explored the effect of infection rates in asymptomatic patients on the transmission of COVID-19 or the effectiveness of control strategies by the basic reproduction number R0. Riyapan et al. [31] and Iaa et al. [32] proposed a modified SIR model considering asymptomatic patients and found that COVID-19 will die out when the basic reproduction number is less than 1 and will continue to exist when the basic reproduction number is bigger than 1. Simulation methods have also been used to study the COVID-19 transmission process, supporting the inference of potentially infected individuals based on the temporal and location information of their activities. Zhan et al. [33] estimated the basic reproduction number for the COVID-19 transmission process on the Diamond Princess by the simulation method. Yu et al. [34] estimated the number of asymptomatic patients and their transmission rate based on a machine learning simulation method. By means of the complex network approach, some scholars studied the COVID-19 transmission process. Chen et al. [35] set up a novel COVID-19 transmission model through considering asymptomatic patients based on the complex network approach to make accurate forecasts of epidemic trends in the presence of incomplete information. Stella et al. [36] developed a complex network model to study the impact of asymptomatic patients and control measures on transmission while capturing the interactions of symptomatic and asymptomatic COVID-19 patients separately.

    To study the transmission process of COVID-19 with asymptomatic patients, the construction of SIR models by using dynamic systems theory is necessary to parameterize the probability of exposure, infection, recovery, death, etc. for representing the change in the health status of the infected person. It is inferior for individual spatial movement and interaction processes, i.e., the heterogeneity of contact exposure among individuals. Therefore, the SIR-based models are only suitable for studying the infectious disease transmission process in a short time [37,38]. The simulation method is closer to the real situation and can support the analysis of temporal transmission mechanisms of infectious diseases and the design of appropriate interventions. However, the acquisition of a large amount of real individual trajectory data as a support is a challenge [39]. Using the complex network to simulate the contact among transmission subjects can reflect the variability of individual contacts in the social network structure. Unfortunately, it cannot represent the spatio-temporal spread of the disease well [38]. To quantitatively study the long-term transmission of infectious diseases with asymptomatic patients, it is necessary to develop a new mathematical model that depicts the characteristics of people with cross-regional movement in real life, that is, an infectious disease transmission process with spatial movement and interaction, and there is no limit to the number of susceptible persons. To our knowledge, there is still a lack of effective methods to discuss the stability of the long-term transmission and conduct system, and to provide performance analysis for the infectious disease transmission process with asymptomatic patients.

    The motivation of this work is to study a long-term transmission process with realistic characteristics that can be used to analyze the mechanism of asymptomatic carriers. To this end, we develop a new two-dimensional Markov process of infectious disease by means of the level-dependent continuous-time quasi-birth-and-death (QBD) process with infinite-state space. This Markov process can characterize interindividual interactions on infinite scales for patient groups containing both symptomatic and asymptomatic patients [40,41]. On the one hand, this Markov process, as a stochastic process, can better reflect stochasticity and suitability for stability analysis compared to deterministic models. On the other hand, the two-dimensional structure can reflect the complex interindividual interactions among symptomatic and asymptomatic patients during spatio-temporal transmission. Therefore, it is suitable for analyzing the long-term transmission process of infectious diseases with asymptomatic patients. Based on the actual situation of infectious diseases and some studies in the literature, the level on the infinite state space is used to represent symptomatic patients and the phase on the finite state space is used to represent asymptomatic patients. Such a level-dependent Markov process can reflect the different stages of the long-term transmission process of infectious diseases. Furthermore, to deal with the level-dependent QBD processes of infectious disease with infinite-state space, we use the RG-factorizations to find a feasible solution. Additionally, we develop two effective RG-factorization algorithms to obtain the feasible solution that can numerically analyze performance measures of this system of infectious disease.

    The main contributions of this study are summarized as follows:

    1) We develop a new two-dimensional continuous-time QBD process with an infinite-state space in the context that infectious diseases spread over a long time. This Markov process represents the status and transfer rate in different types of patients in the long-term and also expresses the spatial movement and interaction of individuals.

    2) We provide a method based on the mean drift technique to obtain a simple sufficient condition for the stability of the proposed Markov process. The stationary probability vector is obtained by using RG-factorization. Additionally, we develop two algorithms to obtain and discuss some key performance measures of this system.

    3) We conduct some numerical analysis to verify the feasibility of the proposed model by using the real data in the context of the COVID-19 epidemic and discuss how some key parameters affect the key performance measures of this system.

    The remainder of this paper is organized as follows. In Section 2, we describe the infectious disease transmission process with asymptomatic patients and introduce the basic parameters. Section 3 establishes a Markov process of infectious disease considering asymptomatic patients. Section 4 obtains the stationary probability vector of the proposed model by using RG-factorizations and computes some key performance measures of the system by means of the stationary probability vector. Section 5 designs two algorithms to get the stationary probability vector and analyzes performance measures of the system. Section 6 provides numerical cases to verify the feasibility of the model and analyzes the impact of some key parameters on performance measures and infectious disease transmission. Section 7 makes conclusions and presents future research.

    In this section, we describe the infectious disease transmission process with asymptomatic patients, and give the relevant symbols and parameters used in the modeling.

    In the infectious disease transmission process with asymptomatic patients, there are four types of patients: symptomatic, asymptomatic, cured and dead patients as shown in Figure 1. Symptomatic patients refer to patients with confirmed infectious diseases who have relevant symptoms. According to the meaning of asymptomatic patients, they include pre-symptomatic patients and patients who always had no relevant symptoms during their illness. Both symptomatic and asymptomatic patients are collectively referred to as infected patients. Asymptomatic patients may be transformed into symptomatic or cured patients. The end of disease progression in symptomatic and asymptomatic patients may be cured or death.

    Figure 1.  Four types of patients in the infectious disease transmission process with asymptomatic patients. The icon defined by the patient type in this figure is represented consistently throughout the next figures.

    In the next, we give the detailed transmission process and assumptions, which are used to set up a Markov process of infectious disease.

    1) The first appearance process of patients: We assume that the interval length of the number of symptomatic patients from 0 to 1 in an observation region is exponentially distributed with mean 1/τA, and the interval length of asymptomatic patients from 0 to 1 is exponentially distributed with mean 1/τB.

    2) The infection process of patients: We assume that the infection process of patients is a Poisson process, which is consistent with Giri [42] and Alawiyah [43]. Both symptomatic and asymptomatic patients have the ability to transmit the disease to others, resulting in the generation of new cases that include both symptomatic and asymptomatic individuals. In the condition that a symptomatic case exists, the infection rates of generating a symptomatic and an asymptomatic patient are denoted by λAA and λAB, respectively. Similarly, in the condition that an asymptomatic case exists, the infection rates of generating a symptomatic and an asymptomatic patient are denoted by λBA and λBB, respectively. The infection process of an infectious disease is shown in Figure 2.

    Figure 2.  The infection process of a Markov process for an infectious disease. In Figure 2(a), the left side of the arrow indicates the situation before infection, the right side indicates the possible situation after infection, and the number on the arrow indicates the infection rate. Figure 2(b) is similar to Figure 2(a).

    3) The disappearance process of patients: The findings of Davies [44,45] showed that the infection period of COVID-19 fits an exponential distribution. In line with the literature, we assume that the interval length of symptomatic patients from being identified to being cured or death is exponentially distributed with mean 1/μA, and the interval length of asymptomatic patients from being identified to self-healing is exponentially distributed with mean 1/μB. The disappearance process of infected patients is shown in Figure 3. It is worth noting that the disappearance of symptomatic patients can be followed by cured patients or dead patients. However, the disappearance of asymptomatic patients can only be followed by cured patients. Realistically speaking, asymptomatic patients generally go through a symptomatic phase before they die.

    Figure 3.  The disappearance process of infected patients. In Figure 3(a), the left side of the arrow indicates the situation before disappearance, the right side indicates the situation after disappearance, and the number on the arrow indicates the disappearance rate. Figure 3(b) is similar to Figure 3(a).

    4) The transformation process of patients: Based on the setting that pre-symptomatic patients are also asymptomatic, some studies have found that asymptomatic patients can transform into symptomatic patients, but symptomatic patients cannot transform into asymptomatic patients [17,46,47]. Pasetto [48] shows that the mean duration that infected patients stay in the incubation period is exponentially distributed. Based on it, we assume that the interval length of transformation from asymptomatic patients to symptomatic patients is exponentially distributed with mean 1/γ.

    5) Limited number of asymptomatic patients: According to the literature [20,21,22], the asymptomatic patients gradually transform to symptomatic patients as the epidemic evolves. Therefore, we reasonably assume that the number of symptomatic and asymptomatic patients takes values in {0,1,2,} and {0,1,2,,M}, respectively. Note that this assumption can simplify the model construction, proof and analysis.

    6) Independence: We assume that all the random variables described above are independent of each other.

    According to He et al. [15] and Chen Yi et al. [16], the infection rates and the interval length differ in different types of patients. Therefore, four different types of infection rates are assumed in Assumption (2), and the two different types of disappearing rates are assumed in Assumption (3). Although in this paper we describe the infectious disease transmission process by Poisson process and exponential distribution, our mathematical modeling can be easily extended, for example, from Poisson process to Markovian arrival processes (MAPS) in Assumption (2) and from exponential distribution to phase type (PH) distribution in Assumption (3).

    Using the key parameters assumed above, in the next section we can develop a new Markov process of infectious disease with infinite-state space.

    In this section, we model the infectious disease transmission process as a Markov process. Specifically, this Markov process is a continuous-time level-dependent QBD process with an infinite state space.

    Let NA(t) and NB(t) be the number of symptomatic patients and asymptomatic patients at time t, respectively. Based on the assumptions, NA(t){0,1,2,} and NB(t){0,1,2,,M} {(NA(t),NB(t)):t0}. Therefore, it is known that {(NA(t),NB(t)):t0} is a continuous-time Markov process with a state space of

    Ω=k=0Level k, (1)

    where the Level k denotes the set of states in level k, for details:

    Level k={(k,0),(k,1),(k,2),,(k,M)},k0. (2)

    Based on the problem description and using the states in all Levels, we can give the state transition relations for this Markov process {(NA(t),NB(t)):t0}, as shown in Figure 4.

    Figure 4.  The state transition relations for the Markov process. In Figure 4, the infection process, the disappearance process and the transformation process are shown by red, green and blue arrows, respectively. If the system is in the state where the arrow starts at a certain moment, it will later transit to one of the states where the arrows end. The symbols on the arrows indicate the transition rates.

    Further, according to the state transition relation given in Figure 4, the infinitesimal generator of this Markov process can be obtained. This infinitesimal generator contains a great deal of the information about our proposed Markov process, and the model analysis and solution in this paper are carried out based on this infinitesimal generator [49]. The infinitesimal generator can be written as

    Q=(A0,0A0,1A1,0A1,1A1,2A2,1A2,2A2,3An1,nAn,nAn,n+1), (3)

    where Ai,j are all square arrays of order M+1. To easily give the details of Ai,j, let an,i=n(μA+λAA+λAB)+i(γ+μB+λBA+λBB), bn,i=nλAB+iλBB and dn,i=nλAA+iλBA. Using an,i, bn,i and dn,i, we can give the details of Ai,j as follows:

    A0,0=((τA+τB)τBμBa0,1b0,12μBa0,2b0,2(M1)μBa0,M1b0,M1MμBa0,M), (4)
    A0,1=(τAγd0,12γd0,2(M1)γd0,M1Mγd0,M). (5)

    When k1,

    Ak,k1=(kμAkμAkμAkμAkμA), (6)
    Ak,k=(ak,0bk,0μBak,1bk,12μBak,2bk,2(M1)μBak,M1bk,M1MμBak,M), (7)
    Ak,1=(dk,0γdk,12γdk,2(M1)γdk,M1Mγdk,M). (8)

    Based on this state transition relation and infinitesimal generator, we can know that the Markov process of infectious disease is a level-dependent QBD process with an infinite state space. Specifically, a state within a given level can only transition to other states in the same level or in the adjacent levels. Moreover, the process exhibits an infinite state space that arises from the infinite number of levels. Finally, the process is level-dependent, as the state transition matrices differ across levels.

    For the level-dependent QBD process we propose, it is easy to prove its stability. When μA>λAA, this QBD process of infectious disease (Q) is irreducible and positive recurrent. Thus, the QBD process of infectious disease is stable. We give details of the proof process in the appendix.

    In this section, we present how to use the RG-factorization method to obtain the stationary probability vector of the QBD process proposed in this paper and provide some important performance measures of the model based on the stationary probability vector. RG-factorization is an efficient way of solving infinite dimensional (or large) systems of linear equations. More details about RG-factorization can be found in the literature [50,51].

    Let Pi,j(t) be the probability that the process is at state (i,j) at time t. We have

    Pi,j(t)=P{NA(t)=i,NB(t)=j}. (9)

    When the QBD process of infectious disease is stable, let πi,j denote the stationary probability that the process is at state (i,j). πi,j can be obtained by

    πi,j=limtPi,j(t), (10)

    After that, we denote the stationary probability vector that the process is at level i by {{\boldsymbol{\pi }}_i} . For any integer i = 0, 1, 2, …, we have

    {{\boldsymbol{\pi }}_i} = \left( {{\pi _{i,0}},{\pi _{i,1}},{\pi _{i,2}}, \ldots ,{\pi _{i,M}}} \right), (11)

    and the stationary probability vector {\boldsymbol{\pi }} of this process is

    {\boldsymbol{\pi }} = \left( {{{\boldsymbol{\pi }}_0},{{\boldsymbol{\pi }}_1},{{\boldsymbol{\pi }}_2}, \ldots } \right). (12)

    Using the RG-factorizations method [49], we can calculate the stationary probability vector of this QBD process. Therefore, we give the U - , R - and G - measures of the QBD process respectively.

    {U_k} = {A_{k,k}} + {A_{k,k + 1}}( - U_{k + 1}^{ - 1}){A_{k + 1,k}},k \geqslant 0, (13)
    {R_k} = {A_{k,k + 1}}( - U_{k + 1}^{ - 1}),k \geqslant 0, (14)
    {G_k} = ( - U_{k + 1}^{ - 1}){A_{k,k - 1}},k \geqslant 1. (15)

    Note that the matrix sequence \left\{ {{R_k}, k \geqslant 0} \right\} is the smallest non-negative solution to the system of nonlinear matrix Eq (16):

    {A_{k,k + 1}} + {R_k}{A_{k + 1,k + 1}} + {R_k}{R_{k + 1}}{A_{k + 2,k + 1}} = 0,k \geqslant 0. (16)

    Similarly, the matrix sequence \left\{ {{G_k}, k \geqslant 1} \right\} is the smallest non-negative solution to the system of nonlinear matrix Eq (17):

    {A_{k,k + 1}}{G_{k + 1}}{G_k} + {A_{k,k}}{G_k} + {A_{k,k - 1}} = 0,k \geqslant 1. (17)

    Once the matrix sequence \left\{ {{R_k}, k \geqslant 0} \right\} or \left\{ {{G_k}, k \geqslant 1} \right\} is given, for k \geqslant 0 we have

    \begin{gathered} {U_k} = {A_{k,k}} + {R_k}{A_{k + 1,k}} \hfill \\ {\text{ }} = {A_{k,k}} + {A_{k,k + 1}}{G_{k + 1}}. \hfill \\ \end{gathered} (18)

    By following the method described in Chapter 1 of Li [49], we can factorize the infinitesimal generator of the level-dependent QBD process as

    Q = (I - {R_U}){U_D}(I - {G_L}), (19)

    where,

    {U_D} = {\text{diag}}({U_o},{U_1},{U_2} \ldots ), (20)
    {R_U} = \left( {\begin{array}{*{20}{c}} 0&{{R_0}}&{}&{} \\ {}&0&{{R_1}}&{} \\ {}&{}&0& \ddots \\ {}&{}&{}& \ddots \end{array}} \right), (21)
    {G_L} = \left( {\begin{array}{*{20}{c}} 0&{}&{}&{} \\ {{G_1}}&0&{}&{} \\ {}&{{G_2}}&0&{} \\ {}&{}& \ddots & \ddots \end{array}} \right). (22)

    The {\text{diag}}({U_o}, {U_1}, {U_2} \ldots) denotes the diagonal matrix composed of {U_o} , {U_1} , {U_2} , ….

    After that, using this UL-type RG-Factorization and R-measure sequence \left\{ {{R_k}, k \geqslant 0} \right\} , we have

    {{\boldsymbol{\tilde \pi }}_k} = {{\boldsymbol{\tilde \pi }}_0}{R_0}{R_1} \cdots {R_{k - 1}},k \geqslant 1, (23)
    {{\boldsymbol{\pi }}_k} = c{{\boldsymbol{\tilde \pi }}_k},k \geqslant 0, (24)

    where c is the normalization constant such that the sum of the stationary probabilities of all states is 1. And the positive constant c is uniquely given by

    c = \frac{1}{{{{{\boldsymbol{\tilde \pi }}}_0}{\bf{e}} + \sum\limits_{k = 1}^\infty {{{{\boldsymbol{\tilde \pi }}}_0}{R_0}{R_1} \cdots {R_{k - 1}}} {\bf{e}}}}. (25)

    In Eq (23), {{\boldsymbol{\tilde \pi }}_k} represents the stationary probabilities vector of level k before normalization. {{\boldsymbol{\tilde \pi }}_0} can be uniquely determined by the system of linear Eq (40).

    \left\{ \begin{gathered} {{{\boldsymbol{\tilde \pi }}}_0}({A_{0,0}} + {R_0}{A_{1,0}}) = 0 \hfill \\ {{{\boldsymbol{\tilde \pi }}}_0}{\bf{e}} = 1 \hfill \\ \end{gathered} \right.. (26)

    In the remainder of this section, we provide performance analysis of the infectious disease transmission process with asymptomatic patients by means of the stationary probability vector of the QBD process of infectious diseases. Therefore, we can use this stationary probability vector {\boldsymbol{\pi }} to provide some useful performance measures as follows:

    i. Average number of symptomatic patients

    \begin{gathered} E[{N_A}] = \sum\limits_{k = 0}^\infty {k{{\boldsymbol{\pi }}_k}{\bf{e}}} \hfill \\ {\text{ }} = \sum\limits_{k = 1}^\infty {kc{{{\boldsymbol{\tilde \pi }}}_0}{R_0}{R_1} \cdots {R_{k - 1}}{\bf{e}}} , \hfill \\ \end{gathered} (27)

    where, {N_A} = \mathop {\lim }\limits_{t \to \infty } {N_A}(t) , a.s.

    ii. Average number of asymptomatic patients

    \begin{gathered} E[{N_B}] = \sum\limits_{k = 0}^\infty {{{\boldsymbol{\pi }}_k}{{\bf{f}}_B}} \hfill \\ {\text{ }} = c{{{\boldsymbol{\tilde \pi }}}_0}{{\bf{f}}_B} + \sum\limits_{k = 1}^\infty {c{{{\boldsymbol{\tilde \pi }}}_0}{R_0}{R_1} \cdots {R_{k - 1}}{{\bf{f}}_B}} , \hfill \\ \end{gathered} (28)

    where, {N_B} = \mathop {\lim }\limits_{t \to \infty } {N_B}(t) , a.s.; {{\bf{f}}_B} = {\left({0, 1, 2, \ldots, M} \right)^T} .

    iii. Average infection rate from symptomatic patients to symptomatic patients

    \begin{gathered} {r_{A \to A}} = {\lambda _{AA}}({\pi _{1,0}} + 2{\pi _{2,0}} + 3{\pi _{3,0}} + \cdots + {\pi _{1,1}} + 2{\pi _{2,1}} + 3{\pi _{3,1}} + \cdots + \hfill \\ {\text{ }}{\pi _{1,2}} + 2{\pi _{2,2}} + 3{\pi _{3,2}} + \cdots + {\pi _{1,M}} + 2{\pi _{2,M}} + 3{\pi _{3,M}} + \cdots ) \hfill \\ {\text{ }} = {\lambda _{AA}}({{\boldsymbol{\pi }}_1}{\bf{e}} + 2{{\boldsymbol{\pi }}_2}{\bf{e}} + 3{{\boldsymbol{\pi }}_3}{\bf{e}} + \cdots ) \hfill \\ {\text{ }} = {\lambda _{AA}}\sum\limits_{k = 1}^\infty {k{{\boldsymbol{\pi }}_k}} \hfill \\ {\text{ }} = {\lambda _{AA}}E[{N_A}]. \hfill \\ \end{gathered} (29)

    iv. Average infection rate from symptomatic patients to asymptomatic patients

    \begin{gathered} {r_{A \to B}} = {\lambda _{AB}}({\pi _{1,0}} + 2{\pi _{2,0}} + 3{\pi _{3,0}} + \cdots + {\pi _{1,1}} + 2{\pi _{2,1}} + 3{\pi _{3,1}} + \cdots + \hfill \\ {\text{ }}{\pi _{1,2}} + 2{\pi _{2,2}} + 3{\pi _{3,2}} + \cdots + {\pi _{1,M}} + 2{\pi _{2,M}} + 3{\pi _{3,M}} + \cdots ) \hfill \\ {\text{ }} = {\lambda _{AB}}({{\boldsymbol{\pi }}_1}{\bf{e}} + 2{{\boldsymbol{\pi }}_2}{\bf{e}} + 3{{\boldsymbol{\pi }}_3}{\bf{e}} + \cdots ) \hfill \\ {\text{ }} = {\lambda _{AB}}\sum\limits_{k = 1}^\infty {k{{\boldsymbol{\pi }}_k}} \hfill \\ {\text{ }} = {\lambda _{AB}}E[{N_A}]. \hfill \\ \end{gathered} (30)

    v. Average infection rate from asymptomatic patients to symptomatic patients

    \begin{gathered} {r_{B \to A}} = {\lambda _{BA}}({\pi _{0,1}} + {\pi _{1,1}} + {\pi _{2,1}} + \cdots + 2{\pi _{0,2}} + 2{\pi _{1,2}} + 2{\pi _{2,2}} + \cdots + \hfill \\ {\text{ 3}}{\pi _{0,3}} + 3{\pi _{1,3}} + 3{\pi _{2,3}} + \cdots + M{\pi _{0,M}} + M{\pi _{1,M}} + M{\pi _{2,M}} + \cdots ) \hfill \\ {\text{ }} = {\lambda _{BA}}({{\boldsymbol{\pi }}_0}{{\bf{f}}_B} + {{\boldsymbol{\pi }}_1}{{\bf{f}}_B} + {{\boldsymbol{\pi }}_2}{{\bf{f}}_B} + \cdots ) \hfill \\ {\text{ }} = {\lambda _{BA}}\sum\limits_{k = 0}^\infty {{{\boldsymbol{\pi }}_k}{{\bf{f}}_B}} \hfill \\ {\text{ }} = {\lambda _{BA}}E[{N_B}]. \hfill \\ \end{gathered} (31)

    vi. Average infection rate from asymptomatic patients to asymptomatic patients

    \begin{gathered} {r_{B \to B}} = {\lambda _{BB}}({\pi _{0,1}} + {\pi _{1,1}} + {\pi _{2,1}} + \cdots + 2{\pi _{0,2}} + 2{\pi _{1,2}} + 2{\pi _{2,2}} + \cdots + \hfill \\ {\text{ 3}}{\pi _{0,3}} + 3{\pi _{1,3}} + 3{\pi _{2,3}} + \cdots + \left( {M - 1} \right){\pi _{0,M}} + \left( {M - 1} \right){\pi _{1,M}} + \left( {M - 1} \right){\pi _{2,M}} + \cdots ) \hfill \\ {\text{ }} = {\lambda _{BB}}({{\boldsymbol{\pi }}_0}{{\bf{g}}_B} + {{\boldsymbol{\pi }}_1}{{\bf{g}}_B} + {{\boldsymbol{\pi }}_2}{{\bf{g}}_B} + \cdots ) \hfill \\ {\text{ }} = {\lambda _{BB}}\sum\limits_{k = 0}^\infty {{{\boldsymbol{\pi }}_k}{{\bf{g}}_B}} , \hfill \\ \end{gathered} (32)

    where, {{\bf{g}}_B} = {\left({0, 1, 2, \ldots, M - 1, 0} \right)^T} .

    vii. Average transition rate from asymptomatic patients to symptomatic patients

    \begin{gathered} {r_{BA}} = \gamma ({\pi _{0,1}} + {\pi _{1,1}} + {\pi _{2,1}} + \cdots + 2{\pi _{0,2}} + 2{\pi _{1,2}} + 2{\pi _{2,2}} + \cdots + \hfill \\ {\text{ 3}}{\pi _{0,3}} + 3{\pi _{1,3}} + 3{\pi _{2,3}} + \cdots + M{\pi _{0,M}} + M{\pi _{1,M}} + M{\pi _{2,M}} + \cdots ) \hfill \\ {\text{ }} = \gamma ({{\boldsymbol{\pi }}_0}{{\bf{f}}_B} + {{\boldsymbol{\pi }}_1}{{\bf{f}}_B} + {{\boldsymbol{\pi }}_2}{{\bf{f}}_B} + \cdots ) \hfill \\ {\text{ }} = \gamma \sum\limits_{k = 0}^\infty {{{\boldsymbol{\pi }}_k}{{\bf{f}}_B}} \hfill \\ {\text{ }} = \gamma E[{N_B}]. \hfill \\ \end{gathered} (33)

    This section provides two efficient algorithms for computing the stationary probability vector of the QBD process of infectious disease. And we use this stationary probability vector to compute the important performance measures of the QBD process.

    Given the infinite state space of the QBD process, it is critical to first identify the maximum truncation level and the R-measure sequence \{ {R_k}, k \geqslant 0\} in order to obtain the stationary probability vector ( {\boldsymbol{\pi }} ). For this purpose, we propose an algorithm based on the method of Bright and Taylor [52], which enables the determination of the maximum truncation level K and all R-measure \{ {R_k}, k \geqslant 0\} .

    Based on Bright and Taylor [52], we derive two key results as follows:

    a) An important UD-type representation of {R_k} .

    {R_k} = \sum\limits_{l = 0}^\infty {U_k^l} \prod\limits_{i = 0}^{l - 1} {D_{k + {2^{l - i}}}^{l - 1 - i}} ,k \geqslant 1, (34)

    where

    U_k^0 = {A_{k,k + 1}}{( - {A_{k + 1,k + 1}})^{ - 1}},k \geqslant 1, (35)
    D_k^0 = {A_{k - 1,k}}{( - {A_{k - 1,k - 1}})^{ - 1}},k \geqslant 1, (36)
    U_k^{l + 1} = U_k^lU_{k + {2^l}}^l{[I - U_{k + {2^{l + 1}}}^lD_{k + 3 \cdot {2^l}}^l - D_{k + {2^{l + 1}}}^lU_{k + {2^l}}^l]^{ - 1}},k \geqslant 0, (37)
    D_k^{l + 1} = D_k^lD_{k - {2^l}}^l{[I - U_{k - {2^{l + 1}}}^lD_{k - {2^l}}^l - D_{k - {2^{l + 1}}}^lU_{k - 3 \cdot {2^l}}^l]^{ - 1}},k \geqslant 0. (38)

    Note that, using this UD-type representation, the rate matrix \left\{ {{R_k}, k \geqslant 1} \right\} is determined directly from the original data matrices {A_{k, k + 1}} , {A_{k + 1, k + 1}} , {A_{k - 1, k}} , and {A_{k - 1, k - 1}} , k \geqslant 1 .

    b) A backward iterative recursive representation of {R_k} .

    As soon as {R_K} is obtained from the UD-type representation, we have

    {R_k} = {A_{k,k + 1}}{( - {A_{k + 1,k + 1}} - {R_{k + 1}}{A_{k + 1,k + 2}})^{ - 1}},k = K - 1,K - 2, \ldots ,1,0. (39)

    Based on a) and b), we propose an algorithm to determine the appropriate maximum truncation level K and obtain the R-measure \left\{ {{R_k}, 0 \leqslant k \leqslant K} \right\} .

    To determine the maximum truncation level K , we iteratively compute the sum of the stationary probabilities of the level K until it is sufficiently small, i.e., {\pi _K}{\bf{e}} < \varepsilon . \varepsilon denotes the controllable accuracy by the user, where we make \varepsilon = {10^{ - 20}} . In addition, we define a sequence of positive integers \{ {\varsigma _k}:k \geqslant 0\} for iterative computation. The algorithm for determining the maximum truncation level and the R-measure is summarized in Algorithm 1.

    Algorithm 1: Calculate K and {R_k}
    Input: Parameter {\lambda _{AA}} , {\lambda _{AB}} , {\lambda _{BA}} , {\lambda _{BB}} , \gamma , {\mu _A} and {\mu _B} .
    Output: K and {R_k}
    1) Step 0: Initialization
    Let {\varsigma _n} = 10n and initialize setting n =1 as the initial value.
    2) Step 1: Calculate the rate matrix {R_K}
    Let K = {\varsigma _n} , and calculate {R_K} according to Eq (34).
    3) Step 2: Determine other rate matrices ( \{ {R_k}, 0 \leqslant k \leqslant K - 1\} ) by backward iteration
    Compute of {R_{K - 1}} , {R_{K - 2}} , …, {R_2} , {R_1} , {R_0} according to Eq (39).
    4) Step 3: Compute the stationary probability vector {{\boldsymbol{\tilde \pi }}_0}
    Compute {{\boldsymbol{\tilde \pi }}_0} by solving the linear equation system (26)
    5) Step 4: Calculate the normalization constant c
    Using the vector {{\boldsymbol{\tilde \pi }}_0} and the R-measure \{ {R_l}, 0 \leqslant l \leqslant K\} , we can compute
    c = \frac{1}{{{{{\boldsymbol{\tilde \pi }}}_0}{\bf{e}} + \sum\limits_{k = 1}^{K + 1} {{{{\boldsymbol{\tilde \pi }}}_0}{R_0}{R_1} \cdots {R_{k - 1}}{\bf{e}}} }}.\;\;\;\;\;(40)
    6) Step 5: Check if the accuracy is satisfied
    If for a positive integer K = {\varsigma _n} , it satisfies
    {{\boldsymbol{\pi }}_{K + 1}}{\bf{e}} < \varepsilon ,\;\;\;\;(41)
    (called stop condition), where {{\boldsymbol{\pi }}_{K + 1}} = c{{\boldsymbol{\tilde \pi }}_0}{R_o}{R_1} \cdots {R_K} , then let K = {\varsigma _n} , and go to Step 6. Otherwise, let n = n + 1 and substitute to calculate {\varsigma _n} = 10(n + 1) , return to Step 1;
    7) Step 6: Output K and \{ {R_k}, 0 \leqslant k \leqslant K\}
    The algorithm stops, and a maximum truncation level K = {\varsigma _n} and a sequence of R-measures \{ {R_k}, 0 \leqslant k \leqslant K\} are obtained as its output.

     | Show Table
    DownLoad: CSV

    Once we have determined the maximum truncation level K and R-measure \{ {R_k}, 0 \leqslant k \leqslant K\} using Algorithm 1, we need to further compute the stationary probability vector and some other important performance measures of the model. Therefore, we design Algorithm 2 as follows.

    Algorithm 2: Compute stationary probability vector {\boldsymbol{\pi }} and performance measures
    Input: Parameter K , \{ {R_l}, 0 \leqslant l \leqslant K\} , {{\boldsymbol{\tilde \pi }}_0} , and c .
    Output: E[{N_A}] , E[{N_B}] , {r_{A \to A}} , {r_{A \to B}} , {r_{B \to A}} , {r_{B \to B}} and {r_{BA}}
    1) Step 1: Calculate the stationary probability vector
    When 1 \leqslant k \leqslant K + 1 , based on Eqs (23) and (24), we can compute
    {{\boldsymbol{\pi }}_k} = c{{\boldsymbol{\tilde \pi }}_0}{R_0}{R_1} \cdots {R_{k - 1}}, (42)
    and according to Eq (26), we can calculate
    {{\boldsymbol{\pi }}_0} = c{{\boldsymbol{\tilde \pi }}_0}. (43)
    2) Step 2: Calculate the important performance measures
    We use Eqs (44)–(50) in Algorithm 2 to calculate the measures E[{N_A}] , E[{N_B}] , {r_{A \to A}} , {r_{A \to B}} , {r_{B \to A}} , {r_{B \to B}} , {r_{BA}} .
    E[{N_A}] = \sum\limits_{l = 0}^K {l{{\boldsymbol{\pi }}_l}{\bf{e}}}, (44)
    E[{N_B}] = \sum\limits_{l = 0}^K {{{\boldsymbol{\pi }}_l}{{\bf{f}}_B}}, (45)
    {r_{A \to A}} = E[{N_A}] \cdot {\lambda _{AA}}, (46)
    {r_{A \to B}} = E[{N_A}] \cdot {\lambda _{AB}}, (47)
    {r_{B \to A}} = E[{N_B}] \cdot {\lambda _{BA}}, (48)
    {r_{B \to B}} = \sum\limits_{k = 0}^K {{{\boldsymbol{\pi }}_k}{{\bf{g}}_B}} \cdot {\lambda _{BB}}, (49)
    {r_{BA}} = E[{N_B}] \cdot \gamma. (50)
    3) Step 3: Output the important performance measures
    Stop algorithm, output these performance measures: E[{N_A}] , E[{N_B}] , {r_{A \to A}} , {r_{A \to B}} , {r_{B \to A}} , {r_{B \to B}} , {r_{BA}} .

     | Show Table
    DownLoad: CSV

    In this section, we give some numerical examples to demonstrate the feasibility of the proposed model in this paper and highlight the impact of key parameters on some important system performance measures.

    We conduct numerical experiments using COVID-19 as an example to verify the feasibility of the model. As obtaining actual data about asymptomatic patients can be challenging from publicly available information, we set the benchmark values of the model parameters based on previous studies. Specifically, He et al. [15] reported that the infection risk ratio between symptomatic and asymptomatic patients is approximately 3.9:1, i.e., ({\lambda _{AA}} + {\lambda _{AB}})/({\lambda _{BA}} + {\lambda _{BB}}) = 3.9/1 . In addition, about 85% of patients infected by symptomatic patients are symptomatic themselves, leading to {\lambda _{AA}}/{\lambda _{AB}} = 0.85/0.15 . Furthermore, roughly 50% of patients infected by asymptomatic patients become symptomatic, i.e., {\lambda _{BA}}/{\lambda _{BB}} = 0.5/0.5 . According to the results of Lauer et al. [18], the average duration for asymptomatic patients to become symptomatic is 5 days, which we make the parameter \gamma = 1/5 . We assume that these values remain constant throughout the observation period, which corresponds to the homogeneous Markov process.

    After consulting relevant experts, referencing the research findings from literature [15,18], and taking into account the system stability conditions, we set the parameter values to be as follows: {\tau _A} = 1.5 , {\tau _B} = 3 , {\lambda _{AA}} = 3.4 , {\lambda _{AB}} = 0.6 , {\lambda _{BB}} = 0.51 , {\lambda _{BA}} = 0.51 , \gamma = 0.2 , {\mu _A} = 6 and {\mu _B} = 1 . Subsequently, we give some numerical experiments to investigate the impact of various parameters on the average number of symptomatic and asymptomatic patients. We also analyze the effect of non-drug interventions and vaccines. Note that, the parameters {r_{A \to A}} , {r_{A \to B}} , {r_{B \to A}} , {r_{B \to B}} and {r_{BA}} are closely related to E[{N_A}] and E[{N_B}] according to Eqs (29)–(33), thus we only observe E[{N_A}] and E[{N_B}] in our numerical experiments.

    1) Case 1: Let's take {\lambda _{BB}} as 0.3, 0.5 and 0.7, respectively, and observe the impact of {\lambda _{AA}} varying in the interval [2.9, 3.9] on E[{N_A}] and E[{N_B}] . The numerical results are shown in Figure 5.

    Figure 5.  The impact of {\lambda _{AA}} on the average number of COVID-19 patients.

    Figure 5 shows the impact of the infection rate of symptomatic patients infecting and generating symptomatic patients ( {\lambda _{AA}} ) on the average number of symptomatic patients and asymptomatic patients when the system is stable. The blue, red and green lines in Figure 5 represent the cases in which the rate of asymptomatic patients infecting and generating asymptomatic patients ( {\lambda _{BB}} ) is set as 0.3, 0.5 and 0.7, respectively. In Figure 5(a), each line shows that the average number of symptomatic patients increases as {\lambda _{AA}} increases. In addition, by comparing the three lines, we can observe that a larger {\lambda _{BB}} leads to a higher average number of symptomatic patients. Similarly, Figure 5(b) shows that the average number of asymptomatic patients also increases as {\lambda _{AA}} and {\lambda _{BB}} increase.

    Based on the results shown in Figure 5, it can be concluded that the proposed method is feasible for modeling the long-term infectious disease transmission process. We can estimate the average number of symptomatic and asymptomatic patients by using our proposed model. The numerical experiments demonstrate that decreasing the infection rate of both symptomatic and asymptomatic patients can decrease the average number of symptomatic and asymptomatic patients. This suggests that certain non-drug interventions, such as implementing social distancing measures, can effectively control the transmission of infectious diseases by reducing the infection rates.

    2) Case 2: Let's take {\mu _B} as 1.1, 1.3 and 1.5, respectively, and observe the impact of {\mu _A} varying in the interval [5,6] on E[{N_A}] and E[{N_B}] . The numerical results are shown in Figure 6.

    Figure 6.  The effect of {\mu _A} on the average number of COVID-19 patients.

    Figure 6 shows the impact of the disappearance rate of symptomatic patients ( {\mu _A} ) on the average number of symptomatic patients and asymptomatic patients. The blue, red and green lines in Figure 6 represent the cases in which the disappearance rate of asymptomatic patients ( {\mu _B} ) is set as 1.1, 1.3 and 1.5, respectively. Each line in Figure 6(a) shows that the average number of symptomatic patients decreases as {\mu _A} increase. In addition, comparing the three lines shows that the larger the value of {\mu _B} , the smaller the average number of symptomatic patients. Similarly, Figure 6(b) shows that the average number of asymptomatic patients also decreases as {\mu _A} and {\mu _B} increase.

    Based on the results in Figure 6, increasing the disappearance rate of symptomatic and asymptomatic patients can result in a decrease in the average number of symptomatic and asymptomatic patients. This shows the effectiveness of non-drug interventions, such as screening and isolation, in controlling the spread of infectious diseases by increasing the disappearance rates.

    3) Case 3: Let's take {\mu _B} as 1.1, 1.3 and 1.5, respectively, and observe the impact of \gamma varying in the interval [0, 1] on E[{N_A}] and E[{N_B}] . The numerical results are shown in Figure 7.

    Figure 7.  The impact of \gamma on the average number of COVID-19 patients.

    Figure 7 shows the trend in the average number of symptomatic patients and asymptomatic patients when the transition rate from asymptomatic patients to symptomatic patients ( \gamma ) increases. The blue, red and green lines in Figure 7 represent the cases where the disappearance rate of asymptomatic patients ( {\mu _B} ) is given as 1.1, 1.3 and 1.5, respectively. Each line in Figure 7(a) indicates that the average number of symptomatic patients increases as \gamma increases. Furthermore, comparing the three lines reveals that the larger the value of {\mu _B} , the smaller the average number of symptomatic patients. In contrast, Figure 7(b) demonstrates that the average number of asymptomatic patients decreases as \gamma and {\mu _B} increase.

    According to the results in Figure 7, increasing the transition rate from asymptomatic patients to symptomatic patients leads to an increase in the average number of symptomatic patients and a decrease in the average number of asymptomatic patients. Due to the higher infection rate of symptomatic patients compared to asymptomatic patients, and the longer duration of infection caused by asymptomatic patients, the impact of increasing the transition rate on controlling the transmission of infectious diseases may vary depending on the specific circumstances.

    4) Case 4: Let's take {\mu _A} as 5, 6 and 7, respectively, and observe the impact of {\mu _B} varying in the interval [0, 1] on E[{N_A}] and E[{N_B}] . The numerical results are shown in Figure 8.

    Figure 8.  The impact of {\mu _B} on the average number of COVID-19 patients.

    Figure 8 illustrates the trend in the average number of symptomatic patients and asymptomatic patients as the disappearance rate of asymptomatic patients ( {\mu _B} ) increases. The blue, red and green lines in Figure 8 represent the cases in which the disappearance rate of asymptomatic patients ( {\mu _A} ) is given as 5, 6 and 7, respectively. Each line in Figure 8(a) shows that the average number of symptomatic patients decreases as {\mu _B} increases. The same trend is observed in Figure 8(b), where the average number of asymptomatic patients decreases as {\mu _B} increases. Furthermore, the decline in the average number of symptomatic and asymptomatic patients is notably faster within the range of 0.4–0.6 compared to other values.

    Based on the findings in Figure 8, increasing the disappearance rate of asymptomatic patients results in a decline in the average number of symptomatic and asymptomatic patients. Notably, there exists a specific range where increasing the disappearance rate of asymptomatic patients can lead to a more substantial reduction in the average number of patients compared to other scenarios. Therefore, it is advisable to identify this range prior to implementing non-drug interventions, as it can contribute to more cost-effective control of the infectious disease transmission process.

    5) Case 5: In order to examine the impact of non-drug intervention strategies and vaccines on infectious disease transmission, we introduce two variables: w and v , representing the effects of non-drug interventions and drug interventions (e.g., vaccines) on the infection rate of disease, respectively. Additionally, let {\beta _1} and {\beta _2} be the weights representing the effects of non-drug interventions and drug interventions on the infection rate of infectious disease, respectively. Obviously, {\beta _1} + {\beta _2} = 1 . Based on this, we have

    \lambda ' = ({\beta _1}\frac{1}{w} + {\beta _2}\frac{1}{v})\lambda (51)

    where \lambda ' denotes the infection rate considering the impact of both non-drug interventions and drug interventions, and \lambda denotes the infection rate without considering them. In this situation, we set the benchmark parameter value for the infection rate without considering the impact of non-drug interventions and drug interventions. Under different situations, we examine the impact of v on the COVID-19 transmission in three cases: {\beta _1} = 0.5 , {\beta _2} = 0.5 ; {\beta _1} = 0.2 , {\beta _2} = 0.8 ; and {\beta _1} = 0.8 , {\beta _2} = 0.2 .

    Figure 9 gives the impact of v on the E[{N_A}] and E[{N_B}] under different w in various weight combination cases. The blue, red and green lines in Figure 9 represent the cases where the effect of non-drug interventions on the infection rate of disease ( w ) is set as 1, 1.5 and 2, respectively. Each line in Figure 9(a) shows that the average number of symptomatic patients decreases as v increases. Furthermore, comparing the three lines reveals that the average number of symptomatic patients decreases as w increases. The same trend is observed in Figure 9(b), where the average number of asymptomatic patients also decreases as v and w increase. According to the results in Figure 9, both non-drug interventions and drug interventions contribute to a decrease in the average number of symptomatic and asymptomatic patients, thus effectively controlling the transmission of infectious disease.

    Figure 9.  The impact of v on the average number of COVID-19 patients.

    Based on the stability and the numerical analysis, we can know that the infectious disease transmission process with asymptomatic patients is stable when {\lambda _{AA}} < {\mu _A} . Controlling the infection rates and disappearance rates of symptomatic and asymptomatic patients can affect the average number of symptomatic and asymptomatic patients when this infectious disease transmission process is stable. Therefore, while paying attention to the control of symptomatic patients, the screening and control of asymptomatic patients should not be neglected. Specifically, the average number of symptomatic patients increases with {\lambda _{AA}} , {\lambda _{AB}} , {\lambda _{BA}} , {\lambda _{BB}} and \gamma ; and decreases with {\mu _A} and {\mu _B} . The average number of asymptomatic patients increases with {\lambda _{AA}} , {\lambda _{AB}} , {\lambda _{BA}} and {\lambda _{BB}} ; and decreases with \gamma , {\mu _A} and {\mu _B} . Therefore, if the control strategies can decrease the parameters {\lambda _{AA}} , {\lambda _{AB}} , {\lambda _{BA}} and {\lambda _{BB}} or increase the parameters {\mu _A} and {\mu _B} , it can be helpful for the scientific control and management of the epidemic. Increasing \gamma is effective for controlling transmission caused by asymptomatic patients, but not for transmission caused by symptomatic patients. Furthermore, there exists an optimal interval in which increasing {\mu _B} can more effectively control the transmission of infectious diseases. Therefore, decisions on how to control \gamma and {\mu _B} should be made based on the specific situation. Both non-drug interventions and drug interventions can influence the parameters of the infectious disease transmission. The above findings are in line with those in the existing related literature [26,27,28,30]. In addition, we consider the impact of transmission from asymptomatic patients on the system. Based on the research results of this paper and the specific situation of the region, it is recommended to regulate the intensity of infectious disease control strategies and the level of drug interventions at an optimal level to achieve effective and cost-efficient control of infectious diseases.

    In this paper, we developed a new Markov process of infectious disease with infinite-state space to study the stable condition of an infectious disease system, which is based on the continuous time level-dependent QBD process. Compared to existing studies, we considered the long-term impact of asymptomatic patients on the transmission process of infectious diseases. We gave a simple sufficient condition for the stability of the model by using the mean drift technique. We obtained the stationary probability vector by means of RG-factorizations. Using this stationary probability vector, we gave some important performance measures of this model. Finally, we verified the feasibility of this paper by using some numerical experiments through taking COVID-19 as an example. The effects of some key parameters on the stationary probability of the model and the role of two types of patients on the transmission of infectious diseases were also quantitatively analyzed based on performance analysis. The modeling technique proposed in this paper is operational and general. It can be applied to other infectious diseases if they consist of symptomatic and asymptomatic patients and both types of patients can infect and generate the two types of patients. Therefore, the methodology and results of this paper can provide some useful theoretical basis and technical support for the long-term scientific prevention and control of infectious diseases.

    In future studies, the proposed model in this paper can be extended in the following ways:

    1) In the framework of the homogeneous Markov process, the patient infection and disappearance processes can be extended from Poisson processes to Markovian arrival processes (MAPS), whose purpose is to reflect the burstiness and short-time clustering characteristics of infectious diseases.

    2) Extending the patient infection and treatment process from a homogeneous Poisson process to a non-homogeneous Poisson process, i.e., the parameters can vary over time.

    3) Develop a data-driven Markov process by deriving the infection, cure and death rates of patients in the model based on the actual data of infectious diseases.

    4) Apply Markov decision process to study optimal control strategies in infectious disease prevention and control, e.g., vaccine input strategy, isolation intensity strategy, economic strategy and etc.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The data used in numerical experiments have been uploaded to the database of Zenodo (https://doi.org/10.5281/zenodo.8280820).

    This work is sponsored by the National Natural Science Foundation of China under Grants 72202010, and the China Postdoctoral Science Foundation (2022M710275).

    The authors declare there is no conflict of interest.

    In this appendix, we use the mean drift technique to provide a simple sufficient condition for the stability of the QBD process. Ensuring model stability is essential for the analysis of the model.

    Theorem A1. When {\mu _A} > {\lambda _{AA}} , this QBD process of infectious disease ( Q ) is irreducible and positive recurrent. Thus, the QBD process of infectious disease is stable.

    Proof. Let {{\bf{A}}_k} = {A_{k,k - 1}} + {A_{k,k}} + {A_{k,k + 1}} , where k denotes the k th level, k \geqslant 1 . Let {h_i} = k{\lambda _{AB}} + i(\gamma + {\mu _B} + {\lambda _{BB}}) , {o_i} = k{\lambda _{AB}} + i{\lambda _{BB}} . Then, {{\bf{A}}_k} can be written as follows:

    {{\bf{A}}_k} = \left( {\begin{array}{*{20}{c}} { - {h_0}}&{{o_0}}&{}&{}&{}&{} \\ {\gamma + {\mu _B}}&{ - {h_1}}&{{o_1}}&{}&{}&{} \\ {}&{2(\gamma + {\mu _B})}&{ - {h_2}}&{{o_2}}&{}&{} \\ {}&{}& \ddots & \ddots & \ddots &{} \\ {}&{}&{}&{(M - 1)(\gamma + {\mu _B})}&{ - {h_{M - 1}}}&{{o_{M - 1}}} \\ {}&{}&{}&{}&{M(\gamma + {\mu _B})}&{ - {h_M}} \end{array}} \right).

    Note that it is easy to see that {{\bf{A}}_k} is a finite state and irreducible birth-death process, so this process must be positive recurrent and stable. Therefore, this birth-death process {{\bf{A}}_k} has the stable state and stationary probability vector.

    Let {{\boldsymbol{\alpha }}^{(k)}} = (\alpha _0^{(k)},\alpha _1^{(k)},\alpha _2^{(k)}, \ldots ,\alpha _M^{(k)}) be the stationary probability vector of {{\bf{A}}_k} . Then, from the characteristics of the stationary probability vector, we can obtain

    {{\boldsymbol{\alpha }}^{(k)}}{{\bf{{A}}}_k} = {{\boldsymbol{\alpha }}^{(k)}}({A_{k,k - 1}} + {A_{k,k}} + {A_{k,k + 1}}) = 0, (A1)
    {{\boldsymbol{\alpha }}^{(k)}}{\bf{e}} = 1, (A2)

    where {\bf{e}} is a column vector of ones with a suitable size.

    According to Eq (A1), we have

    \begin{gathered} {{\boldsymbol{\alpha }}^{(k)}}{{\bf{{A}}}_k} = \left( {\alpha _0^{(k)}( - k{\lambda _{AB}}) + \alpha _1^{(k)}(\gamma + {\mu _B}),\alpha _0^{(k)}k{\lambda _{AB}} + \alpha _1^{(k)}\left[ { - k{\lambda _{AB}} - \left( {{\mu _B} + \gamma + {\lambda _{BB}}} \right)} \right] + \alpha _2^{(k)}2(\gamma + {\mu _B}),} \right. \hfill \\ {\text{ }}\alpha _1^{(k)}\left( {k{\lambda _{AB}} + {\lambda _{BB}}} \right) + \alpha _2^{(k)}\left[ { - k{\lambda _{AB}} - 2\left( {{\mu _B} + \gamma + {\lambda _{BB}}} \right)} \right] + \alpha _3^{(k)}3(\gamma + {\mu _B}), \ldots , \hfill \\ {\text{ }}\alpha _{M - 2}^{(k)}\left[ {k{\lambda _{AB}} + \left( {M - 2} \right){\lambda _{BB}}} \right] + \alpha _{M - 1}^{(k)}\left[ { - k{\lambda _{AB}} - (M - 1)\left( {{\mu _B} + \gamma + {\lambda _{BB}}} \right)} \right] + \alpha _M^{(k)}M(\gamma + {\mu _B}), \hfill \\ {\text{ }}\left. {\alpha _{M - 1}^{(k)}\left[ {k{\lambda _{AB}} + \left( {M - 1} \right){\lambda _{BB}}} \right] + \alpha _M^{(k)}\left[ { - k{\lambda _{AB}} - M\left( {{\mu _B} + \gamma + {\lambda _{BB}}} \right)} \right]} \right). \hfill \\ \end{gathered} (A3)

    From Eq (A1) and the first term in Eq (A3), we have

    \alpha _0^{(k)}( - k{\lambda _{AB}}) + \alpha _1^{(k)}(\gamma + {\mu _B}) = 0,
    \alpha _1^{(k)} = \frac{{k{\lambda _{AB}}}}{{\gamma + {\mu _B}}}\alpha _0^{(k)}. (A4)

    From Eq (A1), the second term in Eqs (A3) and (A4), we have

    \alpha _0^{(k)}k{\lambda _{AB}} + \alpha _1^{(k)}\left[ { - k{\lambda _{AB}} - \left( {{\mu _B} + \gamma + {\lambda _{BB}}} \right)} \right] + \alpha _2^{(k)}2(\gamma + {\mu _B}) = 0,
    \alpha _2^{(k)} = \frac{{k{\lambda _{AB}}\left( {k{\lambda _{AB}} + {\lambda _{BB}}} \right)}}{{2!{{\left( {\gamma + {\mu _B}} \right)}^2}}}\alpha _0^{(k)}.

    Similarly, we can obtain Eqs (A5) and (A6) to get \alpha _l^{(k)} .

    \alpha _{l - 2}^{(k)}\left( {k{\lambda _{AB}} + \left( {l - 2} \right){\lambda _{BB}}} \right) + \alpha _{l - 1}^{(k)}\left[ { - k{\lambda _{AB}} - \left( {l - 1} \right)\left( {{\mu _B} + \gamma + {\lambda _{BB}}} \right)} \right] + \alpha _l^{(k)}l(\gamma + {\mu _B}) = 0,l \geqslant 1, (A5)
    \alpha _l^{(k)} = \frac{{\prod\limits_{i = 0}^{l - 1} {\left( {k{\lambda _{AB}} + i{\lambda _{BB}}} \right)} }}{{l!{{\left( {\gamma + {\mu _B}} \right)}^l}}}\alpha _0^{(k)},l \geqslant 1. (A6)

    After that, depending on Eqs (A2) and (A6), we can obtain

    \alpha _0^{(k)} + \frac{{k{\lambda _{AB}}}}{{\gamma + {\mu _B}}}\alpha _0^{(k)} + \frac{{\prod\limits_{i = 0}^{l - 1} {\left( {k{\lambda _{AB}} + i{\lambda _{BB}}} \right)} }}{{2!{{\left( {\gamma + {\mu _B}} \right)}^2}}}\alpha _0^{(k)} + \cdots + \frac{{\prod\limits_{i = 0}^{M - 1} {\left( {k{\lambda _{AB}} + i{\lambda _{BB}}} \right)} }}{{M!{{\left( {\gamma + {\mu _B}} \right)}^M}}}\alpha _0^{(k)} = 1,
    \alpha _0^{(k)} = \frac{1}{{1 + \sum\limits_{j = 1}^M {\frac{{\prod\limits_{i = 0}^{j - 1} {\left( {k{\lambda _{AB}} + i{\lambda _{BB}}} \right)} }}{{j!{{\left( {\gamma + {\mu _B}} \right)}^j}}}} }},l \geqslant 1.

    Using the aforementioned description, we can obtain the stationary probability vector {{\boldsymbol{\alpha }}^{(k)}} for the birth-death process {{\bf{{A}}}_k} . Subsequently, we can leverage this vector to discuss the stable conditions of the QBD process proposed in this paper. In particular, we analyze two key rates when discussing the stable state: the upward mean drift rate and the downward mean drift rate.

    The upward mean drift rate: mean drift rate from level k to level k + 1

    \begin{gathered} {{\boldsymbol{\alpha }}^{(k)}}{A_{k,k + 1}}{\bf{e}} = {{\boldsymbol{\alpha }}^{(k)}}{(k{\lambda _{AA}},k{\lambda _{AA}} + ({\lambda _{BA}} + \gamma ), \ldots ,k{\lambda _{AA}} + M({\lambda _{BA}} + \gamma ))^T} \hfill \\ {\text{ }} = \alpha _0^{(k)}k{\lambda _{AA}} + \alpha _1^{(k)}\left[ {k{\lambda _{AA}} + \left( {{\lambda _{BA}} + \gamma } \right)} \right] + \cdots + \alpha _M^{(k)}\left[ {k{\lambda _{AA}} + M\left( {{\lambda _{BA}} + \gamma } \right)} \right] \hfill \\ {\text{ }} = k{\lambda _{AA}}(\alpha _0^{(k)} + \alpha _0^{(k)} + \cdots \alpha _0^{(k)}) + \left[ {\alpha _1^{(k)}\left( {{\lambda _{BA}} + \gamma } \right) + \alpha _2^{(k)}2\left( {{\lambda _{BA}} + \gamma } \right) + \cdots + \alpha _M^{(k)}M\left( {{\lambda _{BA}} + \gamma } \right)} \right] \hfill \\ {\text{ }} = k{\lambda _{AA}} + \left[ {\alpha _1^{(k)}\left( {{\lambda _{BA}} + \gamma } \right) + \alpha _2^{(k)}2\left( {{\lambda _{BA}} + \gamma } \right) + \cdots + \alpha _M^{(k)}M\left( {{\lambda _{BA}} + \gamma } \right)} \right]. \hfill \\ \end{gathered}

    The downward mean drift rate: mean drift rate from level k to level k - 1

    \begin{gathered} {{\boldsymbol{\alpha }}^{(k)}}{A_{k,k - 1}}{\bf{e}} = {{\boldsymbol{\alpha }}^{(k)}}{(k{\mu _A},k{\mu _A}, \ldots ,k{\mu _A})^T} \hfill \\ {\text{ }} = \alpha _0^{(k)}k{\mu _A} + \alpha _1^{(k)}k{\mu _A} + \cdots + \alpha _M^{(k)}k{\mu _A} \hfill \\ {\text{ }} = k{\mu _A}(\alpha _0^{(k)} + \alpha _0^{(k)} + \cdots \alpha _0^{(k)}) \hfill \\ {\text{ }} = k{\mu _A}. \hfill \\ \end{gathered} (A7)

    According to the mean drift technique, a QBD process is positive recurrent if its downward mean drift rate is greater than its upward mean drift rate. Therefore, we proceed to calculate

    \begin{gathered} {{\boldsymbol{\alpha }}^{(k)}}{A_{k,k - 1}}{\bf{e}} - {{\boldsymbol{\alpha }}^{(k)}}{A_{k,k + 1}}{\bf{e}} = k{\mu _A} - k{\lambda _{AA}} + \left[ {\alpha _1^{(k)}\left( {{\lambda _{BA}} + \gamma } \right) + \alpha _2^{(k)}2\left( {{\lambda _{BA}} + \gamma } \right) + \cdots + \alpha _M^{(k)}M\left( {{\lambda _{BA}} + \gamma } \right)} \right] \hfill \\ {\text{ }} = k\left( {{\mu _A} - {\lambda _{AA}}} \right) + \left[ {\alpha _1^{(k)}\left( {{\lambda _{BA}} + \gamma } \right) + \alpha _2^{(k)}2\left( {{\lambda _{BA}} + \gamma } \right) + \cdots + \alpha _M^{(k)}M\left( {{\lambda _{BA}} + \gamma } \right)} \right]. \hfill \\ \end{gathered} (A8)

    In Eq (A8), M is a fixed value, so \alpha _1^{(k)}\left( {{\lambda _{BA}} + \gamma } \right) + \alpha _2^{(k)}2\left( {{\lambda _{BA}} + \gamma } \right) + \cdots + \alpha _M^{(k)}M\left( {{\lambda _{BA}} + \gamma } \right) is a finite value. If {\mu _A} - {\lambda _{AA}} > 0 and k tends to positive infinity, k({\mu _A} - {\lambda _{AA}}) will tends to positive infinity. Therefore, when k tends to positive infinity, Eq (A8) must be larger than 0. This means that the downward mean drift rate of the QBD process Q is greater than its upward mean drift rate. Therefore, Q is positive recurrent and stable. This completes the proof.



    [1] S. Dunin-Horkawicz, A. Czerwoniec, M. J. Gajda, et al., MODOMICS: A database of RNA modification pathways, Nucleic Acids Res., 34(2006), D145–D149.
    [2] J. H. Ge and Y. T. Yu, RNA pseudouridylation: New insights into an old modification, Trends Biochem. Sci., 38(2013), 210–218. 2. J. H. Ge and Y. T. Yu, RNA pseudouridylation: New insights into an old modification, Trends Biochem. Sci., 38(2013), 210–218.
    [3] M. Charette and M. W. Gray, Pseudouridine in RNA: what, where, how, and why, IUBMB Life, 49(2010), 341–351. 3. M. Charette and M. W. Gray, Pseudouridine in RNA: what, where, how, and why, IUBMB Life, 49(2010), 341–351.
    [4] D. R. Davis, C. A. Veltri, L. J. J. o. B. S. Nielsen, et al., An RNA model system for investigation of pseudouridine stabilization of the codon-anticodon interaction in tRNALys, tRNAHis and tRNATyr, J. Biomol. Struct. Dyn., 15(1998), 1121–1132.
    [5] A. Basak and C. Query, A pseudouridine residue in the spliceosome core is part of the filamentous growth program in yeast, Cell Reports, 8(2014), 966–973.
    [6] X. Yang, Y. Yang, B. F. Sun, et al., 5-methylcytosine promotes mRNA export-NSUN2 as the methyltransferase and ALYREF as an m5C reader, Cell Res., 27(2017), 606–625.
    [7] M. Frye and F. M. Watt, The RNA methyltransferase Misu (NSun2) mediates Myc-induced proliferation and is upregulated in tumors, Curr. Biol., 16(2006), 971–981.
    [8] X. Wang, Z. Lu, A. Gomez, et al., N6-methyladenosine-dependent regulation of messenger RNA stability, Nature, 505(2014), 117–120.
    [9] C. Roost, S. R. Lynch, P. J. Batista, et al., Structure and thermodynamics of N6-methyladenosine in RNA: A spring-loaded base modification, J. Am. Chem. Soc., 137(2015), 2107–2115.
    [10] T. Chen, Y. J. Hao, Y. Zhang, et al., m6A RNA methylation is regulated by micrornas and promotes reprogramming to pluripotency, Cell Stem Cell, 16(2015), 289–301.
    [11] S. Geula, S. Moshitch-Moshkovitz, D. Dominissini, et al., m6A mRNA methylation facilitates resolution of naive pluripotency toward differentiation, Science, 347(2015), 1002–1006.
    [12] X. Li, X. Xiong, K. Wang, et al., Transcriptome-wide mapping reveals reversible and dynamic N1-methyladenosine methylome, Nat. Chem. Biol., 12(2016), 311.
    [13] S. Nachtergaele and C. J. R. B. He, The emerging biology of RNA post-transcriptional modifications, RNA Biol., 14(2016), 156–163.
    [14] W. Chen, P. M. Feng, H. Tang, et al., RAMPred: Identifying the N-1-methyladenosine sites in eukaryotic transcriptomes, Sci. Rep., 6(2016), 31080.
    [15] W. Chen, H. Tang, J. Ye, et al., iRNA-PseU: Identifying RNA pseudouridine sites, Mol. Ther.-Nucl. Acids, 5(2016).
    [16] J. J. He, T. Fang, Z. Z. Zhang, et al., PseUI: Pseudouridine sites identification based on RNA sequence information, BMC Bioinform., 19(2018), 306.
    [17] W. R. Qiu, S. Y. Jiang, Z. C. Xu, et al., iRNAm5C-PseDNC: Identifying RNA 5-methylcytosine sites by incorporating physical-chemical properties into pseudo dinucleotide composition, Oncotarget, 8(2017), 41178–41188.
    [18] J. W. Li, Y. Huang, X. Y. Yang, et al., RNAm5Cfinder: A web-server for predicting RNA 5-methylcytosine (m5C) sites based on random forest, Sci. Rep., 8(2018).
    [19] P. M. Feng, H. Ding, H. Yang, et al., iRNA-PseColl: Identifying the occurrence sites of different RNA modifications by incorporating collective effects of nucleotides into PseKNC, Mol. Ther.-Nucl. Acids, 7(2017), 155–163.
    [20] W. Chen, P. M. Feng, H. Yang, et al., iRNA-3typeA: Identifying three types of modification at RNA's adenosine sites, Mol. Ther.-Nucl. Acids, 11(2018), 468–474.
    [21] Y. Huang, N. N. He, Y. Chen, et al., BERMP: A cross-species classifier for predicting m6A sites by integrating a deep learning algorithm and a random forest approach, Int. J. Biol. Sci., 14(2018), 1669–1677.
    [22] J. J. Xuan, W. J. Sun, P. H. Lin, et al., RMBase v2.0: Deciphering the map of RNA modifications from epitranscriptome sequencing data, Nucleic Acids Res., 46(2018), D327–D334.
    [23] D. Dominissini, S. Moshitch-Moshkovitz, S. Schwartz, et al., Topology of the human and mouse m6A RNA methylomes revealed by m6A-seq, Nature, 485(2012), U201–U284.
    [24] L. Fu, B. Niu, Z. Zhu, et al., CD-HIT: Accelerated for clustering the next-generation sequencing data, Bioinformatics, 28(2012), 3150–3152.
    [25] W. Z. Li and A. Godzik, Cd-hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences, Bioinformatics, 22(2006), 1658–1659.
    [26] L. Zhu, H. B. Zhang and D. S. J. B. Huang, Direct AUC optimization of regulatory motifs, Bioinformatics, 33(2017), i243.
    [27] H. Zhang, L. Zhu and D. S. J. S. R. Huang, WSMD: Weakly-supervised motif discovery in transcription factor ChIP-seq data, Sci. Rep., 7(2017).
    [28] G. H. Chuai, H. H. Ma, J. F. Yan, et al., DeepCRISPR: Optimized CRISPR guide RNA design by deep learning, Genome Biol., 19(2018).
    [29] Q. Zhang, L. Zhu and D. S. Huang, High-order convolutional neural network architecture for predicting DNA-protein binding sites, IEEE/ACM Transact. Comput. Biol. Bioinform., (2018), 1.
    [30] Q. Zhang, L. Zhu, W. Bao, et al., Weakly-supervised convolutional neural network architecture for predicting protein-DNA binding, IEEE/ACM Transact. Comput. Biol. Bioinform., (2018), 1.
    [31] A. Krizhevsky, I. Sutskever and G. E. Hinton, ImageNet classification with deep convolutional neural networks, NIPS. Curran Assoc. Inc., (2012).
    [32] D. P. Kingma and J. J. C. S. Ba, Adam: A method for stochastic optimization, (2014).
    [33] C. Tan, F. Sun, K. Tao, et al., A survey on deep transfer learning, (2018).
    [34] G. Litjens, T. Kooi, B. E. Bejnordi, et al., A survey on deep learning in medical image analysis, Med. Image Anal., 42(2017), 60–88.
    [35] S. Liang, R. G. Zhang, D. Y. Liang, et al., Multimodal 3D denseNet for IDH genotype prediction in gliomas, Genes, 9(2018).
    [36] L. Zhu, W. L. Guo, C. Lu, et al., Collaborative completion of transcription factor binding profiles via local sensitive unified embedding, IEEE Transact. NanoBiosci., (2016), 1.
    [37] J. X. Wang, L. Chen, Y. Wang, et al., A computational systems biology study for understanding salt tolerance mechanism in rice, Plos One, 8(2013), 177–194.
  • This article has been cited by:

    1. D. Flore, K. Wegener, H. Mayer, U. Karr, C.C. Oetting, Investigation of the high and very high cycle fatigue behaviour of continuous fibre reinforced plastics by conventional and ultrasonic fatigue testing, 2017, 141, 02663538, 130, 10.1016/j.compscitech.2017.01.018
    2. F. Delzendehrooy, M.R. Ayatollahi, A. Akhavan-Safar, L.F.M. da Silva, Strength improvement of adhesively bonded single lap joints with date palm fibers: Effect of type, size, treatment method and density of fibers, 2020, 188, 13598368, 107874, 10.1016/j.compositesb.2020.107874
    3. Joginder Singh, M. R. Tyagi, 2020, Chapter 38, 978-981-15-2646-6, 427, 10.1007/978-981-15-2647-3_38
    4. Wen Chen, Jianzhong Hao, Ming Tang, Improved estimate and accurate measurement of thermal stresses in FRP tendon, 2018, 164, 09500618, 620, 10.1016/j.conbuildmat.2017.12.151
    5. Ranjan K. Behera, S.K. Parida, R.R. Das, Effect of using fibre reinforced epoxy adhesive on the strength of the adhesively bonded Single Lap Joints, 2023, 248, 13598368, 110358, 10.1016/j.compositesb.2022.110358
  • Reader Comments
  • © 2019 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(6266) PDF downloads(928) Cited by(26)

Figures and Tables

Figures(4)  /  Tables(3)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog