projbioplot {GCDkit}R Documentation

Calculates ternary coordinates projected from biotite (and plots the ternary diagram).

Description

The function projbiocoords calculates the coordinates used to define the projection, whereas projbioplot calls projbiocoords and does the other operations needed to define a template that will be used by plotDiagram.

Usage

    projbiocoords(where=WR,add=FALSE)
	projbioplot(mins=FALSE,addWR=FALSE,ticks=FALSE,xmin=-2,xmax=1,ymin=-0.5,ymax=0.5)

Arguments

projbiocoords:

where

A data matrix containing whole rock analyses, to be projected. By default WR

add

Boolean. If TRUE, the results of the calculations will be added to WR as new columns with colnames = c("ms1","fsp","CaAl","bio")

projbioplot:

mins

Boolean. If TRUE, the composition of ideal minerals will be plotted on the diagram, namely "q","fsp","an","an50","cz","Ep","sill","opx","cpx","olv","grs-Gt","Gt","NaCrd","Crd","bio","ms","MgHbl","Edn" and "Pgs".

addWR

Boolean. If TRUE, the results of the calculations will be added to WR as new columns with colnames = c("ms1","fsp","CaAl","bio")

ticks

Boolean. If TRUE, tick marks will be added to the side of the diagram (similar to the ticks option in ternary )

xmin, xmax, ymin, ymax

Numeric. Bounds of the plotting area, in rectangular coordinates. Note that the vertical axis (left side) is at x=0; ms1 = 3 al + 2 (Na + K) is at x=0 and y=0.5, and CaAl = Ca + Al is at x=0 and y=-0.5); the fsp (=Al + (Na + K) ) is at x=√{3}/2 and y=0. Plagioclase an50 is at x=-√{3}/2 and y=0, and most points should fall to the right of it. The defaults are therefore sensible.

Details

The "projection from biotite" of Moyen et al. (2017) is based on four coordinates based on molar proportions:

ms1 = Al - Ca - NK

fsp = -2 Al -2 Ca + 3 NK -1/3 FM

CaAl = Ca

bio = 1/3 FM

where Al = molar Al

Ca = molar Ca

FM = molar Fe + molar Mg

NK = molar Na + molar K

molar proportions are calculated by millications.

They are primarily used for plotting the relevant diagram (ms1-fsp-CaAl, note that bio is dropped during plotting as this is projected from bio).

projbiocoords calculates the coordinates. It also has the side effect of (globally) affecting its results to results, where they are available to addResults, r2clipboard(results), HTMLTableResults, etc.

projbioplot is used purely for defining a Figaro template. The values of the four coordinates are calculated (and (globally) assigned to results). If called with plotDiagram, the diagram is plotted, which is the real use of the function.

If used as a plugin, a GUI function .projbioGUI is also supplied (and linked to the menu item), in charge of gathering the missing arguments and calling plotDiagram.

projbio.png

Value

For projbioplot, nothing. This function is meant only to be called via plotDiagram.

projbiocoords returns a matrix of 4 columns containing the new coordinates, (ms1,fsp,CaAl,bio).

Author(s)

Jean-François Moyen jfmoyen@gmail.com

References

Moyen JF, Laurent O, Chelle-Michou C, Couzinie S, Vanderhaeghe O, Zeh A, Villaros A, Gardien V, 2017. Collision vs. subduction-related magmatism: two contrasting ways of granite formation and implications for crustal growth. Lithos 277:154-177. doi: 10.1016/j.lithos.2016.09.018

and in particular the supplementary item SE4 "Multivariate statistics and projection for granitic rocks".

See Also

plotDiagram millications

Examples

    data(sazava)
    accessVar("sazava")
    
    projbiocoords()
    plotDiagram("projbioplot",FALSE,TRUE)

    plotDiagram("projbioplot",FALSE,TRUE,mins=TRUE,xmin=-2)

[Package GCDkit version 6.1 Index]