====== Continental Curst: Middle, Rudnik and Gao, 2003 ======

=== Details ===

| Reference | <ref:rudnick2003> |
| Data | bulk |
| Input | ppm |

<use:spdgrams.ree_continental_curst_middle>
<test:spdgrams.ree_continental_curst_middle>

<WRAP center round todo 60%>
  * <del>заголовок</del>
  * краткое описание
  * <del>статья</del>
  * <del>ссылка</del>
  * <del>карточка</del>
  * подробное описание
  * <del>код</del>
  * тест
</WRAP>


=== Source ===

<Code linenums lang-r>
require('ggplot2');

cmp <- read.csv("constants_continental_curst_rudnick2003.csv", sep=";", stringsAsFactors=F);

spdgrams.ree_continental_curst_middle <- function(ppm_data, el_order = NULL)
{
ns <- el_order;
 
if (is.null(el_order))
{
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)$middle);

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, Middle');

return(ret);
}
</Code>

{{tag>spdgrams continental}}