Advanced Search

Question: How to compute sensitivity analysis on basic reproduction number?

Question:How to compute sensitivity analysis on basic reproduction number?

Yee Voon 20 Maple
sensitivity
March 04 2018
0

Basic_Reproduction_Number.mw
Hello, anybody knows how to compute sensitivity analysis for basic reproduction number using Maple coding? By taking the partial derivatives of all the parameters involved in the basic reproduction number using normalised forward sensitivity index formula. Thanks.

> 

restart
> 

unprotect(Pi)
> 

R[0] := k*tau*[(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q))]*(Pi*(-mu*p+mu+phi)/(mu*(vartheta+mu+phi))+epsilon*Pi*(mu*p+vartheta)/(mu*(vartheta+mu+phi)))

k*tau*[(rho*(Upsilon*(mu+alpha+eta)+chi)+(1-rho)*(Upsilon*(1-q)*eta+mu+beta+chi))/((mu+beta+chi)*(mu+alpha+eta)-chi*eta*(1-q))]*(Pi*(-mu*p+mu+phi)/(mu*(vartheta+mu+phi))+epsilon*Pi*(mu*p+vartheta)/(mu*(vartheta+mu+phi)))

 (1)
> 

``


 

Download Basic_Reproduction_Number.mw

 
Please Wait...