From 485ebd6d5b34ff2d8f6761f2dda44ea3e2b32aa6 Mon Sep 17 00:00:00 2001 From: Ben Green Date: Sun, 30 Aug 2015 20:59:39 -0400 Subject: wrote script to find background rate using nls --- R Scripts/sim-analysis.R | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) (limited to 'R Scripts/sim-analysis.R') 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)) -- cgit v1.2.3-70-g09d2