load("bengal.dat") # data available on Moodle page bengal # look at event times u <- (bengal-1877)/101 # rescale them to (0,1) and plot their empirical CDF plot(u,c(1:141)/141,type="s",panel.first=abline(0,1,col="grey")) ks.test(u,y="dunif") # Kolmogorov-Smirnov test of uniformity