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