diff options
| author | Ben Green <bgreen@g.harvard.edu> | 2015-08-30 20:59:39 -0400 |
|---|---|---|
| committer | Ben Green <bgreen@g.harvard.edu> | 2015-08-30 20:59:39 -0400 |
| commit | 485ebd6d5b34ff2d8f6761f2dda44ea3e2b32aa6 (patch) | |
| tree | 96e610b53c18f73a155974964e17f98bd1ade871 /R Scripts/sim-analysis.R | |
| parent | d1e9e8cfabfee86fd7ac9e108fab0d463da24ce8 (diff) | |
| download | criminal_cascades-485ebd6d5b34ff2d8f6761f2dda44ea3e2b32aa6.tar.gz | |
wrote script to find background rate using nls
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)) |