oce (version 0.8-10)

plotTaylor: Plot a model-data comparison diagram.

Description

Plot a model-data comparison diagram as described by Taylor (2001).

Usage

plotTaylor(x, y, scale, pch, col)

Arguments

x

a vector of reference values of some quantity, e.g. measured over time or space.

y

a matrix whose columns hold values of values to be compared with those in x. (If y is a vector, it is converted first to a one-column matrix).

scale

optional scale, interpreted as the maximum value of standard deviation.

pch

list of characters to plot, one for each column of y.

col

list of colours for points on the plot, one for each column of y.

Details

Creates a diagram as described by Taylor (2001). The graph is in the form of a semicircle, with radial lines and spokes connecting at a focus point on the flat (lower) edge. The radius of a point on the graph indicates the standard deviation of the corresponding quantity, i.e. x and the columns in y. The angle connecting a point on the graph to the focus provides an indication of correlation coefficient with respect to x. The ``east'' side of the graph indicates $R=1$, while $R=0$ is at the ``north edge'' and $R=-1$ is at the ``west'' side. The x data are indicated with a bullet on the graph, appearing on the lower edge to the right of the focus at a distance indicating the standard deviation of x. The other points on the graph represent the columns of y, coded automatically or with the supplied values of pch and col.

The example shows two tidal models of the Halifax sealevel data, computed with tidem with just the M2 component and the S2 component; the graph indicates that the M2 model is much better than the S2 model.

References

Taylor, Karl E., 2001. Summarizing multiple aspects of model performance in a single diagram, J. Geophys. Res., 106:D7, 7183--7192.

Examples

library(oce)
data(sealevelHalifax)
x <- sealevelHalifax[["elevation"]]
yM2 <- predict(tidem(sealevelHalifax, constituents="M2"))
yS2 <- predict(tidem(sealevelHalifax, constituents="S2"))
plotTaylor(x, cbind(yM2, yS2))