N-MORB

<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.007), Rb = c(0.56), Ba = c(6.3), Th = c(0.12),
                              U = c(0.047), Nb = c(2.33), Ta = c(0.132), La = c(2.5),
                              Ce = c(7.5), Pb = c(0.3), Pr = c(1.32), Sr = c(90),
                              Nd = c(7.3), Zr = c(74), Hf = c(2.05), Sm = c(2.63),
                              Eu = c(1.02), Gd = c(3.68), Tb = c(0.67), Dy = c(4.55),
                              Y = c(28), Ho = c(1.01), Er = c(2.97), Tm = c(0.456),
                              Yb = c(3.05), Lu = c(0.455));

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);
}