barCodes(
  x,
  yrng = c(0, 0.03),
  width = 20,
  col = rgb(0, 0, 0, 25, maxColorValue = 255),
  border = NA,
  ...
)

## Not run: 
#Load EUROEVOL Data
data(euroevol)

#Subset Danish Dates
denmark <- subset(euroevol,Country=="Denmark")

#Calibrate and Bin
denmarkDates <- calibrate(x=denmark$C14Age,errors=denmark$C14SD) 
denmarkBins <- binPrep(sites=denmark$SiteID,ages=denmark$C14Age,h=200) #200 years bin size

#Compute median date for each bin
bm <- binMed(x=denmarkDates,bins=denmarkBins)

#Compute median date for each date
dm <- medCal(denmarkDates)

#Compute SPD 
denmarkSPD <- spd(x=denmarkDates,bins=denmarkBins,timeRange=c(10000,4000))

#Plot SPD and barCodes of median dates
plot(denmarkSPD,runm=200)
barCodes(dm,yrng=c(0,0.01)) 

#Plot SPD and barCodes of median bins in BC/AD
plot(denmarkSPD,runm=200,calendar="BCAD")
barCodes(BPtoBCAD(bm),yrng=c(0,0.01)) 

## End(Not run)