1 files changed, 7 insertions, 5 deletions
diff --git a/R Scripts/sim-analysis.R b/R Scripts/sim-analysis.R
index e56f575..51659f5 100755
--- a/R Scripts/sim-analysis.R
+++ b/R Scripts/sim-analysis.R
@@ -68,18 +68,20 @@ colnames(simdata)=c('mean.time','med.time','mean.30','mean.100')
 ####### plot simulation data #######
 # mean=662.6, med=488, m30=0.076, m60=0.119, m100=0.167
 d = c(662.6,488,0.076,0.167)
-xlabs = c('Mean time between infections',
-          'Median time between infections',
+xlabs = c('Mean days between infections',
+          'Median days between infections',
           'Infections within 30 days',
           'Infections within 100 days')
+xlims = matrix(c(660,765,475,625,0.03,0.08,0.1,0.18),ncol=2,byrow=T)
 par(mfrow=c(2,2))
 for(i in 1:4){
   data = d[i] 
   sdata = simdata[,i]
-  xl = c(0.99*min(min(sdata),data),1.01*max(max(sdata),data))
+  xl = xlims[i,]
   h = hist(sdata,50,xlim=xl,col='#1f78b4',freq=T,border=NA,axes=F,
-       xlab=xlabs[i],ylab='Relative Frequency',main=NULL)
-  axis(1,at=pretty(xl,5))
+       xlab=xlabs[i],ylab='Relative frequency',main=NULL)
+  if(i %in% c(1,2)) axis(1,at=pretty(xl,5))
+  if(i %in% c(3,4)) axis(1,at=pretty(xl,5),lab=paste(pretty(xl,5)*100,'%',sep=''))
   axis(2,at=pretty(h$counts,3),lab=pretty(h$counts/n,3))
   abline(v=data,lwd=4,col='#e41a1c')
   box(lwd=1.1)