Fig3-Histogram

#################################################################################################################################################
##### R script for producing Figure 3 of Garcia-Diaz et al. 'A guide to developing and using quantitative models in conservation management' ####
#################################################################################################################################################

### Define the number of values to simulate (1000 in this case)
sample.size<-1000

### Simulate 1000 values from a Poisson distibution with mean = 5
sim1<-rpois(sample.size, lambda=5)

###  Produce 1000 values, with values spanning the same range as 'sim1'. 

sim2<-c(rep(min(sim1), 150), rep(min(sim1)+2, 150), rep(5, 200), rep(5+2, 100), rep(10, 100), rep(max(sim1)-2, 100), rep(max(sim1), 200))
mean(sim2)

length(sim2)

###  Plotting the histrograms in a multipanel figure

### Define the number and position of the two panels ('mfrow'), and the margins of the figures

par(mfrow=c(1, 2), mar=c(5, 7, 5, 6))

#### Histogram of the Poisson distributed values with custom figure specifications (refer to the R documentatoin and help for details)
hist(sim1, axes=FALSE, xlab="", ylab="", main="")

axis(1, at =  round(seq(from=min(sim1), to=max(sim1), length.out=5), digits=0), cex.axis= 1.5, font=1, family='serif', las=1)
axis(2, at = c (0, 50, 100, 150, 200), cex.axis = 1.5, las=2, font=1, family='serif')
mtext("Frequency", side=2, line=4, cex=2, font=1, family='serif')
mtext("Poisson distribution", side=3,line=1, cex=3, font=1, family='serif')

box()

#### Histogram,with custom figure specifications, of the other 1000 values not following a distribution  (refer to the R documentation and help for details)
hist(sim2, axes=FALSE, xlab="", ylab="", main="")

axis(1, at =  round(seq(from=min(sim2), to=max(sim2), length.out=5), digits=0), cex.axis= 1.5, font=1, family='serif', las=1)
axis(2, at = c (0, 50, 100, 150, 200), cex.axis = 1.5, las=2, font=1, family='serif')
mtext("Frequency", side=2, line=4, cex=2, font=1, family='serif')
mtext("No probability distribution", side=3,line=1, cex=3, font=1, family='serif')

box()