Generalized Immunity changes

HeLab · Apr 12, 2022 · 8c4bff2 · 8c4bff2
1 parent 679f3e2
commit 8c4bff2
Showing 1 changed file with 10 additions and 9 deletions.
diff --git a/batchRuns/ODESGeneralizedImmunity.R b/batchRuns/ODESGeneralizedImmunity.R
@@ -771,15 +771,16 @@ allPara<-read.csv(paramFile,sep=",",header=T)
 
 currentPara<-allPara[allPara$No==No,]
 
-t2 = seq(from=0,to=30*365,by=1) 
-N0 = c(PW_NU=100,PW_ND=0,PW_SU=0,PW_SD=0,PW_AU=0,PW_AD=0,
-       PR_NU=0,PR_ND=0,PR_SU=0,PR_SD=0,PR_AU=0,PR_AD=0,
-       NU=10000,ND=0,SU=0,SD=0,AU=0,AD=0,
-       b=allPara[1,]$bavgnull*(1-allPara[1,]$seasonality))
-for (i in 1:42){
-  out9 = ode(y=N0,times=t2,func=SEAI_p,parms=allPara[i,])
-  out9_freqRes[i] = (out9[,8]+out9[,9]+out9[,10]+out9[,11]+out9[,12]+out9[,13])/(out9[,2]+out9[,3]+out9[,4]+out9[,5]+out9[,6]+out9[,7]+out9[,8]+out9[,9]+out9[,10]+out9[,11]+out9[,12]+out9[,13]) 
-}
+#t2 = seq(from=0,to=30*365,by=1) 
+#N0 = c(PW_NU=100,PW_ND=0,PW_SU=0,PW_SD=0,PW_AU=0,PW_AD=0,
+#       PR_NU=0,PR_ND=0,PR_SU=0,PR_SD=0,PR_AU=0,PR_AD=0,
+#       NU=10000,ND=0,SU=0,SD=0,AU=0,AD=0,
+#       b=allPara[1,]$bavgnull*(1-allPara[1,]$seasonality))
+#for (i in 1:42){
+#  out9 = ode(y=N0,times=t2,func=SEAI_p,parms=allPara[i,])
+#  out9_freqRes[i] = (out9[,8]+out9[,9]+out9[,10]+out9[,11]+out9[,12]+out9[,13])/(out9[,2]+out9[,3]+out9[,4]+out9[,5]+out9[,6]+out9[,7]+out9[,8]+out9[,9]+out9[,10]+out9[,11]+out9[,12]+out9[,13]) 
+#}
+