# NOT RUN {
data("SeymourIslandBins")
#Compare bins from Macellari and Witts before the K/Pg extinction
plot(SeymourIslandBins$Macellari.1984.Section.D.,
rep(0,length(SeymourIslandBins$Macellari.1984.Section.D.)),
ylim=c(0,1.1),xlab='Stratigraphic Height [m]',xlim=c(930,1007.5))
points(SeymourIslandBins$Witts.et.al.2016.Section.A,
rep(0.2,length(SeymourIslandBins$Witts.et.al.2016.Section.A)),pch=19)
lines(c(1007.5,1007.5),c(0,2))
legend('topleft',pch=c(1,19,NA),lty=c(NA,NA,1),legend=c('Macellari','Witts','K/Pg boundary'))
# }
Run the code above in your browser using DataLab