diff --git a/batchRuns/ODESGeneralizedImmunity.R b/batchRuns/ODESGeneralizedImmunity.R index e51265c..73edf1a 100644 --- a/batchRuns/ODESGeneralizedImmunity.R +++ b/batchRuns/ODESGeneralizedImmunity.R @@ -108,12 +108,12 @@ SEAI_p<-function(t,y,p){ NU2ND = (Bu*(1-rNU/capacity)*infN*d1) #naive U class getting enough # of infections to move to SU - NU2SU = (Bu*(1-rNU/capacity)*rho1*1/((1/(infN*muWU))+(1/(1-infN)))) + NU2SU = (Bu*(1-rNU/capacity)*rho1*1/((infN/(muWU))+((1-infN)))) #naive U class die from infections NUInfDie = (Bu*(1-rNU/capacity)*infN*muWU*muNU) #naive D class getting enough # of infections to move to SD - ND2SD = (Bd*(1-rND/capacity)*rho1*1/((1/(infN*muRU))+(1/(1-infN)))) + ND2SD = (Bd*(1-rND/capacity)*rho1*1/((infN/(muRU))+((1-infN)))) #naive D class die from infections NDInfDie = (Bd*(1-rND/capacity)*infN*muRD*muND) @@ -122,10 +122,10 @@ SEAI_p<-function(t,y,p){ SU2SD = (Bu*(1-rSU/capacity)*infS*d2) #SU class getting enough # of infections to move to AU - SU2AU = (Bu*(1-rSU/capacity)*rho2*1/((1/(infS*muWU))+(1/(1-infS)))) + SU2AU = (Bu*(1-rSU/capacity)*rho2*1/((infS/(muWU))+((1-infS)))) #SD class getting enough # of infections to move to AD - SD2AD = (Bd*(1-rSD/capacity)*rho2*1/((1/(infS*muRD))+(1/(1-infS)))) + SD2AD = (Bd*(1-rSD/capacity)*rho2*1/((infS/(muRD))+((1-infS)))) #AU class getting infected, show symptoms, and treated AU2AD = (Bu*(1-rAU/capacity)*infA*omega*d3)