## Composite charcoal record for North America:
ID=pfSiteSel(id_region=="WNA0",l12=1)
## Transform data
res3=pfTransform(ID,method=c("MinMax","Box-Cox","Z-Score"),BasePeriod=c(200,4000))
## Composite
comp=pfComposite(res3,bins=seq(0,5000,200))
plot(comp,type="density",smoothing=TRUE,spar=0.3)
Run the code above in your browser using DataLab