| attribute {verification} | R Documentation |
An attribute plot illustrates the reliability, resolution and uncertainty of a forecast with respect to the observation. The frequency of binned forecast probabilities are plotted against proportions of binned observations. A perfect forecast would be indicated by a line plotted along the 1:1 line. Uncertainty is described as the vertical distance between this point and the 1:1 line. The relative frequency for each forecast value is displayed in parenthesis.
## Default S3 method:
attribute(x, obar.i, prob.y, obar = NULL, titl =
NULL, ...)
## S3 method for class 'prob.bin':
attribute(x, ...)
x |
A vector of forecast probabilities or a ``prob.bin''
class object produced by the verify function. |
obar.i |
A vector of observed relative frequency of forecast bins. |
prob.y |
Relative frequency of forecasts of forecast bins. |
obar |
Climatological or sample mean of observed events. |
titl |
Plot title. |
... |
Graphical parameters |
Matt Pocernich <pocernic@rap.ucar.edu>
Hsu, W. R., and A.H. Murphy, 1986: The attributes diagram: A geometrical framework for assessing the quality of probability forecasts. Int. J. Forecasting 2, 285-293.
Wilks, D. S. (1995) Statistical Methods in the Atmospheric Sciences Chapter 7, San Diego: Academic Press.
## Data from Wilks, table 7.3 page 246.
y.i <- c(0,0.05, seq(0.1, 1, 0.1))
obar.i <- c(0.006, 0.019, 0.059, 0.15, 0.277, 0.377, 0.511,
0.587, 0.723, 0.779, 0.934, 0.933)
prob.y<- c(0.4112, 0.0671, 0.1833, 0.0986, 0.0616, 0.0366,
0.0303, 0.0275, 0.245, 0.022, 0.017, 0.203)
obar<- 0.162
attribute(y.i, obar.i, prob.y, obar, titl = "Sample Attribute Plot")
## Function will work with a ``prob.bin'' class objects as well.
## Note this is a random forecast.
obs<- round(runif(100))
pred<- runif(100)
A<- verify(obs, pred, frcst.type = "prob", obs.type = "binary")
attribute(A, titl = "Alternative plot", xlab = "Alternate x label" )