SiO2 - P2O5 : Harker

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

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

classifydiag.harker_si_p <- 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$P2O5_wtp*100/norm_sum);
 
labels <- c("SiO2", "P2O5");
 
ret <- c();
ret$type <- "diagram";
ret$points <- pts;
ret$labels  <- labels;
ret$lines <- data.frame();
ret$marks <- data.frame();
 
return(ret);
}