Estimating arrival times of people in a shop using R. Part 3. Full article at: http://www.firsttimeprogrammer.blogspot.com/2015/07/estimating-arrival-times-of-people-in.html

arrivalTimes3.R

#-------------------------------------------------------------------------------
#                           Compare simulated to obseved data
#-------------------------------------------------------------------------------

# Let's compare observed data with simulated data

# Comparative boxplot
boxplot(interarrivals,simulated.min,xlab='Minutes',col=c('cyan','chartreuse'),
        border=c('blue','seagreen'),names=c('Observed','Simulated'),
        main='Interarrival times',notch=TRUE,horizontal=TRUE)

# Comparative histograms
hist(simulated.min,col='chartreuse',breaks=breaks.points1,xlab=x.lab,ylab=y.lab,
     main=main.lab,border='seagreen',density=20,angle=30)
hist(interarrivals,col='cyan',breaks=breaks.points1,border='blue',
     density=20,angle=160,add=TRUE)

# Add a legend
legend("topright",c("Simulated","Observed"), col=c("seagreen", "blue"), lwd=10)

# Comparative histograms of number of people per arrival
hist(people.per.arr,xlab='No. of people',breaks=break.points2,
     main='Simulated number of people per arrival',col='chartreuse',
     border='seagreen',density=30,angle=30)
hist(data$pers,breaks=break.points2,col='red',border='red4',density=30,
     angle=160,add=TRUE)

# Add a legend
legend("topright",c("Simulated","Observed"), col=c("seagreen", "red"), lwd=10)