Классификация слюд в системе FeAl/MgLi

Классификация по статье True and brittle micas: composition and solid-solution series. G. TISCHENDORF, H.-J. FORSTER2, B. GOTTESMANN, M. RIEDER

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

classifydiag.micas_tischendorf <- function (wtp_data, annotate=T, labels=F, ...)
{  
font_size=3;
line_height=0.7;
ds <- selectNames(recalc.ms_tpf_realc0(wtp_data),c('Mg_pfu','Li_pfu','Fe_pfu','Al_c4_pfu'));

end_members <- data.frame(x=c(0,-1,-2,0,0,2,3,0), y=c(3,2,-1,-3,-1.8,-1,0,1), lab=c('ann','','pol','mus','hyp-mus','eas','phl','sid'));

p <- ggplot(ds, aes(x=Mg_pfu-Li_pfu,y=Fe_pfu-Al_c4_pfu)) +
  geom_point() +
  geom_point(data=end_members,aes(x=x,y=y),shape = 25) +
  geom_text(data=end_members,aes(x=x,y=y,label=lab), size=font_size, lineheight=line_height, hjust=-1, vjust=0) +
  theme_bw() +
  xlim(c(-2,3)) +
  ylim(c(-3,3)) + 
  xlab("Li <-> Mg")+
  ylab("Al_iv <-> Fe_tot");
  


if (labels)
{
    p <- p + geom_text(aes(label=Name),hjust=0, vjust=0, size=font_size, lineheight=line_height);
}
  
  
return(p);
}