1 van 1

Model infectieziekten

Geplaatst: za 21 mar 2009, 13:14
door PhilipVoets
Geachte heer/mevrouw,

Voor mijn Geneeskundescriptie doe ik onderzoek naar infectieziekten. Ik probeer voor de aardigheid een model op te stellen dat het aantal geïnfecteerde personen
\(I\)
als functie van tijd
\(t\)
beschrijft. Hiervoor zijn onder meer het basis-reproductiegetal
\(R\)
en de ziektevatbaarheid
\(S\)
van belang. De betekenissen hiervan - voor zover relevant - zijn respectievelijk het gemiddeld aantal secundaire besmettingen door een geïnfecteerd individu en de fractie personen van een populatie die bevattelijk is voor de ziekte. Om
\(R\)
secundaire besmettingen te veroorzaken is gemiddeld gezien de tijd
\(T\)
nodig. In het irreële geval van een oneindig grote populatie leek mij dat het model er simpelweg als volgt uit moest ziet:
\( I(t) = I(0)(RS)^{t/T} = I(0)e^{ln(RS)t/T} \)


Nu is er echter in werkelijkheid een maximale populatieomvang
\(N\)
en vanzelfsprekend geldt
\(I(t) < N\)
. Ik was niet zeker hoe ik bovenstaande formule zodanig uit moest breiden dat de term
\(N\)
hierin meegenomen werd en er dus een realistischer beeld van de situatie ontstaat. Mijn poging was de volgende.

Omdat bovenstaande vergelijking in zijn laatste vorm gezien kan worden als de oplossing van de differentiaalvergelijking:
\( dI/dt = ln(RS)I/T \)
Leek mij dat ik door deze diffentiaalvergelijking uit te breiden naar een vorm voor begrensde groei, oftewel:
\( dI/dt = ln(RS)I(N-I)/T \)
En hier de oplossing voor te zoeken (Verhulst), het probleem zou hebben opgelost. Mijn vraag is echter of ik dit uitbreiden van de differentiaalvergelijking ongestraft mag doen. Ik kom dan uit op de volgende oplossing (enigszins vereenvoudigd):
\( I(t) = I(0)N/(I(0) + (N-I(0))(RS)^{-Nt/T}) \)
Op zich voldoet deze formule aan een aantal criteria dat er aan gesteld wordt (startwaarde
\(I(0)\)
, maximale waarde
\(N\)
, groei afhankelijk van het product
\(RS\)
e.d.), maar toch heb ik mijn twijfels.

Suggesties worden zeer op prijs gesteld.

Vriendelijke groet en bij voorbaat dank,

Philip

Re: Model infectieziekten

Geplaatst: zo 22 mar 2009, 09:33
door PeterPan
Je definitie voor
\(R\)
R = het gemiddeld aantal secundaire besmettingen door een geïnfecteerd individu,

kan niet juist zijn.

Het moet zonder twijfel zijn:
\(R = \)
het gemiddeld aantal mogelijke secundaire besmettingen door een geïnfecteerd individu.

(Als
\(S\approx 0\)
, dan zou in jouw definitie
\(R\approx 0\)
gelden. R en S behoren onafhankelijk te zijn.)
\(I(t) < N\)
is te grof. Er geldt
\(I(t) < SN\)
.
\( dI/dt = ln(RS)I(SN-I)/T \)
is een eenvoudige uitbreiding van de formule die een functie
\(I\)
oplevert die de grens
\(I=SN\)
niet overschrijdt, maar hem wel willekeurig dicht benadert.

Dat is een heel bijzondere eigenschap. Je mag proberen andere aanpassingen aan de beginformule te maken, maar dat zal tot geen resultaat lijden, omdat dit asymptotische gedrag heel uitzonderlijk is.

De eenvoud is het kenmerk van het ware.

Re: Model infectieziekten

Geplaatst: zo 22 mar 2009, 14:39
door PhilipVoets
Ja, je hebt gelijk. Ik ben iets te overhaast te werk gegaan. Mijn herziene suggestie luidt:
\( dI/dt = ln(RS)I(1-I/SN)/T \)
Hetgeen leidt tot de oplossing:
\( I(t) = I(0)SN/(I(0)+(SN-I(0))(RS)^{-t/T}) \)
Kun jij je hier ook in vinden?

Dit model houdt er echter geen rekening mee dat de ziektevatbaarheid ook in het verloop van de tijd verandert...

Re: Model infectieziekten

Geplaatst: zo 22 mar 2009, 15:23
door PeterPan
Als ik het zo snel zie is
\( I(t) = I(0)SN/(I(0)+(SN-I(0))(RS)^{-t/T}) \)
oplossing van de differentiaalvergelijking
\( dI/dt = \ln(RS)I(SN-I)/T \)


en
\( dI/dt = \ln(RS)I(1-I/SN)/T \)
levert de oplossing
\( I(t) = I(0)SN/(I(0)+(SN-I(0))(RS)^{-SNt/T}) \)
en dat lijkt me niet correct.