Ushbu tadqiqot uchun ishlab chiqilgan Ebola virusi kasalligi yuqish modeli, inson populyatsiyasini sezgir \(S\), ta’sirlangan \(E\), yuqumli \(I\), kasalxonaga yotqizilgan \(H\), sog‘aygan \(R\), vafot etgan, ammo hali dafn qilinmagan \(D\) va dafn qilingan \(B\) bo‘limlariga ajratadigan; atrof-muhitdagi Ebola virusi yuklamasini esa \(P\) o‘zgaruvchisi bilan ifodalovchi sakkiz holat o‘zgaruvchili, deterministik va chiziqli bo‘lmagan oddiy differensial tenglamalar (ODE) modelidir. Model yuqumli shaxslar, kasalxonadagi bemorlar, o‘lik jasadlar, dafn qilingan jasadlar va kontaminatsiyalangan atrof-muhit orqali yuqishni bitta matematik tuzilma ichida birlashtiradi. Tadqiqot shuningdek kasalliksiz va endemik muvozanatlarni, asosiy reproduksiya sonini \(\mathcal R_0\), mahalliy va global barqarorlikni, bifurkatsiyani, parametr sezgirligini, DRC kuzatuv ma’lumotlari bilan parametr moslashtirishni va uchta optimal boshqaruv o‘zgaruvchisini o‘rganadi.
Mualliflarning asosiy ilmiy maqsadi Ebola yuqishida insondan insonga kontakt bilan bir qatorda, o‘limdan keyingi kontakt, dafn jarayonlari, uy sharoitidagi yaqin parvarish va atrof-muhitdagi virus yuklamasi kabi yo‘llarning nisbiy rolini miqdoriy model doirasida ko‘rsatishdir. Biroq tadqiqotdagi “African tradition practices” iborasi barcha Afrika madaniyatlarini bitta xulq-atvor qolipi sifatida ifodalamasligi kerak. Matematik model ayniqsa infeksiyalangan shaxsga yaqin parvarish, o‘lik jasad bilan kontakt, dafn jarayoni va bu jarayonlarda atrof-muhitga patogen sochilishi kabi muayyan xulq-atvor yo‘llarini soddalashtirilgan parametrlar bilan ifodalaydi; etnografik dala tadqiqoti yoki madaniy xulq-atvor prevalensiyasini o‘lchash amalga oshirilmagan.
Modelning sezgirlik jadvalida eng katta musbat normallashtirilgan sezgirlik indekslaridan biri atrof-muhit yuqish koeffitsiyenti \(\beta_p\) uchun \(+7.878\times10^{-1}\), kasalxonada o‘limga o‘tish koeffitsiyenti \(\eta\) uchun \(+5.032\times10^{-1}\), kasalxonadan atrof-muhitga patogen sochilishi koeffitsiyenti \(\eta_1\) uchun \(+5.031\times10^{-1}\) va yuqumli shaxslardan atrof-muhitga sochilish koeffitsiyenti \(\sigma\) uchun \(+2.707\times10^{-1}\) deb xabar qilingan. Bunga qarshi atrof-muhit patogenining parchalanish koeffitsiyenti \(\alpha\), kasalxonada sog‘ayish koeffitsiyenti \(\xi\), dafn qilingan sinfdan chiqish koeffitsiyenti \(b\), o‘limdan dafn etishga o‘tish koeffitsiyenti \(d\) va kasalxonaga yotqizish koeffitsiyenti \(\gamma_1\) manfiy sezgirlik ko‘rsatadi.
Sonli natijalarda mualliflar \(u_1\) deb atagan “izolyatsiya” strategiyasi yakka qo‘llanilgan boshqaruvlar orasida eng kuchli kamaytiruvchi ta’sirga ega deb talqin qilingan. \(u_1+u_2\) izolyatsiya va davolash kombinatsiyasi juft strategiyalar orasida eng kuchli variant sifatida taqdim etilgan; \(u_1+u_2+u_3\) kombinatsiyasi esa yuqumli shaxslar, kasalxonaga yotqizilgan shaxslar va atrof-muhitdagi patogen yuklamasini birgalikda eng ko‘p bostiruvchi ssenariy sifatida xabar qilingan. Manba grafiklarida bu taqqoslash taxminan 20 kunlik simulyatsiya oynasida ko‘rsatilgan.
Bu natijalar real jamiyatda uchta aralashuvning o‘lchangan dala samaradorligi emas. Natijalar mualliflarning matematik modeliga, tanlangan parametr qiymatlariga, boshlang‘ich shartlarga, parametrlarni moslashtirish usuliga va boshqaruv xarajatlarining \(A_1=A_2=A_3=c_1=c_2=c_3=1\) deb olinishi bilan bog‘liq simulyatsiya natijalaridir. Model bir jinsli aralashishni faraz qiladi; makoniy geterogenlikni, inson harakatchanligini va an’anaviy xulq-atvorlarning vaqt davomida o‘zgarishini ifodalamaydi.
Tadqiqot muammosi va adabiyotdagi bo‘shliq
Avvalgi Ebola matematik modellarida atrof-muhit orqali yuqish, kasalxonaga yotqizish, karantin, qo‘rquv xulq-atvori, davolash quvvati yoki fraksional hosilali tuzilmalar o‘rganilgan. Ushbu tadqiqot mualliflari dafn va yaqin parvarish amaliyotlari ko‘plab modellarda bevosita o‘lchanadigan alohida yuqish mexanizmlari sifatida emas, balki umumiyroq epidemiologik bo‘limlar ichida bilvosita ifodalanganini ta’kidlaydi.
Tadqiqot bu bo‘shliqni vafot etgan va dafn qilingan shaxslarni alohida holat o‘zgaruvchilari sifatida saqlash va yuqumli, kasalxonadagi, vafot etgan va dafn qilingan shaxslardan atrof-muhitga virus sochilishini alohida koeffitsiyentlar bilan belgilash orqali ko‘rib chiqadi.
Bu yerda “madaniy xulq-atvorni miqdorlashtirish” iborasini ehtiyotkorlik bilan o‘qish kerak. Masalan, dafn qilingan jasadlardan atrof-muhitga sochilish koeffitsiyenti \(b_1\) va vafot etgan shaxslardan sochilish koeffitsiyenti \(d_1\) madaniy kontakt bilan bog‘liq yo‘llarning matematik ifodasiga kiradi; ammo 2-jadvalda \(b_1\) va \(d_1\) bevosita madaniy dala ma’lumotlaridan baholanmagan, assumed, ya’ni faraz qilingan parametrlar sifatida keltirilgan.
Bu Model Ebola Yuqishini Qaysi Yo‘llar Orqali Ifodalaydi?
Model beshta yuqish manbasini bitta infeksiya kuchida birlashtiradi: yuqumli odamlar \(I\), kasalxonaga yotqizilgan odamlar \(H\), hali dafn qilinmagan vafot etganlar \(D\), dafn qilingan shaxslar \(B\) va atrof-muhit patogen yuklamasi \(P\); bu manbalarning sezgir shaxslarni \(S\) ta’sirlangan \(E\) sinfiga o‘tkazish ta’siri mos ravishda \(\beta,\beta_1,\beta_2,\beta_3,\beta_p\) koeffitsiyentlari bilan og‘irlashtiriladi.
Modelning holat o‘zgaruvchilari
| Belgi | Manbadagi ma’nosi | Modeldagi roli |
|---|---|---|
| \(S\) | Susceptible — sezgir | Hali infeksiyalanmagan, infeksiyaga moyil odamlar. |
| \(E\) | Exposed — ta’sirlangan | Infeksiyalangan, ammo hali yuqumli sinfga o‘tmagan shaxslar. |
| \(I\) | Infectious — yuqumli | Jamiyatda yuqishni yuzaga keltira oladigan yuqumli shaxslar. |
| \(H\) | Hospitalized — kasalxonaga yotqizilgan | Sog‘liqni saqlash muassasasidagi yuqumli shaxslar. |
| \(R\) | Recovered — sog‘aygan | Kasalxonadan sog‘ayib chiqqan shaxslar. |
| \(D\) | Deceased — vafot etgan | EVD sababli vafot etgan, ammo hali \(B\) sinfiga o‘tmagan shaxslar. |
| \(B\) | Buried — dafn qilingan | Dafn qilingan shaxslar; model bu bosqichda ham atrof-muhitga sochilish davom etishi mumkinligini faraz qiladi. |
| \(P\) | Environmental pathogen load | Kontaminatsiyalangan sirt, suyuqlik yoki atrof-muhitdagi Ebola virusi/patogen yuklamasi. |
Infeksiya kuchi \(\Psi\)
Tadqiqotning asosiy yuqish tenglamasi:
\[ \Psi=\beta I+\beta_1H+\beta_2D+\beta_3B+\beta_pP \tag{1} \]
Bu formula nimani hisoblaydi? \(\Psi\), sezgir shaxsga modeldagi barcha yuqish manbalaridan keladigan umumiy infeksiya bosimini ifodalaydi. Yangi ta’sirlanish oqimi keyin \(\Psi S\) ko‘rinishida ODE tizimiga kiradi.
O‘zgaruvchilar: \(\beta\) yuqumli shaxslardan; \(\beta_1\) kasalxonadagi shaxslardan; \(\beta_2\) vafot etgan shaxslardan; \(\beta_3\) dafn qilingan shaxslardan va \(\beta_p\) atrof-muhit patogen rezervuaridan yuqishni og‘irlashtiruvchi koeffitsiyentlardir.
Birliklar: Manba koeffitsiyentlarning fizik birliklarini alohida ko‘rsatmaydi. Ammo \(\Psi S\) kishi/vaqt birligida yangi ta’sirlanish hosil qilishi kerakligi sababli \(\Psi\) ning o‘lchami matematik jihatdan \(1/\text{zaman}\) bo‘lishi kerak. Odamlar sonidan bevosita foydalanadigan massa-ta’sir shaklida \(\beta,\beta_1,\beta_2,\beta_3\) koeffitsiyentlarining birliklari odamlar soni masshtabiga, \(\beta_p\) ning birligi esa \(P\) qanday o‘lchanishiga bog‘liq.
Sakkizta chiziqli bo‘lmagan ODE
Modelning asosiy tizimi:
\[ \begin{aligned} \dot S &= \Delta-\Psi S-\mu S,\\ \dot E &= \Psi S-(\mu+\delta)E,\\ \dot I &= \delta E-(\gamma+\gamma_1+\mu)I,\\ \dot H &= \gamma_1I-(\eta+\xi+\mu)H,\\ \dot R &= \xi H-\mu R,\\ \dot D &= \gamma I+\eta H-dD,\\ \dot B &= dD-bB,\\ \dot P &= \phi_p+\sigma I+\eta_1H+d_1D+b_1B-\alpha P. \end{aligned} \tag{2} \]
\(\dot S\): sezgir aholi \(\Delta\) orqali tug‘ilish/migratsiya bilan tizimga kiradi; \(\Psi S\) orqali infeksiyaga ta’sirlanadi va \(\mu S\) orqali tabiiy o‘lim sababli chiqadi.
\(\dot E\): yangi ta’sirlanishlar \(\Psi S\); ta’sirlangan shaxslar \(\delta\) tezlikda yuqumli sinfga o‘tadi va \(\mu\) tabiiy o‘limiga duchor bo‘ladi.
\(\dot I\): ta’sirlangan shaxslar \(\delta E\) bilan yuqumli bo‘ladi; \(\gamma\) kasallik bilan bog‘liq o‘lim, \(\gamma_1\) kasalxonaga yotqizish va \(\mu\) tabiiy o‘lim chiqishlaridir.
\(\dot H\): kasalxonaga yotqizish oqimi \(\gamma_1I\). Kasalxonadagi shaxslar \(\eta\) bilan EVD sababli vafot etishi, \(\xi\) bilan sog‘ayishi yoki \(\mu\) bilan tabiiy o‘limga uchrashi mumkin.
\(\dot R\): sog‘ayish faqat kasalxona sinfidan \(\xi H\) oqimi bilan yuz beradi; model sog‘aygan shaxslar \(\mu\) tabiiy o‘limiga duchor bo‘lishini faraz qiladi.
\(\dot D\): jamiyatdagi yuqumli shaxslardan \(\gamma I\) va kasalxonadagilardan \(\eta H\) o‘lim oqimi keladi; \(dD\) bilan dafn qilingan sinfga o‘tiladi.
\(\dot B\): \(dD\) bilan dafn sinfiga kirish bor, \(bB\) ushbu sinfdan chiqish hadidir. Manba \(b\) ning fizik talqinini batafsil tushuntirmaydi.
\(\dot P\): atrof-muhitga tashqi/asosiy patogen kirishi \(\phi_p\), yuqumli shaxslardan \(\sigma I\), kasalxonadagilardan \(\eta_1H\), vafot etganlardan \(d_1D\) va dafn qilinganlardan \(b_1B\) bilan patogen kirishi yuz beradi; \(\alpha P\) patogenning atrof-muhitda parchalanishi/tozalanishini ifodalaydi.
Model nega musbatlik va chegaralanganlik isbotini amalga oshiradi?
Epidemiologik modelda odamlar soni yoki patogen zichligi manfiy bo‘lishi mumkin emas. Shu sababli mualliflar nol yoki musbat boshlang‘ich shartlardan boshlanadigan yechimlar barcha \(t>0\) davomida manfiy bo‘lmasligini ko‘rsatuvchi musbatlik teoremasini tuzgan.
Faol inson populyatsiyasi:
\[ N_A=S+E+I+H+R \]
deb ta’riflangan va Grönwall tengsizligi yordamida:
\[ N_A(t)\leq \frac{\Delta}{\mu} \tag{5} \]
chegarasi olingan. Xuddi shunday \(D\), \(B\) va \(P\) uchun yuqori chegaralar chiqarilgan. Uning matematik roli yechimning chekli biologik sohada qolishini kafolatlashdir.
Asosiy Reproduksiya Soni \(\mathcal R_0\) Ushbu Modelda Nimani Anglatadi?
\(\mathcal R_0\), to‘liq sezgir populyatsiyaga kirgan tipik infeksiya modeldagi inson, kasalxona, o‘limdan keyingi, dafn qilingan jasad va atrof-muhit yuqish yo‘llarining barchasi orqali o‘rtacha nechta yangi infeksiya hosil qila olishini umumlashtiradigan o‘lchovsiz chegara kattaligidir; biroq ushbu tadqiqotda klassik kasalliksiz muvozanat va \(\mathcal R_0\leq1\) barqarorlik talqini tashqi atrof-muhit patogen kirishi \(\phi_p=0\) bo‘lgan shart bilan bog‘liq.
Kasalliksiz muvozanat
Infeksiya mavjud bo‘lmaganda manba kasalliksiz muvozanatni:
\[ \mathrm{DFE}= \left( \frac{\Delta}{\mu}, 0,0,0,0,0,0,0 \right) \tag{8} \]
ko‘rinishida beradi.
Bu muvozanat nuqtasida \(E,I,H,R,D,B,P=0\). Ammo atrof-muhitga mustaqil patogen kirishi \(\phi_p>0\) bo‘lsa, \(\dot P\) tenglamasida nol atrof-muhit patogen yuklamasi saqlanmaydi, shu sababli mualliflar keyingi barqarorlik tahlilida \(\phi_p=0\) shartidan ochiq foydalanadi.
\(\mathcal R_0\) ni chiqarish
Mualliflar next-generation matrix yondashuvida infeksiyalangan tizim o‘zgaruvchilarini \(E,I,H,D,B,P\) sifatida olgan. \(K=FV^{-1}\) matritsasi tuzilgan va:
\[ \mathcal R_0=\rho(K) \]
deb ta’riflangan. Bu yerda \(\rho(K)\), \(K\) matritsasining spektral radiusi, ya’ni mutlaq qiymati eng katta xos qiymatidir.
Manba:
\[ a_1=\delta+\mu,\qquad a_2=\gamma+\gamma_1+\mu,\qquad a_3=\eta+\xi+\mu \]
ta’riflaridan keyin \(\mathcal R_0\) ni ikki asosiy yuqish zanjiri yig‘indisi ko‘rinishida beradi:
\[ \begin{aligned} \mathcal R_0={}& \left[ \beta+ \frac{\beta_2\gamma}{d}+ \frac{\beta_3\gamma}{b}+ \frac{\beta_p\sigma}{\alpha}+ \frac{\beta_p\gamma}{\alpha} \left(\frac{d_1}{d}+\frac{b_1}{b}\right) \right] \frac{\Delta\delta}{\mu a_1a_2}\\ &+ \left[ \beta_1+ \frac{\beta_2\eta}{d}+ \frac{\beta_3\eta}{b}+ \frac{\beta_p\eta_1}{\alpha}+ \frac{\beta_p\eta}{\alpha} \left(\frac{d_1}{d}+\frac{b_1}{b}\right) \right] \frac{\Delta\delta\gamma_1}{\mu a_1a_2a_3}. \end{aligned} \tag{10} \]
Birinchi katta had, yuqumli \(I\) holatidan boshlanadigan bevosita, o‘limdan keyingi va atrof-muhit yuqish yo‘llarining hissalarini jamlaydi.
Ikkinchi katta had, kasalxonaga yotqizish \(\gamma_1\) orqali davom etuvchi yuqish, kasalxona o‘limi va atrof-muhitga sochilish yo‘llarini qo‘shadi.
Matematik talqin: Musbat parametr \(\mathcal R_0\) ni oshirsa, epidemiyaning davom etish salohiyati ortadi; ayrim o‘tish/tozalanish parametrlarining oshishi yuqumli davrni yoki atrof-muhitdagi saqlanishni qisqartirgani uchun \(\mathcal R_0\) ni kamaytirishi mumkin.
Mahalliy va global barqarorlik
Manba \(\phi_p=0\) shartida kasalliksiz muvozanat nuqtasi \(\mathcal R_0\leq1\) bo‘lganda mahalliy asimptotik barqaror ekanini bildiradi.
Global tahlil uchun chiziqli hadlar bilan Volterra turidagi Lyapunov funksiyasi birlashtirilgan. Mos \(g_1,\ldots,g_5\) koeffitsiyentlari tanlangach hosila:
\[ \dot V_0= -\frac{\mu}{S}(S-S_0)^2 +a_1(\mathcal R_0-1)E \tag{14} \]
ko‘rinishiga keltirilgan.
\(\mathcal R_0\leq1\) bo‘lsa ikkinchi had musbat emas va birinchi had ham noldan katta bo‘la olmaydi. LaSalle invariantlik prinsipi yordamida kasalliksiz holatning global asimptotik barqarorligi asoslanadi.
Endemik muvozanat nimani ko‘rsatadi?
Mualliflar tizimning o‘ng tomonlarini nolga tenglashtirib endemik muvozanatni chiqargan. Dastlabki komponentlar:
\[ S^*=\frac{\Delta}{\mu\mathcal R_0} \]
\[ E^*=\frac{\Delta(\mathcal R_0-1)}{a_1\mathcal R_0} \]
\[ I^*=\frac{\Delta\delta(\mathcal R_0-1)} {a_1a_2\mathcal R_0} \]
\[ H^*=\frac{\Delta\delta\gamma_1(\mathcal R_0-1)} {a_1a_2a_3\mathcal R_0} \]
ko‘rinishida. \(\mathcal R_0-1\) omili tashqi atrof-muhit kirishi bo‘lmagan standart chegara holatida \(\mathcal R_0>1\) bo‘lganda infeksiyalangan muvozanat komponentlarining musbat bo‘lishini matematik tarzda ko‘rsatadi.
Bifurkatsiya tahlili
Manba markaziy ko‘pxillik nazariyasidan foydalanib \(\mathcal R_0=1\) atrofida bifurkatsiya yo‘nalishini o‘rganadi. Hisoblangan markaziy-ko‘pxillik koeffitsiyentlari:
\[ \mathbf a<0,\qquad \mathbf b>0 \]
bo‘lgani sababli mualliflar \(\phi_p=0\) holatida oldinga (forward) bifurkatsiya xulosasiga keladi. Oldinga bifurkatsiya klassik holatda kasalliksiz muvozanat \(\mathcal R_0<1\) tomonida barqaror qolishi va endemik tarmoq chegara oshirilganda paydo bo‘lishini anglatadi.
Biroq manbada bu yerda ifoda nomuvofiqligi bor. 3-shakl va matn \(\phi_p=0\) uchun oldinga bifurkatsiyani muhokama qilgan holda, shu paragraf oxirida “\(\phi_p\geq0\) uchun bifurkatsiya yo‘q” degan jumla mavjud. Bundan tashqari, mualliflarning avvalgi matematik izohiga ko‘ra \(\phi_p>0\) bo‘lganda nol atrof-muhit patogenli kasalliksiz muvozanat allaqachon mavjud emas. Verianla bu jumlani jimlik bilan tuzatib, boshqa tengsizlik yaratmaydi.
Normallashtirilgan sezgirlik indeksi
Parametrning \(\mathcal R_0\) ga nisbiy ta’siri:
\[ S_p^{\mathcal R_0} = \frac{\partial \mathcal R_0}{\partial p} \frac{p}{\mathcal R_0} \tag{30} \]
bilan hisoblangan.
Musbat indeks: parametr oshganda, boshqalari o‘zgarmasa \(\mathcal R_0\) mahalliy ravishda oshishga moyil.
Manfiy indeks: parametr oshganda \(\mathcal R_0\) mahalliy ravishda kamayadi.
Mutlaq kattalik: parametrdagi nisbiy o‘zgarishga \(\mathcal R_0\) qanchalik sezgirligini ko‘rsatadi. Bu mahalliy sezgirlik tahlilidir; parametrlarning global noaniqlik taqsimotlarini birgalikda namunaga oladigan global sensitivity tahlili emas.
| Parametr | Sezgirlik indeksi | Belgi | Manba tenglamasidagi asosiy roli |
|---|---|---|---|
| \(\beta_p\) | \(+7.878\times10^{-1}\) | Musbat | Atrof-muhit patogen yuklamasining yuqish koeffitsiyenti |
| \(\eta\) | \(+5.032\times10^{-1}\) | Musbat | Kasalxonadagi EVD o‘limlari |
| \(\eta_1\) | \(+5.031\times10^{-1}\) | Musbat | Kasalxonadan atrof-muhitga patogen sochilishi |
| \(\sigma\) | \(+2.707\times10^{-1}\) | Musbat | Yuqumli shaxslardan atrof-muhitga sochilish |
| \(\beta_1\) | \(+9.994\times10^{-2}\) | Musbat | Eq. (1) ga ko‘ra kasalxona manbali yuqish |
| \(\beta\) | \(+6.914\times10^{-2}\) | Musbat | Eq. (1) ga ko‘ra yuqumli inson kontakti |
| \(\beta_2\) | \(+4.320\times10^{-2}\) | Musbat | Vafot etgan shaxslardan yuqish |
| \(\gamma\) | \(+2.129\times10^{-2}\) | Musbat | Yuqumli sinfdan EVD o‘limi |
| \(b_1\) | \(+5.639\times10^{-3}\) | Musbat | Dafn qilingan shaxslardan atrof-muhitga sochilish |
| \(\beta_3\) | \(+4.547\times10^{-4}\) | Musbat | Dafn qilingan shaxslardan bevosita yuqish |
| \(\delta\) | \(+3.445\times10^{-4}\) | Musbat | Ta’sirlangandan yuqumli holatga o‘tish |
| \(\gamma_1\) | \(-3.608\times10^{-2}\) | Manfiy | Kasalxonaga yotqizish |
| \(d\) | \(-5.144\times10^{-2}\) | Manfiy | Vafot etgan sinfdan dafn qilingan sinfga o‘tish |
| \(b\) | \(-6.093\times10^{-2}\) | Manfiy | Dafn qilingan sinfdan chiqish |
| \(\xi\) | \(-6.233\times10^{-2}\) | Manfiy | Kasalxonada sog‘ayish |
| \(\alpha\) | \(-7.878\times10^{-2}\) | Manfiy | Atrof-muhit patogenining parchalanishi |
Modelda \(d_1\) parametri mavjud bo‘lishiga qaramay, manba 1-jadvalining sezgirlik ro‘yxatida \(d_1\) uchun alohida indeks xabar qilinmagan. Shu sababli Verianla \(d_1\) uchun sezgirlik qiymatini yaratmaydi.
Sezgirlik natijalarini talqin qilishdagi manba ichidagi nomlash muammosi
Manba model ta’rifida \(\beta\) yuqumli shaxs, \(\beta_1\) kasalxona, \(\beta_2\) vafot etgan shaxs, \(\beta_3\) dafn qilingan shaxs va \(\beta_p\) atrof-muhit yuqish koeffitsiyentidir. Bunga qarshi “The effect of contact rates” bo‘limida mualliflar yuqumli shaxslar uchun \(\beta_1\), kasalxonadagilar uchun \(\beta_2\), vafot etganlar uchun \(\beta_3\) ifodalaridan foydalanadi.
Shu sababli, masalan, natijalar bo‘limidagi “\(\beta_1\) 0 dan 0,0007 gacha ko‘tarilganda \(\mathcal R_0\) 0,62 dan 6,40 gacha oshdi” jumlasini Eq. (1) dagi qaysi kontakt koeffitsiyentiga moslashtirish kerakligi manba ichida to‘liq izchil emas. Verianla bu egri chiziqni o‘zicha \(\beta\) yoki \(\beta_1\) deb qayta belgilamaydi.
Xuddi shuningdek \(\sigma\), Eq. (2) da yuqumli shaxslardan atrof-muhitga patogen sochilishi koeffitsiyentidir; manbaning sezgirlik muhokamasidagi bir jumlada esa bu koeffitsiyent bevosita simptomatik inson kontakti bilan bog‘langan. Bu ikki mexanizm bir xil emas va manba ichidagi talqin farqi saqlanishi kerak.
Atrof-muhit Yuqishi O‘zi Yolg‘iz Epidemiyani Davom Ettiradimi?
Manba simulyatsiyalarida o‘rganilgan \(\beta_p\) oralig‘ida atrof-muhit yo‘li o‘zi yolg‘iz \(\mathcal R_0\) ni 1 dan yuqoriga chiqarmagan va mualliflar infeksiya uchun yuqori atrof-muhit patogen yuklamasi kerakligini bildirgan; biroq bu natija faqat modelning tanlangan parametr oralig‘iga tegishli va real dunyoda atrof-muhit orqali yuqish hech qanday sharoitda epidemiyani davom ettira olmasligini isbotlamaydi.
Manbada \(\beta_p\) 0 bilan 0,0007 oralig‘ida o‘zgartirilganda \(\mathcal R_0<1\) bo‘lib qolgan. Atrof-muhitga tashqi patogen kirishi \(\phi_p\) oshirilganda yuqumli shaxslar soni biroz ko‘paygan; mualliflar atrof-muhit yo‘li bevosita insondan insonga va o‘limdan keyingi kontaktga nisbatan ikkilamchi bo‘lib qolganini talqin qiladi.
Shu bilan birga, xuddi shu maqolaning mahalliy sezgirlik jadvalida \(\beta_p\) eng katta musbat sezgirlik indeksiga ega. Bu ikki natija bevosita zid bo‘lishi shart emas: biri ma’lum boshlang‘ich nuqtadagi mahalliy hosila sezgirligi, ikkinchisi tanlangan parametrik simulyatsiya oralig‘idir. Ammo natijalarni “atrof-muhit ahamiyatsiz” deb soddalashtirish to‘g‘ri bo‘lmaydi.
Tadqiqot Usuli va Natijalari
DRC epidemiyasi ma’lumotlari bilan parametrlarni baholash
Mualliflar DRC Ebola kuzatuv ma’lumotlaridagi kumulyativ infeksiya va o‘lim qatorlaridan foydalanib modelning ko‘plab parametrlarini baholaydi.
Yangi holat va o‘lim o‘sishlari:
\[ \Delta I_t=I_t-I_{t-1}, \qquad \Delta D_t=D_t-D_{t-1} \]
deb ta’riflangan.
Manba xabar qilingan o‘sishlarda:
- yangi infeksiyalarning o‘rtachasini 4,33, dispersiyasini 40,73, variatsiya koeffitsiyentini 1,47 va qiyshiqligini 2,30;
- yangi o‘limlarning o‘rtachasini 1,81, dispersiyasini 4,46, variatsiya koeffitsiyentini 1,17 va qiyshiqligini 2,18;
- nol hisoblar ulushini infeksiyalarda %33,3, o‘limlarda %28,6
deb xabar qiladi.
Besh hisobot intervalidan iborat siljiydigan oynada dispersiya o‘rtacha bilan oshgani kuzatilgani uchun mualliflar hisob ma’lumotlari geteroskedastik ekanini va log-masshtabdagi moslashtirish katta kuzatuvlarning ta’sirini kamaytirishini ta’kidlaydi.
Log-o‘zgartirilgan og‘irlikli eng kichik kvadratlar
Parametr vektori \(\theta\) uchun ishlatilgan maqsad funksiyasi:
\[ \min_{\theta,I(0),D(0)} \sum_{t=0}^{T} \left[ w_I \left( \log(\hat I_t(\theta)+1)-\log(I_t+1) \right)^2 + \left( \log(\hat D_t(\theta)+1)-\log(D_t+1) \right)^2 \right] \tag{31} \]
ko‘rinishida.
\(\hat I_t(\theta)\) va \(\hat D_t(\theta)\) model bashorat qilgan kumulyativ infeksiya va o‘lim sonlari; \(I_t,D_t\) kuzatuvlardir.
\(+1\) dan foydalanish matematik jihatdan nol hisoblarda \(\log(0)\) aniqlanmasligini oldini oladi.
\(w_I>1\), infeksiya qatorining umumiy optimallashtirishdagi og‘irligini oshirish uchun ishlatilgan. Ammo manbaning asosiy matni \(w_I\) ning aniq sonli qiymatini xabar qilmaydi.
Optimallashtirish bounded nonlinear least squares bilan amalga oshirilgan, ODE tizimi sonli barqarorlik uchun BDF (Backward Differentiation Formula) turidagi stiff solver bilan yechilgan. Manba qaysi dasturiy/paket implementatsiyasi ishlatilganini bildirmaydi.
Parametrlarni baholash va ma’lumotlarni qayta ishlashda vaqt birligi sifatida hafta ishlatilgan, keyingi sonli simulyatsiyalarda esa vaqt birligi kun bo‘lgan. Manba parametrlarni haftadan kunga o‘tkazish tafsilotini alohida o‘zgartirish jadvali bilan bermaydi.
Manba 2-jadvalidagi parametr qiymatlari
| Parametr | Qiymat | Manba holati | Eq. (1)–(2) dagi ma’nosi |
|---|---|---|---|
| \(\beta\) | \(2.508\times10^{-1}\) | Estimated | Yuqumli inson kontakti koeffitsiyenti |
| \(\delta\) | \(2.243\times10^{-1}\) | Estimated | Ta’sirlangandan yuqumli holatga o‘tish |
| \(\gamma\) | \(7.242\times10^{-2}\) | Estimated | Yuqumli sinfdan EVD o‘limi |
| \(\gamma_1\) | \(4.685\times10^{-2}\) | Estimated | Kasalxonaga yotqizish |
| \(\eta\) | \(2.018\times10^{-1}\) | Estimated | Kasalxona sinfidan EVD o‘limi |
| \(\sigma\) | \(1.666\times10^{-5}\) | Estimated | Yuqumli shaxslardan atrof-muhitga sochilish |
| \(\beta_1\) | \(2.937\times10^{-2}\) | Estimated | Kasalxona manbali yuqish |
| \(\beta_2\) | \(2.292\times10^{-3}\) | Estimated | Vafot etgan shaxslardan yuqish |
| \(\beta_3\) | \(1.428\times10^{-3}\) | Estimated | Dafn qilingan shaxslardan yuqish |
| \(\beta_p\) | \(2.317\times10^{-7}\) | Assumed | Atrof-muhit patogen yuqish koeffitsiyenti |
| \(\eta_1\) | \(1.820\times10^{-3}\) | Estimated | Kasalxonadan atrof-muhitga sochilish |
| \(\xi\) | \(1.000\times10^{-10}\) | Estimated | Kasalxonada sog‘ayish |
| \(d\) | \(1.000\) | Assumed | Vafot etgandan dafn qilingan sinfga o‘tish |
| \(b\) | \(2.431\times10^{-1}\) | Assumed | Dafn qilingan sinfdan chiqish |
| \(b_1\) | \(9.785\times10^{-1}\) | Assumed | Dafn qilingan shaxslardan atrof-muhitga sochilish |
| \(d_1\) | \(2.645\times10^{-1}\) | Assumed | Vafot etgan shaxslardan atrof-muhitga sochilish |
| \(\alpha\) | \(7.594\times10^{-1}\) | Estimated | Atrof-muhit patogenining parchalanishi |
Bu jadval muhim talqin chegarasini yaratadi: modelning madaniy/atrof-muhit orqali uzatilishiga kuchli ta’sir qiladigan \(\beta_p\), \(b_1\), \(d_1\), \(d\) va \(b\) parametrlarining bir qismi ma’lumotdan baholanmagan, balki faraz qilingan. Ayniqsa eng yuqori sezgirlik indeksiga ega \(\beta_p\) ning 2-jadvalda Assumed bo‘lishi bu koeffitsiyent real dala ta’sirining bevosita o‘lchangan kattaligi emasligini anglatadi.
Model moslashuv grafigining talqini
Manba 6-shakl modelning kumulyativ infeksiya va o‘lim qatorlarini kuzatuvlar bilan taqqoslaydi. Vizual ravishda model ikkala qatorning umumiy ko‘tarilish tendensiyasini kuzatadi; ammo ayniqsa infeksiya qatorida ayrim hisobot intervallarida kuzatuv nuqtalari bilan model egri chizig‘i orasida ko‘rinadigan og‘ishlar mavjud.
Manbaning asosiy matnida RMSE, MAE, \(R^2\), parametr ishonch oralig‘i, bootstrap noaniqligi yoki mustaqil validatsiya to‘plami xabar qilinmagan. Shuning uchun “model ma’lumotlarga mukammal mos tushdi” degan xulosa chiqarib bo‘lmaydi.
Kalibratsiya ma’lumot manbasidagi bibliografik nomuvofiqlik
Matn ishlatilgan ma’lumotni Kongo Demokratik Respublikasidagi Ebola epidemiyasi deb ta’riflaydi va 6-shakl izohida uni 2022 DRC Ebola epidemiyasi deb ataydi. Biroq adabiyotlar ro‘yxatidagi [57] yozuvining sarlavhasi “Surveillance data for the Ebola outbreak by subcounty, 2020” ko‘rinishida.
Bundan tashqari, matn DRC uchun \(N=47.312.719\) aholi qiymatini [58] bilan bog‘laydi; ammo adabiyotlar ro‘yxatidagi [58] sarlavhasi “Uganda – World Health Organization country profile”.
Bu nomuvofiqliklar kalibratsiyada qaysi ma’lumot yili va qaysi aholi manbasi ishlatilganini mustaqil ravishda qayta ishlab chiqishni qiyinlashtiradi. Verianla bu bibliografik moslikni manbadan tashqari ma’lumot bilan jimlikda tuzatmaydi.
An’anaviy dafn va parvarish yo‘llarining modeldagi ta’siri
Manba \(b_1\), \(d_1\) va \(\sigma\) oshishi atrof-muhit patogen yuklamasini va modelning yuqish salohiyatini oshirishini xabar qiladi.
\(b_1\), dafn qilingan shaxslar bilan bog‘liq atrof-muhitga sochilishni; \(d_1\), hali dafn qilinmagan vafot etganlardan atrof-muhitga sochilishni; \(\sigma\) esa yuqumli shaxslardan atrof-muhitga sochilishni ifodalaydi.
Mualliflar bu yo‘llarni janozada o‘lik jasadni yuvish/teginish va uyda infeksiyalangan shaxsga tana suyuqliklari bilan yaqin kontaktni o‘z ichiga olishi mumkin bo‘lgan parvarish amaliyotlari bilan bog‘laydi. Bu natija model mexanizmiga tegishli; tadqiqot real jamoalarda bu xulq-atvorlarning tez-tez uchrashini alohida dala so‘rovi bilan o‘lchamagan.
Optimal Boshqaruv Muammosida \(u_1\), \(u_2\) va \(u_3\) Nimani Anglatadi?
Mualliflar \(u_1\) ni izolyatsiya va xavfsiz dafn bilan, \(u_2\) ni klinik davolash bilan, \(u_3\) ni atrof-muhit dekontaminatsiyasi bilan bog‘laydi; ammo boshqariladigan ODE tizimidagi matematik joylashuvi bu og‘zaki yorliqlar bilan to‘liq mos kelmagani sababli, ayniqsa \(u_1\) va \(u_2\) ni qayta qo‘llashda manba tenglamalari bilan muallif izohini birgalikda tekshirish kerak.
Optimal boshqaruvning maqsad funksiyasi
Minimallashtirilishi ko‘zlangan umumiy xarajat:
\[ \mathcal J(u_1,u_2,u_3) = \int_0^T \left[ A_1I+A_2H+A_3P+ \frac{c_1}{2}u_1^2+ \frac{c_2}{2}u_2^2+ \frac{c_3}{2}u_3^2 \right]dt \tag{32} \]
ko‘rinishida.
\(A_1I\), yuqumli shaxslar yukini; \(A_2H\), kasalxona yukini; \(A_3P\), atrof-muhit patogen yuklamasini jazolaydi.
\(\frac{c_i}{2}u_i^2\) hadlari esa boshqaruv intensivligini xarajat/sa’y-harakat sifatida jazolaydi. Kvadratik jazo matematik jihatdan boshqaruvni “xarajatsiz har doim 1” qilib tanlashga to‘sqinlik qiladigan optimallashtirish tuzilmasini beradi.
Manba boshqaruv simulyatsiyalarida:
\[ A_1=A_2=A_3=1,\qquad c_1=c_2=c_3=1 \]
dan foydalangan. Shu sababli tadqiqot real sog‘liqni saqlash iqtisodiyoti xarajat ma’lumotlari bilan bajarilgan cost-effectiveness tahlili emas.
Boshqariladigan model tizimi
Manbaning Eq. (33) tenglamasi:
\[ \begin{aligned} \dot S &= \Delta-\Psi S-\mu S,\\ \dot E &= \Psi S-\mu E-\delta(1-u_1)E,\\ \dot I &= \delta(1-u_1)E-(\gamma+\mu)I-\gamma_1(1-u_2)I,\\ \dot H &= \gamma_1(1-u_2)I-(\eta+\xi+\mu)H,\\ \dot R &= \xi H-\mu R,\\ \dot D &= \gamma I+\eta H-dD,\\ \dot B &= dD-bB,\\ \dot P &= \sigma I+\eta_1(1-u_3)H+d_1(1-u_3)D+b_1(1-u_3)B-\alpha P. \end{aligned} \tag{33} \]
Bu yerda \(0\leq u_i(t)\leq1\) boshqaruv intensivligi oralig‘idir.
Boshqaruv yorliqlari bilan tenglamalar orasidagi muhim farq
Manba bayonida \(u_1\) “isolation” sifatida baholanadi va keyin yuqumli shaxslarni jamiyatdan ajratish hamda EVD o‘limlarini xavfsiz dafn bilan boshqarish deb talqin qilinadi. Ammo Eq. (33) da \(u_1\) bevosita \(D\) yoki \(B\) tenglamasiga emas, \(\delta(1-u_1)E\) hadiga qo‘llangan; ya’ni matematik jihatdan ta’sirlangan \(E\) sinfdan yuqumli \(I\) sinfga o‘tishni kamaytiradi.
Xuddi shuningdek \(u_2\) matnda “medical treatment” deb ataladi. Shunga qaramay, Eq. (33) da \(u_2\), \(\gamma_1(1-u_2)I\) orqali yuqumli sinfdan kasalxona sinfiga o‘tishni kamaytiradi; sog‘ayish koeffitsiyenti \(\xi\) ni bevosita oshiruvchi had emas.
\(u_3\) esa Eq. (33) da \(H,D,B\) manbalaridan atrof-muhitga sochilish hadlarini \((1-u_3)\) bilan kamaytiradi va og‘zaki “environmental decontamination” ta’rifi bilan ko‘proq bevosita mos keladi. Shu bilan birga asosiy model Eq. (2) da mavjud tashqi atrof-muhit kirishi \(\phi_p\), boshqariladigan Eq. (33) da yo‘q; manba bu model o‘zgarishining sababini alohida tushuntirmaydi.
Pontryagin Maksimum Prinsipi
Optimal boshqaruv muammosi Pontryagin Maksimum Prinsipi yordamida ko‘rib chiqiladi. Holat o‘zgaruvchilariga mos \(\lambda_S,\lambda_E,\lambda_I,\lambda_H,\lambda_R,\lambda_D,\lambda_B,\lambda_P\) costate/adjoint o‘zgaruvchilari ta’riflangan va terminal shartlar:
\[ \lambda_i(T)=0 \]
deb berilgan.
Umumiy mantiq Hamiltonianning boshqaruv o‘zgaruvchilari bo‘yicha hosilalarini nolga tenglashtirib, har bir vaqtda kasallik yuki bilan aralashuv xarajati orasidagi optimal boshqaruv intensivligini topishdir.
Eq. (36) bilan Eq. (38) orasidagi ishora nomuvofiqligi
Manba Eq. (36) da proyeksiya formulasi ichida \(u_1\) uchun:
\[ \frac{\delta E^*(\lambda_E-\lambda_I)}{c_1} \]
va \(u_2\) uchun:
\[ \frac{\gamma_1 I^*(\lambda_I-\lambda_H)}{c_2} \]
ifodalarini ko‘rsatadi.
Ammo darhol keyin optimality condition yechilayotganda Eq. (38) da:
\[ \hat u_1= -\frac{\delta E(\lambda_E-\lambda_I)}{c_1} \]
va:
\[ \hat u_2= -\frac{\gamma_1I(\lambda_I-\lambda_H)}{c_2} \]
deb yozilgan.
Demak, dastlabki ikki boshqaruv uchun Eq. (36) va Eq. (38) orasida ishora farqi bor. Bu kodlash nuqtayi nazaridan muhim farq. Manba kodi yoki muallif tasdig‘isiz Verianla bu ikki shakldan birini “to‘g‘ri” deb tanlamaydi.
Optimal boshqaruv simulyatsiyasining boshlang‘ich shartlari
| Holat | Boshlang‘ich qiymat |
|---|---|
| \(S(0)\) | 47.312.419 |
| \(E(0)\) | 0 |
| \(I(0)\) | 151 |
| \(H(0)\) | 11 |
| \(R(0)\) | 88 |
| \(D(0)\) | 50 |
| \(B(0)\) | 0 |
| \(P(0)\) | 0 |
Inson holatlarining yig‘indisi \(47.312.719\) bo‘lib, manbada ishlatilgan umumiy aholi qiymatiga teng.
Yakka, juft va uchlik boshqaruv natijalari
| Strategiya | Manbada xabar qilingan nisbiy natija | Talqin chegarasi |
|---|---|---|
| \(u_1\) yakka | Mualliflar yakka boshqaruvlar ichida eng kuchli deb talqin qilgan ssenariy | \(u_1\) ning og‘zaki “isolation” yorlig‘i bilan Eq. (33) dagi matematik o‘rni to‘liq mos kelmaydi. |
| \(u_2\) yakka | Infeksiya, kasalxona yuki va patogen yukini kamaytiradi; \(u_1\) dan kuchsizroq | Real klinik davolash samaradorligi emas. |
| \(u_3\) yakka | Atrof-muhit patogen yukini kamaytiradi; umumiy ta’sir \(u_1\) dan kuchsizroq | Atrof-muhitni tozalashning o‘zi yetarli deb hisoblanmaydi. |
| \(u_1+u_2\) | Juft strategiyalar orasida eng kuchli deb talqin qilinadi | Manba grafiklarining model shartlariga xos natija. |
| \(u_1+u_3\) | Sezilarli kamayish; inson va atrof-muhit yo‘lini birga nishonga oladi | Asosiy matnda sonli foiz kamayishi berilmagan. |
| \(u_2+u_3\) | Samarali, ammo \(u_1\) ni o‘z ichiga olgan juft strategiyalardan kuchsizroq | Modelda ochiq izolyatsiya komponenti yo‘qligi bilan bog‘lanadi. |
| \(u_1+u_2+u_3\) | O‘rganilgan ssenariylar orasida infeksiya, kasalxona va atrof-muhit yukida eng kuchli bostirish | Taxminan 20 kunlik model simulyatsiyasi; dala tajribasi emas. |
Uchta boshqaruv birgalikda qo‘llanganda model nimani ko‘rsatadi?
Manba 12-shaklda aralashuvsiz ssenariyda yuqumli shaxslar, kasalxona yuki va atrof-muhit patogen konsentratsiyasi 20 kunlik simulyatsiya davomida tez ko‘tarilgani, uchta boshqaruv birgalikda qo‘llanganda esa barcha egri chiziqlar pastroq darajalarda qolganini ko‘rsatadi.
Mualliflarning mexanistik talqini bir nechta yuqish zanjirini bir vaqtda nishonga olish qayta aloqa ta’sirini yaratishidir: yuqumli shaxslar kamroq bo‘lsa atrof-muhitga kamroq patogen chiqaradi; atrof-muhit patogen yuklamasi pastroq bo‘lsa yangi infeksiya bosimi ham kamayadi.
Bu izoh modelning o‘z tenglamalari ichida izchil; ammo real epidemiyada aynan shunday kattalikdagi ta’sir yuz berishini kafolatlamaydi.
Tadqiqot qo‘llab-quvvatlaydigan natijalar
- Inson, kasalxona, o‘limdan keyingi, dafn qilingan jasad va atrof-muhit yo‘llari bitta Ebola ODE modelida matematik tarzda ifodalanishi mumkin.
- Manba modelida \(\mathcal R_0\) bevosita va bilvosita ko‘plab yuqish yo‘llaridan tashkil topgan.
- \(\phi_p=0\) shartida \(\mathcal R_0\leq1\) kasalliksiz muvozanat barqarorligi bilan bog‘langan.
- Manba sezgirlik jadvalida \(\beta_p,\eta,\eta_1,\sigma\) kuchli musbat mahalliy sezgirlik ko‘rsatadi.
- \(\alpha,\xi,b,d,\gamma_1\) manfiy sezgirlik ko‘rsatadi.
- Tanlangan atrof-muhit parametr ssenariylarida atrof-muhit yo‘li o‘zi yolg‘iz inson kontakti yo‘llaridan kuchsizroq ko‘rinadi.
- Mualliflar \(u_1\) deb belgilagan boshqaruv yakka ssenariylar ichida eng kuchli deb xabar qilingan.
- Uchta boshqaruvning birgalikda qo‘llanishi simulyatsiyada o‘rganilgan ssenariylar ichida eng kuchli umumiy bostirishni hosil qilgan.
Tadqiqot qo‘llab-quvvatlamaydigan natijalar
- Model Afrikadagi har qanday jamiyatning barcha an’anaviy amaliyotlari Ebola yuqishini oshirishini isbotlamaydi.
- Model parametrlari madaniy xulq-atvor prevalensiyasining bevosita dala o‘lchovi emas.
- Simulyatsiya natijalari real epidemiyada aynan shu aralashuv ta’siri albatta yuz berishini isbotlamaydi.
- \(\beta_p\) ning yuqori sezgirligi atrof-muhit yuqishi barcha sharoitlarda real dunyodagi dominant yuqish yo‘li ekanini isbotlamaydi.
- 20 kunlik optimal boshqaruv grafiklari klinik yoki jamoa asosidagi randomizatsiyalangan aralashuv tadqiqoti emas.
- Model makoniy tarqalishni, shaharlararo harakatni yoki chegaralararo mobillikni hisoblamaydi.
- Model an’anaviy xulq-atvorlarning epidemiya davomida o‘zgarishini dinamik modellashtirmaydi.
- Manba barcha parametrlar uchun ishonch oralig‘i yoki identifiability tahlilini taqdim etmaydi.
- Manba kodi berilmagani sababli Eq. (36)/Eq. (38) ishora nomuvofiqligi amalda qanday hal qilingani tasdiqlanmaydi.
- Modelning moslashuv grafigi o‘zi yolg‘iz tashqi namunadagi bashorat muvaffaqiyatini ko‘rsatmaydi.
Tadqiqotning asosiy cheklovlari
Bir jinsli aralashish: Har bir shaxsning kontakt tuzilmasi bitta epidemiologik hovuz ichida bo‘lgandek soddalashtirilgan.
Makoniy geterogenlik yo‘q: Turli qishloq, shahar, sog‘liq hududi yoki mamlakat orasidagi geografik yuqish modellashtirilmagan.
Mobillik yo‘q: Inson harakatchanligi ochiq o‘zgaruvchi emas.
Madaniy xulq-atvorlar doimiy parametrlar bilan ifodalanadi: jamiyat xulq-atvorining vaqt bilan o‘zgarishi, aralashuvni qabul qilish yoki turli jamoalarning geterogenligi alohida dinamikalar emas.
Ba’zi markaziy parametrlar faraz qilingan: ayniqsa \(\beta_p,b,d,b_1,d_1\) manba 2-jadvalda Assumed sifatida keltirilgan.
Kalibratsiya manbalari orasida bibliografik nomuvofiqlik mavjud: ma’lumot yili, DRC/Uganda aholi manbasi va 6-shakl izohi to‘liq mos emas.
Fit noaniqligi yetishmaydi: parametr ishonch oraliqlari, profile likelihood, Bayesian posterior yoki boshqa formal identifiability/uncertainty tahlili xabar qilinmagan.
Boshqaruv xarajatlari real iqtisodiy ma’lumot emas: \(A_i\) va \(c_i\) qiymatlari simulyatsiyada 1 deb olingan.
Boshqaruv tenglamasi bilan og‘zaki mexanizm to‘liq mos kelmaydi: ayniqsa \(u_1,u_2\) ni qayta qo‘llash ehtiyotkorlikni talab qiladi.
Optimal boshqaruv formulalarida ishora nomuvofiqligi bor: Eq. (36) bilan Eq. (38) dagi dastlabki ikki boshqaruv ifodasi bir xil ishoraga ega emas.
Atrof-muhit tashqi kirishi boshqaruv modelida yo‘qoladi: Eq. (2) dagi \(\phi_p\), Eq. (33) da mavjud emas va manba bu o‘zgarish sababini alohida bermaydi.
Manba va Usul Haqida Izoh
Original sarlavha: Modeling the role of the environment and African tradition practices in the transmission dynamics of Ebola virus disease
Mualliflar: Ijuptil Kwajighu Joseph; Georg Bader; Nicholas Kwasi-Do Ohene Opoku.
Mas’ul mualliflar: Ijuptil Kwajighu Joseph va Nicholas Kwasi-Do Ohene Opoku manbada yulduzcha bilan corresponding author sifatida ko‘rsatilgan.
Teng hissa: Ko‘rib chiqilgan manbada teng birinchi mualliflik/teng hissa bayonoti mavjud emas.
Muassasa: African Institute for Mathematical Sciences, Accra, Ghana.
Muharrir: Mohamed Fathy El-Amin Mousa.
Jurnal: Scientific African.
Jild: 32.
Maqola raqami: e03392.
Nashr davri: Rasmiy ScienceDirect yozuvi Volume 32, June 2026 deb ko‘rsatadi.
DOI: 10.1016/j.sciaf.2026.e03392.
Qabul qilingan sana: 12 Dekabr 2025.
Qayta ko‘rib chiqilgan sana: 23 Mart 2026.
Qabul qilingan sana: 6 May 2026.
Onlayn mavjudlik: 11 May 2026.
Nashriyot: Elsevier B.V.
Platforma: ScienceDirect.
Manba turi: Matematik epidemiologiya tadqiqoti; deterministik chiziqli bo‘lmagan ODE modellashtirish, ma’lumot kalibratsiyasi, barqarorlik/bifurkatsiya, sezgirlik va optimal boshqaruv tahlili.
Taqriz holati: Scientific African jurnalida chop etilgan, qabul qilish–qayta ko‘rib chiqish–qabul qilish sanalari mavjud jurnal maqolasi; PDF batafsil taqriz hisobotlarini o‘z ichiga olmaydi.
Litsenziya: Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International — CC BY-NC-ND 4.0.
Mualliflik huquqi/vizual izoh: ND sharti sababli manba shakllarining vizual kompozitsiyasi tarjima qilinib yoki qayta uslublashtirilib moslashtirilmagan. 1-shakl shuningdek uchinchi tomon [7] manbasiga asoslanganligi sababli qayta ishlatilmagan. Maqoladagi SVG faqat tenglamalarda ochiq belgilangan model bog‘lanishlaridan butunlay mustaqil kompozitsiyada yaratilgan.
Model o‘zgaruvchilari: \(S,E,I,H,R,D,B,P\).
Asosiy usullar: musbatlik va boundedness tahlili; next-generation matrix; \(\mathcal R_0\); Jacobian mahalliy barqarorlik; Lyapunov–LaSalle global barqarorlik; endemik muvozanat; center-manifold bifurkatsiya tahlili; normalize forward sensitivity; log-o‘zgartirilgan og‘irlikli bounded nonlinear least squares; BDF stiff ODE solver; Pontryagin Maximum Principle.
Ma’lumot konteksti: Maqola DRC Ebola kuzatuv ma’lumotlaridan foydalanilganini bildiradi. 6-shakl izohi moslashuvni “2022 Ebola outbreak in the Democratic Republic of the Congo” deb ta’riflaydi. Biroq [57] havolasining sarlavhasida 2020 ifodasi mavjud. Bu nomuvofiqlik manbada hal qilinmagan.
Aholi manbasi nomuvofiqligi: Matn \(N=47.312.719\) qiymatini DRC aholisi sifatida [58] bilan bog‘laydi; adabiyotlar ro‘yxatidagi [58] yozuvining sarlavhasi “Uganda – World Health Organization country profile”.
Vaqt birligi: Parametrlarni baholash va ma’lumotlarni qayta ishlash haftalarda, sonli simulyatsiyalar kunlarda olib borilgan. O‘tkazish tafsiloti alohida jadval yoki algoritm sifatida xabar qilinmagan.
Fitting og‘irligi: \(w_I>1\) ekani bildirilgan; aniq \(w_I\) qiymati asosiy matnda berilmagan.
Model kodi: Ko‘rib chiqilgan PDFda manba kod ombori yoki bajariladigan model implementatsiyasiga havola xabar qilinmagan.
Parametr noaniqligi: Asosiy matnda baholangan parametrlar uchun ishonch oralig‘i, posterior taqsimot yoki formal identifiability tahlili xabar qilinmagan.
Manba ichidagi parametr nomi nomuvofiqligi: Eq. (1) da \(\beta,\beta_1,\beta_2,\beta_3,\beta_p\) mos ravishda yuqumli, kasalxona, vafot etgan, dafn qilingan va atrof-muhit manbalariga berilgan bo‘lsa-da, natijalar bayonida indekslar boshqacha belgilangan.
Manba ichidagi mexanizm talqini: \(\sigma\) matematik modelda yuqumli shaxslardan atrof-muhitga sochilish koeffitsiyenti; sezgirlik muhokamasida esa bevosita simptomatik kontakt dalili sifatida talqin qilingan. Bu ikki mexanizm teng emas.
Manba ichidagi boshqaruv ta’rifi nomuvofiqligi: \(u_1\) og‘zaki ravishda isolation/safe burial bilan, \(u_2\) medical treatment bilan izohlangan bo‘lsa-da, boshqariladigan Eq. (33) da bu boshqaruvlar mos ravishda \(E\rightarrow I\) va \(I\rightarrow H\) o‘tish koeffitsiyentlarini kamaytiradi.
Manba ichidagi optimal boshqaruv ishora nomuvofiqligi: Eq. (36) bilan Eq. (38) da \(u_1\) va \(u_2\) uchun berilgan ifodalarning ishoralari bir-biriga mos emas.
Manba ichidagi atrof-muhit modeli farqi: Asosiy Eq. (2) da \(\phi_p\) tashqi atrof-muhit patogen kirishi mavjud bo‘lsa, optimal boshqaruv Eq. (33) da bu had olib tashlangan; manba sababini alohida tushuntirmaydi.
Bifurkatsiya matni: Manba \(\phi_p=0\) uchun oldinga bifurkatsiya ko‘rsatganini aytgan holda xuddi shu bo‘limda “\(\phi_p\geq0\) uchun bifurkatsiya yo‘q” ifodasini ham ishlatadi. Bu ifoda avvalgi bayon bilan to‘liq mos emas va Verianla tomonidan jimlik bilan tuzatilmagan.
Sezgirlik jadvali: Modelda mavjud \(d_1\) parametri uchun 1-jadvalda normallashtirilgan sezgirlik indeksi ko‘rinmaydi.
CRediT — Ijuptil Kwajighu Joseph: Collected the data; Contributed data or analysis tools; Performed the analysis; Wrote the paper; Other contribution.
CRediT — Georg Bader: Conceived and designed the analysis; Other contribution.
CRediT — Nicholas Kwasi-Do Ohene Opoku: Conceived and designed the analysis; Other contribution.
Moliyalashtirish: Mualliflar tadqiqot uchun moliyalashtirish olmaganini bildirgan.
Manfaatlar to‘qnashuvi: Mualliflar tadqiqotga ta’sir qilishi mumkin bo‘lgan ma’lum raqobat qiluvchi moliyaviy manfaat yoki shaxsiy munosabat mavjud emasligini bildirgan.
Minndorchilik: Birinchi muallif ayniqsa African Institute for Mathematical Sciences ga tadqiqotni boshlash uchun berilgan asos va imkoniyat sababli minnatdorlik bildiradi.
Ma’lumotlar mavjudligi: Manba natijalarni qo‘llab-quvvatlovchi ma’lumotlar tadqiqot ichida tegishli tarzda havola qilinganini bildiradi.
Sun’iy intellekt bayonoti: Ko‘rib chiqilgan asosiy PDFda generativ sun’iy intellektdan foydalanish bo‘yicha alohida bayonot mavjud emas.
Asosiy talqin chegarasi: Ushbu tadqiqot Ebolaning real dunyodagi barcha yuqish mexanizmlarini bevosita o‘lchamaydi. Matematik model ma’lum holat bo‘limlari, kalibrlangan/faraz qilingan parametrlar va bir jinsli aralashish farazi ostida inson, o‘limdan keyingi va atrof-muhit yuqish yo‘llarining xatti-harakatini o‘rganadi. Optimal boshqaruv natijalari simulyatsiya natijasidir; real epidemiyada o‘lchangan aralashuv samaradorligi emas.
Izoh qoldiring
E-pochta manzilingiz chop etilmaydi. Majburiy maydonlar * bilan belgilangan