Primitive Mantle

<data plot> Reference [wiki]: Data: bulk Input: ppm </data>

require('ggplot');

spdgrams.ree_primitive_mantle <- function(ppm_data)
{
primitive_mantle = data.frame(Cs = c(0.032), Rb = c(0.635), Ba = c(6.989), Th = c(0.085),
                              U = c(0.021), Nb = c(0.713), Ta = c(0.041), La = c(0.687),
                              Ce = c(1.775), Pb = c(0.185), Pr = c(0.276), Sr = c(21.1),
                              Nd = c(1.354), Zr = c(11.2), Hf = c(0.309), Sm = c(0.444),
                              Eu = c(0.168), Gd = c(0.596), Tb = c(0.108), Dy = c(0.737),
                              Y = c(4.55), Ho = c(0.164), Er = c(0.48), Tm = c(0.074),
                              Yb = c(0.493), Lu = c(0.074));

ns <- names(primitive_mantle);

ds <- selectNames(ppm_data,c(paste0(ns,'_ppm')));

names(ds) <- ns;

ret <- c();

values <- cbind(ppm_data$Name, rowApply(`/`,ds,primitive_mantle[1,]));
names(values)[1] <- 'Name';

ret$values <- values;
ret$grades <- ns;
ret$labels <- c('REE', 'Sample/Primitive mantle');

return(ret);
}