====== Alk-MgO-FeO: Thoeliitic - Calc-Alkaline ======

Thoeliitic - Calc-Alkaline for bulk composition.

=== Details ===

| Reference | | 
| Data | |
| Input | K2O_wtp, Na2O_wtp, MgO_wtp, FeO_wtp, Fe2O3_wtp |

<use:classifydiag.amf>
<test:classifydiag.amf>

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

=== Source ===

<Code linenums lang-r>
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);
}
</Code>

{{tag>classifydiag}}