Basalts Geodinamic Settings Pearce

Basalts Geodinamic Settings Pearce.

require('ggtern');

classifydiag.basalts_pearce <- function (wtp_data, annotate=T, labels=F, ...)
{  
font_size=3;
line_height=0.7;
ds <- selectNames(wtp_data,c('Name','Ti_gpt','Zr_gpt','Y_gpt'));

ds$X <- ds$Zr_gpt;
ds$Y <- ds$Ti_gpt/100;
ds$Z <- ds$Y_gpt*3;

p <- ggtern(ds, aes(x=X, y=Y, z=Z)) +
  geom_point() +
  labs( x = "Zr_gpt", y = "Ti_gpt/100", z = "Y_gpt*3" ) + 
  theme_bw();
  
  if (labels)
{
    p <- p + geom_text(aes(label=Name),hjust=0, vjust=0, size=font_size, lineheight=line_height);
}

return (p);
}