TAS Classification

TAS

require('ggplot2');

classifydiag.tas <- function (wtp_data, ...){

ds <- selectNames(wtp_data, c("Name", "SiO2_wtp", "K2O_wtp", "Na2O_wtp"));

pts <- data.frame(Name=ds$Name,
                  X=ds$SiO2_wtp,
                  Y=ds$K2O_wtp + ds$Na2O_wtp);

marks  <- data.frame(Label=c("Picro-\nbasalts",
                             "Foidite",
                             "Tephrite-\nBasanite",
                             "Basalts",
                             "Trachy-\nbasalts",
                             "Phono-\ntephrite",
                             "Basaltic\nandesite",
                             "Basaltic\ntrachy-\nandesite",
                             "Tephri-\nphonolite",
                             "Andesite",
                             "Trachy-\nandesite",
                             "Phonolite",
                             "Dacite",
                             "Trachy-\ndacite",
                             "Trachyte",
                             "Rhyolite"),
                     X=c( 43,  39,  45,  48,  49,  50,  55,  54,   54,  60,  58,   58,  67,  66,   65,  73),
                     Y=c(2.0, 8.5, 6.5, 3.0, 6.0, 9.0, 3.0, 6.0, 12.0, 4.0, 8.0, 14.0, 4.5, 8.0, 11.0, 8.0));

limits <- c(35,77,0,16);

lines <- data.frame(X1 = c(53,  57,  63, 53,  53,  53),
                    Y1 = c(0,   0,   0,  0.5, 1.6, 2.6),
                    X2 = c(53,  57,  63, 63,  63,  63),
                    Y2 = c(2.6, 3.15, 4,  1.0, 2.3, 4));

labels <- c("SiO2", "K2O+Na2O");

ret <- c();

ret$type <- "diagram";
ret$points <- pts;
ret$limits <- limits;
ret$marks <- marks;

ret$labels  <- labels;

return(ret);
}