Basalts Geodinamic Settings Meschede

Basalts Geodinamic Settings Meschede.

require('ggtern');

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

ds$X <- ds$Zr_gpt/4;
ds$Y <- ds$Nb_gpt*10;
ds$Z <- ds$Y_gpt;

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

return (p);
}