diff options
Diffstat (limited to 'R Scripts/sim-analysis.R')
| -rwxr-xr-x | R Scripts/sim-analysis.R | 8 |
1 files changed, 4 insertions, 4 deletions
diff --git a/R Scripts/sim-analysis.R b/R Scripts/sim-analysis.R index e368419..af280cc 100755 --- a/R Scripts/sim-analysis.R +++ b/R Scripts/sim-analysis.R @@ -31,8 +31,8 @@ for (i in 1:n.infections){ n = 10000 mean.time = rep(NA,n) med.time = rep(NA,n) +mean.30 = rep(NA,n) mean.50 = rep(NA,n) -mean.100 = rep(NA,n) n.vicpairs = rep(NA,n) hyp_lcc = upgrade_graph(hyp_lcc) @@ -56,16 +56,16 @@ for(q in 1:n){ n.vicpairs[q] = length(vic.time) mean.time[q] = mean(vic.time) med.time[q] = median(vic.time) + mean.30[q] = mean(vic.time<=30) mean.50[q] = mean(vic.time<=50) - mean.100[q] = mean(vic.time<=100) # decay.rate[sim,q] = getDecayRate(vic.time) #get lm of histogram. or maybe just cor? } print(proc.time()-ptm) -# save(mean.time,med.time,mean.50,mean.100,file='Results/sims-Aug28.RData') +# save(mean.time,med.time,mean.30,mean.50,file='Results/sims-Aug28.RData') ####### plot simulation data ####### -data = 0.1179039 # 778.2056, 638, 0.07248908, 0.1179039 +data = 0.1179039 # 778.2056, 638, 0.05105811, 0.07248908 simdata = mean.100 xl = c(0.95*min(min(simdata),data),1.05*max(max(simdata),data)) yl = c(0,max(hist(simdata,25)$density)) |