Fixed Gen Immunity Script

HeLab · Apr 12, 2022 · 98ba5cc · 98ba5cc
1 parent abd4b7a
commit 98ba5cc
Showing 1 changed file with 4 additions and 4 deletions.
diff --git 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)