====== Y-Nb: Granitic Rocks, Pearce, 1984  ======

Y-Nb discriminant diagrams for syn-collision granites (syn-COLG), volcanic arc granites (VAG), within plate granites (WPG) and ocean ridge granites (ORG).

=== Details ===

| Reference | <ref:pearce1984> |
| Data | bulk |
| Input | Nb_ppm, Y_ppm |

<use:classifydiag.granites_y_nb_pearce>
<test:classifydiag.granites_y_nb_pearce>

=== Source ===

<Code>

require('ggplot2');
 
classifydiag.granites_y_nb_pearce <- function (wtp_data, ...)
{  
ds <- selectNames(wtp_data,c('Name','Nb_ppm','Y_ppm'))
 
pts <- data.frame(Name=ds$Name,
                  X = ds$Y_ppm,
                  Y = ds$Nb_ppm);
 
labels <- c("Y", "Nb");
 
lines <- data.frame(X1 = c(40, 1000,    3, 500),
                    Y1 = c(1,    95,  500, 250),
                    X2 = c(50,   50,   50,  30),
                    Y2 = c(11,   11,   11,  25));
 
marks  <- data.frame(Label=c("ORG", "WPG", "VAG +\n syn-COLG"),
                         X=c(  300,    20,                  4),
                         Y=c(    3,    90,                 40));

ret <- c();
ret$type <- "diagram";
ret$points <- pts;
ret$labels  <- labels;
ret$marks  <- marks;
ret$lines  <- lines;
ret$limits <- c(1,1000,1,2000);
ret$scale <- c('log','log');
 
return(ret);
}

</Code>

{{tag>classifydiag wholerock}}