NbZrY {GCDkit}R Documentation

Meschede (1986) Zr/4 - 2Nb - Y LDA/QDA)

Description

Linear and quadratic discrimination-analysis-based versions of the ternary diagram Zr/4 - 2Nb - Y (Meschede 1986)) for geotectonic discrimination of igneous rocks of Vermeesch 2006.

Usage

NbZrY(Zr = WR[, "Zr", drop = FALSE], Nb = WR[, "Nb", drop = FALSE], 
    Y = WR[, "Y", drop = FALSE], type = c("Meschede", "LDA", "QDA"), ternary = TRUE, 
    show.labels = TRUE, short = TRUE, ...)

Arguments

Zr

vector of Zr concentrations (ppm)

Nb

vector of Nb concentrations (ppm)

Y

vector of Y concentrations (ppm)

type

either 'LDA' for linear discriminant analysis, 'QDA' for quadratic discriminant analysis, or 'Meschede' for the nominal decision boundaries of Meschede (1986).

ternary

logical. If FALSE, produces a logratio plot. Only used if type is LDA or QDA

show.labels

logical. If TRUE, labels the discrimination fields

short

logical. uses abbreviated labels for the discrimination fields

...

additional arguments for the generic points function.

Details

This plot comes from the GeoplotR collection available online at https://isoplotr.es.ucl.ac.uk/geoplotr.

NbZrY.png

Value

A list with tectonic affinities.

Note

This plot is not Figaro-compatible.

Author(s)

Pieter Vermeesch, p.vermeesch@ucl.ac.uk,

GCDkit implementation by Vojtěch Janoušek, vojtech.janousek@geology.cz

References

Meschede M (1986) A method of discriminating between different types of mid-ocean ridge basalts and continental tholeiites with the Nb-Zr-Y diagram. Chem Geol 56:207-218 doi: 10.1016/0009-2541(86)90004-5

Vermeesch P (2006) Tectonic discrimination diagrams revisited. Geochem Geophys Geosyst 7:Q06017 doi: 10.1029/2005GC001092

See Also

Meschede

Examples

sampleDataset("atacazo")

out <- NbZrY(type="QDA")
print(out$class)

print(out)

[Package GCDkit version 6.3.0 Index]