1。介绍
结核病患者中占主导地位的人类免疫缺陷病毒(HIV)根据世界卫生组织(
1 )报告。因此,结核病是死亡率的主要原因个人容易受到艾滋病毒(
2 ]。因此,两种疾病的综合效应被认为是极其致命的个人或边际效应相比每个个体潜在的疾病。此外,(
3 合同)证实,大多数艾滋病患者肺结核作为第一个艾滋病的证据,随着全球大约33.3%的3860万名艾滋病毒阳性患者,也感染了肺结核,都是承包的风险完全结核病。其他可能的乙型肝炎病毒感染传播,占主导地位的模式的高传染性疾病传播的感染下的社区研究(库马西,阿散蒂地区)及其环境通常是通过输血感染,无保护的性行为,使用未经消毒的叶片或修整设备,在交付(母婴传播
4 ]。另一方面,艾滋病毒传播艾滋病毒阳性的个体与性有关的风险行为,因此把感染艾滋病毒感染的风险。
根据(
5 ),3600万感染艾滋病毒的人,大约有400万被发现有慢性乙肝病毒。
的免疫系统是人体的自然防御系统对病原体的感染,疾病是由一些特殊的细胞称为CD4细胞。这些底层细胞是白细胞,反对各种形式的人类感染系统,从而提供免疫系统的性能的概述。因此,CD4细胞计数的一个人在任何时候决定他或她的力量水平的免疫力。然而,CD4细胞计数的大多是受到病毒和细菌感染。这就解释了为什么这些潜在的任意两个的疾病(艾滋病、结核病、和HB)在个人是极其危险的
6 ]。因此,这些疾病在考虑必须被视为一个公共卫生问题在加纳和整个世界。成年人和青少年在加纳,一线艺术方案的规定由nonthymine核苷逆转录酶抑制剂(替诺福韦
+
emtricitabine或替诺福韦
+
拉米夫定)和一个nonnucleoside逆转录酶抑制剂通常依法韦伦[
7 ]。
提出了一些数学和标准的流行病学模型在研究传染病动力学。这些模型有助于我们理解疾病的传播感染。然而,大多数的这些模型不能估计其他相关疾病指标如第一次感染的概率和恢复以及预计感染和恢复时间敏感和感染者。大多数标准的流行病学模型用于估算过渡概率(TPs)不能概括过渡的估计疾病结果离散时间为未来的预测步骤。本文旨在解决上述问题采用离散时间马尔可夫链模型提出了(
8 ]。在这项研究中,
n th-step转移概率矩阵为每个疾病决心使用一阶马尔可夫链模型。这些发现将有助于决策者确定哪些疾病的研究最致命传染病和这样的措施可以减少或降低其发病率。
2。材料和方法
2.1。数据来源
研究使用的数据收集从一个地区性医院在加纳,因为它作为一个主要的转诊中心。数据获得定期跟进艾滋病毒、结核病、乙型肝炎患者从2016年1月至2016年12月。每个病人之后一年不管他们的条目次课程。结核病复发或再感染患者以及患者证实合并感染或其他医疗并发症(失访)被排除在外或撤回就医。对象之间的敏感(感染),这些疾病经过体检,项发现那些被感染和死亡(或免疫疾病可能)的研究期间(直到最后合格主题之后一年)被记录下来。类似的,在研究过程中受试者感染(医学筛查后每个疾病),项恢复(TB)的人,感染,和死也记录下来。
2.2。S-I-R建模通过马尔可夫链过程
2.2.1。模型开发
考虑三种离散状态:敏感(0)状态,感染(状态1)和删除(状态2)状态。如果(
X
我
,
我
=
0、1、2
)代表个体的数量在任何国家在任何时候从基础疾病
t ,那么很明显,
X
我
是一个随机过程,状态0、1和2。
因此,一阶time-homogeneous马尔可夫依赖统计可以建模为
(1)
P
X
n
=
我
n
X
n
−
1
=
我
n
−
1
,
…
,
X
1
=
我
1
,
X
0
=
我
0
=
P
X
n
=
我
n
X
n
−
1
=
我
n
−
1
。
然后,转移概率
P
我
j
为
我
,
j
=
0、1、2
用矩阵形式表示吗
(2)
P
我
j
=
P
00
P
01
P
02
P
10
P
11
P
12
0
0
1
,
在哪里
(3)
∑
j
=
0
2
P
我
j
=
1
,
我
=
0、1、2
。
2.2.2。州S-I-R模型的定义
易感状态(S):由个人没有暴露前和个人从感染中恢复过来。
感染状态(我):它由感染者和疾病的携带者。
删除状态(R):由个人死于疾病或免疫恢复后的疾病(艾滋病除外)在研究过程中。
2.2.3。参数马尔可夫链(过渡的可能性)
P
我
我
:剩余的状态的概率
我 。
P
我
j
:从状态转移概率
我 州
j ,
我
≠
j
。
言论。
的参数
P
01
主要是将文献中称为离散时间的感染力量。此外,元素
P
02
和
P
12
分别表示未受感染的死亡率和感染者
P
10
经济复苏或背叛概率(
9 ]。死亡是一个吸收态自成为敏感或感染的概率是零。时间步长单位,确保从一个状态转换到另一个原因是每年来衡量。
2.3。模型的假设
当前状态的个体仅依赖个体的状态在前一个时间步。
没有单独的删除状态可以感染或感染。
转移概率是独立于时间和随时间保持不变或研究期间。
连续转换或复发证实合并感染的疾病或其他医疗并发症没有考虑或不符合资格标准的研究。
删除状态由受试者死于疾病或发现免疫后恢复。
唯一的假设所需的有关损失和取款是他们未来有相同的经验,那些剩余的观察。
2.4。估计转换概率
最大似然估计(标定)被用来估计每个疾病的跃迁概率与各自的标准错误。表
1 在研究期间显示个人的数量在任何状态(S-I-R)为每个队列研究的三种疾病(艾滋病、结核病、和HB)。
表1
个人在任何状态队列的数量。
组
易受影响
受感染的
死
敏感的人
X
00
X
01
X
02
感染者
X
10
X
11
X
12
X
00
:易感个体数量在研究结束的时期保持敏感;
X
01
:被感染的易感个体人数在研究期间的结束;
X
02
:数量的易感个体死亡或免疫恢复后仍在研究期间的结束;
X
10
:恢复的感染者人数在研究期间的结束;
X
11
:最后的感染者人数仍然受感染的研究时期;
X
12
:许多被感染的个体死亡或免疫恢复后仍在研究结束的时期。
言论。
感染者的数量恢复第二阶段的末尾或研究期间的群组研究艾滋病和乙肝会零因为艾滋病毒和慢性乙型肝炎不治愈而不是肺结核,但可以治疗。存在一个有限的治疗,尤其是对很多情况下乙肝。
过渡事件是相互独立的(定义的马尔可夫原则);转移概率的可能性,
P
我
j
之前,二项模型:
(4)
l
P
我
j
N
,
x
=
N
我
x
我
j
P
我
j
x
我
j
1
−
P
我
j
N
我
−
x
我
j
,
在哪里
N
我
j
是观察到的转换,从状态的数量吗
我 来
j 和
(5)
∑
j
P
我
j
=
1。
从方程(
4 ),在此期间不断跃迁概率的假设,转移概率矩阵估计多项分布给出
(6)
P
^
我
j
=
x
我
j
∑
j
x
我
j
=
x
我
j
N
我
,
为
我
,
j
=
0 1
抽样分布的,与标准错误的ML估计
(7)
年代
^
。
e
P
我
j
=
p
^
我
j
1
−
p
^
我
j
N
我
。
2.5。估计疾病指标
易得病的概率之间的首次个人都会被感染
米
−
1
和
米 时间步的状态
我
,
j
=
0 1
从转移概率矩阵(S-I-R)给出
(8)
f
01
米
=
P
X
n
+
米
=
1
,
X
n
+
米
−
1
=
0
,
…
,
X
n
+
1
=
0
X
n
=
0
=
P
00
米
−
1
P
01
。
同样的,个人之间第一次恢复被感染的概率
米
−
1
和
米 给药时间步骤
(9)
f
10
米
=
P
X
n
+
米
=
0
,
X
n
+
米
−
1
=
1
,
…
,
X
n
+
1
=
1
X
n
=
1
=
P
11
米
−
1
P
10
。
预期的时间来感染和恢复有一个封闭的解,计算
(10)
E
τ
我
j
1
=
∑
米
=
1
∞
米
f
我
j
米
P
r
我
⟶
j
=
1
1
−
P
我
我
,
为
我
,
j
=
0 1
,
我
≠
j
分子,
∑
米
f
我
j
米
是首次通过时间的期望值
我 州
j 和分母
(11)
P
r
我
⟶
j
=
P
我
j
1
−
P
我
我
总体概率或终生的概率过渡状态
我 州
j (
8 ]。
预期寿命(
W
我
,
我
=
0 1
敏感和感染者估计也可以使用以下方程:
(12)
W
=
我
−
问
−
1
×
1
1
,
在哪里
我 是一个
2
×
2
单位矩阵,
(13)
问
=
P
00
P
01
P
10
P
11
。
2.6。估计< inline-formula > < mml:数学xmlns: mml = " http://www.w3.org/1998/Math/MathML " id = " M51 " > < mml: mrow > < mml: msup > < mml: mrow > < mml: mi > P < / mml: mi > < / mml: mrow > < mml: mrow > < mml: mi > n < / mml: mi > < / mml: mrow > < / mml: msup > < / mml: mrow > < / mml:数学> < / inline-formula >过渡矩阵
估算的方法
n
th
一步一步转移概率矩阵为每个疾病使用Bhat提出的特征值和特征向量的方法是(
10 ]。的
P
我
j
n
,
我
,
j
=
0、1、2
转移概率矩阵估计为每个疾病使用分解方法,需要特征值及其对应的特征向量。因此,它可以用下面的分解:
(14)
P
n
=
问
Λ
n
问
−
1
,
在哪里
问 是
3
×
3
满秩矩阵
X
0
,
X
1
,
X
2
和
X
j
,
j
=
0、1、2
特征值所对应特征向量是正确的吗
λ
j
j
=
0、1、2
。因此,
(15)
P
X
j
=
λ
j
X
j
,
(16)
Λ
n
=
λ
0
n
0
0
0
λ
1
n
0
0
0
λ
2
n
。
3所示。结果与讨论
3.1。马尔可夫链S-I-R模型
假设的离散状态的马尔可夫链模型的疾病(艾滋病、结核病、和HB)敏感(0)状态,感染(1),和删除或死亡(状态2)状态。让
X
我
,
我
=
0、1、2
代表个人的数量在任何国家在任何时候从基础疾病
t 从方程(满足一阶time-homogeneous马尔可夫的依赖
1 )。很明显,
X
我
满足马尔可夫链状态空间模型
年代
=
0、1、2
。
表
2 - - - - - -
4 展示个人在研究人口的数量状态0,1和2艾滋病,结核病和HB,分别。
表2
许多个人在任何国家艾滋病毒的最后时期。
组
总
易受影响的国家
感染状态
删除状态
易受影响
8132年
7598年
502年
23
受感染的
502年
0
478年
24
表3
许多个人在任何国家结核病的结尾段。
组
总
易受影响的国家
感染状态
删除状态
易受影响
669年
552年
99年
18
受感染的
99年
81年
18
0
表4
许多个人在任何国家HB的结尾段。
组
总
易受影响的国家
感染状态
删除状态
易受影响
2929年
2281年
623年
25
受感染的
623年
0
600年
23
3.2。估计每个疾病的转移概率
采用最大似然估计为每个疾病转移概率矩阵。表
5 - - - - - -
7 显示转换概率的估计和置信区间为艾滋病,结核病和HB。
表5
毫升的跃迁概率估计艾滋病毒。
参数
估计
SE
99% conf.区间
P
00
0.93
0.00283
0.923 - -0.937
P
01
0.06
0.00263
0.053 - -0.067
P
02
0.01
0.00110
0.007 - -0.013
P
11
0.95
0.00973
0.925 - -0.975
P
12
0.05
0.00973
0.025 - -0.075
表6
毫升结核病的跃迁概率的估计。
参数
估计
SE
99% conf.区间
P
00
0.83
0.01452
0.793 - -0.867
P
01
0.15
0.01381
0.114 - -0.186
P
02
0.02
0.00541
0.006 - -0.034
P
10
0.82
0.03861
0.721 - -0.919
P
11
0.18
0.03861
0.081 - -0.279
表7
ML估计HB的转移概率。
参数
估计
SE
99% conf.区间
P
00
0.78
0.00765
0.760 - -0.800
P
01
0.21
0.00752
0.191 - -0.229
P
02
0.01
0.00184
0.005 - -0.015
P
11
0.96
0.00785
0.940 - -0.980
P
12
0.04
0.00785
0.020 - -0.060
疾病的转移概率矩阵,分别介绍
3.3。州用图算法的分类模型
图算法被用来代表每个疾病的跃迁概率,容易把美国分成复发,瞬态或吸收。从图
1 ,可以得出结论:{0}和{1}都是瞬态状态的转换对艾滋病毒和肝炎b。这是因为,在一个人离开这些州,没有积极的返回其原始状态的概率。然而,国家{0,1}是一个瞬态类的转换对于肺结核因为两个国家有相同的等价类以及沟通,但是,在一个人离开这个类,没有积极的返回该类的概率。状态{2}也称为删除状态是任何疾病的吸收状态转换从一个人在这个国家仍死一定概率的一个。因此,可以推断出从国家分类,马尔可夫链不是遍历。尽管链周期,这不是不可约以来所有的州不属于同一等价类。图
1 显示一个图形表示的转换概率,分别为每一个疾病。
图1
导演油印的从一个状态过渡到另一个地方。
3.4。估计疾病指标
3.4.1。第一个转变的概率
易感个体成为第一个感染的概率估计的潜在疾病使用方程(
8 )。此外,受感染的个人第一的概率变得敏感或恢复估计只有肺结核使用方程(
9 ),因为艾滋病毒和乙肝病毒感染没有复苏的情况下。图
2 显示了一个情节首先过渡概率在不同时间的步骤从1到50。
图2
第一次感染的概率和恢复。
从图
2 ,它可以观察到三个潜在疾病的第一次感染的概率在研究人口相对较低(低于0.20)。这可能是由于干预(或药物)管理研究人口控制这些疾病的传播(
11 ]。此外,它也可以推断出,在第十次步骤(在10年之后,研究期间),首次乙肝感染的概率易感个体相对较高,其次是结核感染和艾滋病毒感染。相反,第十次后感染的第一步的概率相对较高在艾滋病患者与结核病和HB患者。这一发现是一致的,在获得
12 ),大多数人接触到艾滋病毒首先获得感染后10年的接触没有治疗。同时,第一个康复肺结核的概率是非常高的低次步骤但拒绝大幅增加时间步长。因此,它只是表明,如果一个病人诊断结核病的早期或潜伏阶段,持续控制措施可以实施第一次从疾病中恢复过来。这是因为大多数活动性结核病患者接受适当的治疗至少2周后不再是会传染的。然而,如果一个结核病患者感染的时间较长,还有一个相对很小的概率从疾病中恢复在未来时间的步骤。
3.5。在任何时间内感染的概率
感染的概率在任何时候从三个潜在疾病估计使用的总合第一感染易感个体的概率。图
3 是一个概率随时间累积的情节过渡。
图3
多年来为感染者感染的概率。
感染的概率的估计从第一过渡概率如图
3 发现乙肝似乎更多的传染性研究人口中随着时间的推移,其次是结核病和艾滋病毒感染,分别。这一发现同意从其他健康调查的结果。谁(
12 ]从他们多年来的研究还发现,尽管艾滋病毒和乙肝的传播方式是类似的,乙肝的传染性比艾滋病毒强50到100倍。这是因为乙肝病毒与HIV甚至可以住外主机至少一个星期,还导致易感个体感染以及病毒孵化率较高。此外,与乙型肝炎病毒相关血液水平相对高于艾滋病毒,从而导致这种病毒容易传播的情况下如孩子交付(母婴),甚至通过体液,不能在艾滋病毒传播的情况发生。另一方面,乙肝还发现以来比肺结核传染乙肝感染的概率相对高于任何时候结核病传播。细菌感染,如结核病主要由抗生素治愈这对病毒感染无效,他们感染动力学甚至可以控制与乙型肝炎在任何阶段的风险。然而,艾滋病似乎感染后第一次感染相对较少,这可能是由于个人的病毒首先攻击免疫系统通常需要大约5到10年它完全妥协,即使不及时治疗(
12 ]。因此,个人与很强的免疫系统发展感染以相对较低的速度与结核病和HB感染也会影响免疫系统。感染的概率也透露,这三种疾病的增加而增加。这进一步表明,艾滋病、结核病和HB传染性如果放任不管增加时间。这占这三种疾病的致命感染由于其高速率。这也证实了为什么这些疾病的综合效应的任意两个个体可以是极其致命的比个人影响的人患有艾滋病,结核病,或者HB。
3.6。其他估计疾病指标
3.6.1。总体的概率感染和恢复
其他疾病感染和恢复的总体概率等指标估计使用方程(
11 )。的总体概率感染艾滋病、结核病、HB被发现大约0.86,0.88,和0.95,分别。这些发现符合结果的累积金额首次感染概率图
3 随着时间的推移,HB是传染性结核病和艾滋病。因此,可以推断,乙型肝炎是最传染病在研究期间在地区医院。研究还发现,某些治疗结核与艾滋病毒和乙肝无法治愈的感染或慢性阶段。
操作。评估的总体概率的统计显著性感染
估计总体的统计显著性概率感染每个疾病的马尔可夫模型测量使用引导区间估计技术在1%水平的意义。因此,它可以得出有99%的信心水平的总体概率感染每个疾病会落在置信区间估计总结如表
8 。
表8
估计引导置信区间。
总体的概率感染
马尔可夫链预测
SE
99% conf.区间
艾滋病毒
0.86
0.00385
0.84981 - -0.86977
肺结核
0.88
0.01218
0.84753 - -0.91226
乙型肝炎
0.95
0.00397
0.93979 - -0.96005
3.7。预计感染和恢复时间
此外,预计感染和恢复时间估计方程(
10 每个疾病)。预期时间感染易感个体被发现14.29年,5.88年和4.55年艾滋病,肺结核,乙肝感染,群内分别考虑。即个体暴露于艾滋病毒将开始经历感染后平均约14.29年即使的时间长度可以是个体之间的差异很大。这一发现肯定,获得在
12 ),个人暴露于艾滋病毒充分发展感染的症状和体征通常10至15年。另一方面,个体暴露于结核病的平均时间将成为完全感染暴露后如果不及时治疗被估计为5.88年,而HB被估计为4.55年。因此,可以推断,乙型肝炎感染平均发展速度比艾滋病毒和结核病感染。
相反,结核病被发现有一个预期恢复时间1.22年,而预期的时间在艾滋病毒和乙肝的症状和体征是减少估计20年和25年,分别。然而,将恢复从艾滋病毒和乙型肝炎为本研究定义正好与各自的预期寿命的估计感染者使用马尔可夫链模型。可以得出结论,因此,它的寿命受感染的个人关于艾滋病毒和乙肝确实超出了受感染的个人复苏时期,要么死亡。
3.8。健康的预期寿命和感染者
的寿命,这是一个统计测量的平均时间预计个人生活或生存,估计使用方程(
12 )对健康和受感染的个人。因此,健康的人的预期寿命在HB、艾滋病毒、结核病是估计为28.41年,31.43年和59.15年。这些估计进一步透露,在人口乙型肝炎流行,任何健康个体的平均寿命相对较小而获得更高的估计人口感染艾滋病毒和结核病,分别。然而,平均寿命为个人感染了艾滋病毒,乙肝,和结核病被发现20年,25年,分别为60.37年。这意味着艾滋病患者的平均寿命相对低于HB和结核病患者即使他们被发现是非常比艾滋病毒感染。这是因为当艾滋病患者的免疫系统被破坏随着时间的推移,它创造了一种利用削弱免疫系统的各种感染结核病和其他并发症等影响的受害者。表
9 介绍了其他疾病指标的估计(感染的概率,预计感染/时间复苏,预期寿命和健康/感染者)艾滋病,肺结核,乙肝。
表9
估计每个疾病的其他疾病的指标。
指标
艾滋病毒
肺结核
乙型肝炎
总体的概率感染
0.86
0.88
0.95
整体复苏的概率
- - - - - -
1.00
- - - - - -
预计感染时间(年)
14.29
5.88
4.55
预期恢复时间(年)
20.00
∗
1.22
25.00
∗
健康的人的预期寿命(年)
31.43
59.15
28.41
感染者的预期寿命(年)
20.00
60.37
25.00
∗
预计恢复时间从艾滋病和结核病感染者也称为寿命。
3.9。估计< inline-formula > < mml:数学xmlns: mml = " http://www.w3.org/1998/Math/MathML " id = " M81 " > < mml: mrow > < mml: msup > < mml: mrow > < mml: mi > P < / mml: mi > < / mml: mrow > < mml: mrow > < mml: mi > n < / mml: mi > < / mml: mrow > < / mml: msup > < / mml: mrow > < / mml:数学> < / inline-formula >转移概率矩阵
的
P
n
转移概率矩阵估计方程(
14 )。链是周期但不是不可约自删除状态(2)是一个吸收状态。的
P
n
转移概率矩阵预测每个疾病在任何时间步转移概率。
P
n
转移概率矩阵为艾滋病毒估计
(18)
P
H
我
V
n
=
0.93
n
3
0.95
n
−
0.93
n
1
−
3
0.95
n
+
2
0.93
n
0
0.95
n
1
−
0.95
n
0
0
1
。
很明显,
n
=
1
的估计
P
n
转换矩阵给出了实际先转换矩阵
P
H
我
V
。因此,在任何时间步转移矩阵
n
≥
1
可以生成的安装吗
P
n
矩阵。
另外,估计
n 一步一步转移矩阵给出了结核病
(19)
P
结核病
n
=
0.8
0.9831
n
+
1.6
0.0269
n
0.16
0.9831
n
−
0.0269
n
1
−
0.9831
n
−
0.0032
0.0269
n
0.86
0.9831
n
−
0.0269
n
0.16
0.9831
n
+
0.84
0.0269
n
0
0
0
1
。
很明显,
n
=
1
的估计
P
n
转移矩阵(2一般。)给出了实际第一过渡矩阵
P
T
B
。
最后,
n 一步一步转移概率给出了乙肝
(20)
P
乙肝
n
=
0.78
n
1.17
0.96
n
−
0.78
n
1
−
1.17
0.96
n
+
0.17
0.78
n
0
0.96
n
1
−
0.96
n
0
0
1
。
同时,考虑到
n
=
1
的估计
P
n
(在2一般。)转移矩阵给出了第一个过渡矩阵
P
乙肝
。