Continental Curst, Rudnik and Gao, 2004 (Cs-Rb-Ba...)

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

require('ggplot');
 
cmp <- read.csv("constants_continental_curst_compositions.csv", sep=";", stringsAsFactors=F);
 
spdgrams.ree_continental_curst <- function(ppm_data)
{
ns <- c('Cs', 'Rb', 'Ba', 'Th', 'U', 'Nb', 'Ta', 'La', 'Ce', 'Pb', 'Pr',
        'Sr', 'Nd', 'Zr', 'Hf', 'Sm', 'Eu', 'Gd', 'Tb', 'Dy', 'Y', 'Ho',
        'Er', 'Tm', 'Yb', 'Lu');
 
pmsc <- c();
 
for (n in ns)
  pmsc <- c(pmsc, subset(cmp, component == n)$total);
 
ds <- selectNames(ppm_data,c(paste0(ns,'_ppm')));
 
names(ds) <- ns;
 
ret <- c();
 
values <- cbind(ppm_data$Name, rowApply(`/`,ds,pmsc));
names(values)[1] <- 'Name';
 
ret$values <- values;
ret$grades <- ns;
ret$labels <- c('REE', 'Sample/Continental curst, Total');
 
return(ret);
}