A explicação do script a seguir está no artigo R na Prática – Parte 26.
setwd("C:\\geoStats\\dados\\compositional")
cdata30=read.csv("cdata30.csv",sep=";",
header=TRUE)
library(compositions)
cdata=acomp(cdata30)
#setwd("C:\\geoStats\\figuras\\chap3")
#jpeg("Fig3_15.jpeg",width=5,height=5,
#units="in",res=300)
barplot(cdata,col=c("darkblue","aquamarine4",
"floralwhite"))
#dev.off()