FeO/MgO-SiO2: Tholeiitic-Calc-alkaline basalts

Tholeiitic-Calc-alkaline basalts.

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

require('ggplot2');

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

ds <- selectNames(wtp_data, c("Name", "SiO2_wtp", "MgO_wtp", "FeO_wtp", "Fe2O3_wtp"));

pts <- data.frame(Name=ds$Name,
                  X=ds$SiO2_wtp,
                  Y=(ds$FeO_wtp+ds$Fe2O3_wtp)/ds$MgO_wtp);

marks  <- data.frame(Label=c("Calc-alkaline", "Tholeiitic"),
                     X=c(55, 55),
                     Y=c(1, 3));

limits <- c(48,65,0,5);

lines <- data.frame(X1 = c(48),
                    Y1 = c(0.8),
                    X2 = c(65),
                    Y2 = c(3.5));

labels <- c("SiO2", "FeO*/MgO");

ret <- c();

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

return(ret);
}