data(Namib)
HMcomponents <- c("zr","tm","rt","TiOx","sph","ap","ep",
"gt","st","amp","cpx","opx")
am <- amalgamate(Namib$PTHM,feldspars=c("KF","P"),
lithics=c("Lm","Lv","Ls"),heavies=HMcomponents)
plot(ternary(am))
Run the code above in your browser using DataLab