Diagramas ternários são muito comuns com dados geológicos, no R eles podem ser feitos de forma simples com o pacote Ternary (Smith 2017). Vamos usar os dados da tabela ‘ternario.csv’.
> library(Ternary)
> data<-read.csv("ternario.csv")
Primeiro fazemos um plot simples, sem os dados, podemos editá-lo de diversas formas:
> TernaryPlot(alab="Quartzo", blab="Feldspato", clab="Plagioclásio", atip="Q", btip="A", ctip="P", point="up", lab.cex=0.8, grid.minor.lines = 0, grid.lty="solid", col=rgb(0.9, 0.9, 0.9), grid.col="white", axis.col="grey", ticks.col="grey", padding=0.08)
> HorizontalGrid() #adiciona linhas de grid
Agora adcionamos os dados das amostras, nossa matriz de dados ‘data’ possui três rochas, veja como elas são plotadas separadamente:
> TernaryPoints(data[1:7,], cex=1.5, bg="yellow", pch=21)
> TernaryPoints(data[8:13,], cex=1.5, bg="blue", pch=21)
> TernaryPoints(data[14:21,], cex=1.5, bg="green", pch=21)
A função TernaryPolygon() desenha linhas conectando as amostras:
> TernaryPolygon(data[1:7,])
> TernaryPolygon(data[8:13,])
Por fim, vamos desenhar uma legenda para o gráfico usando a função legend():
> legend('topright', pch=21, pt.cex=1.5, pt.bg=c("yellow","blue","green"), legend=c('Rocha A', 'Rocha B', 'Rocha C'), 0.8, bty='n')