====== Yb+Ta-Rb: Granitic Rocks, Pearce, 1984  ======

Yb+Ta-Rb 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 | Yb_ppm, Ta_ppm, Rb_ppm |

<use:classifydiag.granites_ybta_rb_pearce>
<test:classifydiag.granites_ybta_rb_pearce>

=== Source ===

<Code>
require('ggplot2');
 
classifydiag.granites_ybta_rb_pearce <- function (wtp_data, ...)
{  
ds <- selectNames(wtp_data,c('Name','Yb_ppm','Ta_ppm','Rb_ppm'))
 
pts <- data.frame(Name=ds$Name,
                  X = ds$Yb_ppm + ds$Ta_ppm,
                  Y = ds$Rb_ppm);
 
labels <- c("Yb + Ta", "Rb");
 
lines <- data.frame(X1 = c(  6, 200,   0,   40),
                    Y1 = c(  0, 450, 105, 2000),
                    X2 = c(  6,   6,   6,    6),
                    Y2 = c(200,   8, 200,  200));
 
marks  <- data.frame(Label=c("VAG", "ORG", "syn-COLG", "WPG"),
                         X=c( 1.05,    40,       1.15,    80),
                         Y=c(   5 ,     7,       1000,   700));

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

</Code>

{{tag>classifydiag wholerock}}