====== K2O-SiO2: Basic-Acid basalts ======

Basic-Acid basalts

=== Details ===

| Reference | |
| Data | bulk |
| Input | SiO2_wtp, K2O_wtp |

<use:classifydiag.basalts_k2o_sio2>
<test:classifydiag.basalts_k2o_sio2>

<WRAP center round todo 60%>
  * <del>заголовок</del>
  * краткое описание
  * статья
  * ссылка
  * <del>карточка</del>
  * подробное описание
  * <del>код</del>
  * тест
</WRAP>

=== Код метода ===

<Code linenums lang-r>
require('ggplot2');

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

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

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

marks  <- data.frame(Label=c("Low-K", "Medium-K", "High-K", "Acid", "Basic"),
                     X=c(60, 60, 60, 60, 55),
                     Y=c(0.2, 1.2, 2.5, 3.8, 3.2));

limits <- c(50,65,0,4);

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");

ret <- c();

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

return(ret);
}
</Code>

{{tag>classifydiag wholerock}}