ThTaHf {GCDkit}R Documentation

Wood (1980) Th - Hf/3 - Ta LDA/QDA)

Description

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

Usage

ThTaHf(Th = WR[, "Th", drop = FALSE], Hf = WR[, "Hf", drop = FALSE], 
    Ta = WR[, "Ta", drop = FALSE], type = c("Wood", "LDA", "QDA"), ternary = TRUE, 
    show.labels = TRUE, short = TRUE, ...)

Arguments

Th

vector of Th concentrations (ppm)

Hf

vector of Hf concentrations (ppm)

Ta

vector of Y concentrations (ppm)

type

either 'LDA' for linear discriminant analysis, 'QDA' for quadratic discriminant analysis, or 'Wood' for the nominal decision boundaries of Wood (1980)

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

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

Wood DA (1980) The application of a Th-Hf-Ta diagram to problems of tectonomagmatic classification and to establishing the nature of crustal contamination of basaltic lavas of the British Tertiary volcanic province. Earth Planet Sci Lett 50: 11-30 doi: 10.1016/0012-821X(80)90116-8

See Also

Wood

Examples

sampleDataset("sazava")

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

print(out)

[Package GCDkit version 6.3.0 Index]