Q-J diagram Morimoto and Kitamura, 1983

Classifiaction of pyroxene in Q-J system.

require('ggplot2');
source('recalc_cpx_sturm.r');

classifydiag.cpx_morimoto_kitamura <- function (wtp_data, annotate=T ...)
{  
ds <- selectNames(recalc.cpx_sturm(wtp_data),c('Name','Ca_pfu','Mg_pfu','Fe_M2_pfu','Na_pfu'));

ds$J <- 2*ds$Na_pfu;
ds$Q <- ds$Ca_pfu + ds$Fe_M2_pfu + ds$Mg_pfu;

txts <- data.frame(Label=c('Quad','Ca-Na','Na','Other','Jo,Ka,\nPe,Es','En,Fs,Df,\nHd,Wo'),
                     J=c(0.1,0.9,1.5,0.3,0,0.0),
                     Q=c(1.5,0.9,0.2,0.5,1,2));

p <- ggplot(ds, aes(x=J,y=Q)) +
  geom_point() +
  xlim(c(0,2)) +
  ylim(c(0,2)) + 
  geom_abline(intercept = 0, slope = 0.25) +
  geom_abline(intercept = 0, slope = 4) +
  geom_abline(intercept = 1.5, slope = -1) +
  geom_abline(intercept = 2.0, slope = -1) +
  xlab("J")+
  ylab("Q");
  
  if (annotate)
  {
  p <- p + geom_text(data=txts, aes(label=Label));
  }
  
  return(p);
}