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$2*Nb_gpt;
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();

return (p);
}