Alk-MgO-FeO: Thoeliitic - Calc-Alkaline

Thoeliitic - Calc-Alkaline for bulk composition.

Использовать Проверить

require('ggtern');

classifydiag.amf <- function (wtp_data, ...)
{  
ds <- selectNames(wtp_data,c('Name','K2O_wtp','Na2O_wtp','MgO_wtp','FeO_wtp','Fe2O3_wtp'));

pts <- data.frame(Name=ds$Name,
                  X = ds$K2O_wtp + ds$Na2O_wtp,
                  Y = ds$FeO_wtp + ds$Fe2O3_wtp,
                  Z = ds$MgO_wtp);
 
labels <- c("Alk", "FeO*", "MgO");

lines <- data.frame(X1 = c(0.68, 0.63, 0.59, 0.46, 0.37, 0.30, 0.25, 0.22, 0.20, 0.19, 0.18, 0.17),
                    Y1 = c(0.29, 0.32, 0.35, 0.44, 0.50, 0.53, 0.55, 0.55, 0.55, 0.52, 0.47, 0.40),
                    Z1 = c(0.03, 0.05, 0.06, 0.10, 0.15, 0.17, 0.20, 0.23, 0.25, 0.29, 0.36, 0.43),
                    X2 = c(0.63, 0.59, 0.46, 0.37, 0.30, 0.25, 0.22, 0.20, 0.19, 0.18, 0.17, 0.15),
                    Y2 = c(0.32, 0.35, 0.44, 0.50, 0.53, 0.55, 0.55, 0.55, 0.52, 0.47, 0.40, 0.35),
                    Z2 = c(0.05, 0.06, 0.10, 0.15, 0.17, 0.20, 0.23, 0.25, 0.29, 0.36, 0.43, 0.50));

marks  <- data.frame(Label=c("Tholeiitic", "Calc-Alcaline"),
                     X=c(0.18, 0.57),
                     Y=c(0.60, 0.15),
                     Z=c(0.22, 0.28));

ret <- c();
ret$type <- "ternogram";
ret$points <- pts;
ret$labels  <- labels;
ret$lines <- lines;
ret$marks <- marks;
 
return(ret);
}