Различия

Здесь показаны различия между двумя версиями данной страницы.

--- classifydiag:basalts_k2o_sio2 [2016/03/04 14:29]
admin создано
+++ classifydiag:basalts_k2o_sio2 [2018/11/12 15:07] (текущий)
@@ Строка 1: / Строка 1: @@
 ====== 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}}