classifydiag:harker_si_k [pin.ifz.ru]

Эта страница только для чтения. Вы можете посмотреть её исходный текст, но не можете его изменить. Сообщите администратору, если считаете, что это неправильно.

====== SiO2 - K2O : Harker ======

Harker's diagramm SiO2 - K2O. Nomolized by oxides sum.

=== Details ===

| Reference | <ref:harker1909> |
| Data | bulk |
| Input | Na2O_wtp, MgO_wtp, Al2O3_wtp, SiO2_wtp, P2O5_wtp, K2O_wtp, CaO_wtp, TiO2_wtp, MnO_wtp, Fe2O3_wtp |

<use:classifydiag.harker_si_k>
<test:classifydiag.harker_si_k>

=== Source ===

<Code linenums lang-r>
classifydiag.harker_si_k <- function (wtp_data, ...)
{  
ds <- selectNames(wtp_data, c('Name','Na2O_wtp', 'MgO_wtp', 'Al2O3_wtp', 
                              'SiO2_wtp', 'P2O5_wtp', 'K2O_wtp', 'CaO_wtp',
                              'TiO2_wtp', 'MnO_wtp', 'Fe2O3_wtp'));

norm_sum <- rowSums(ds[,c('Na2O_wtp', 'MgO_wtp', 'Al2O3_wtp', 
                          'SiO2_wtp', 'P2O5_wtp', 'K2O_wtp', 'CaO_wtp',
                          'TiO2_wtp', 'MnO_wtp', 'Fe2O3_wtp')]);

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

labels <- c("SiO2", "K2O");

ret <- c();
ret$type <- "diagram";
ret$points <- pts;
ret$labels  <- labels;
ret$limits <- c(45,85,0,10);
 
return(ret);
}
</Code>
{{tag>classifydiag harker}}