R: Rank histogram for ensemble forecasts

Rankhist {SpecsVerification} | R Documentation

## Rank histogram for ensemble forecasts

### Description

Calculate the rank histogram for an archive of ensemble forecasts and their corresponding verifying observations.

### Usage

    Rankhist(ens, obs, reduce.bins = 1, handle.na = "na.fail")

### Arguments

`ens` | matrix of dimension (N,K). An archive of K-member ensemble forecasts for N time instances.
`obs` | vector of length N. The corresponding verifying observations.
`reduce.bins` | number of adjacent bins that will be merged into one bin; has to be a divisor of K+1
`handle.na` | how should missing values in ensemble and observation data be handled; possible values are 'na.fail' (fails if any data is missing) and 'use.complete' (only uses times where all ensemble members and obs are available); default: 'na.fail'

### Value

a vector of length (K+1)/reduce.bins containing the rank counts

### References

Anderson J.L. (1996). A Method for Producing and Evaluating Probabilistic Forecasts from Ensemble Model Integrations. J. Climate, 9, 1518–1530. Hammill T.M. (2001). Interpretation of Rank Histograms for Verifying Ensemble Forecasts. Mon. Wea. Rev., 129, 550–560.

### See Also

PlotRankhist, TestRankhist

### Examples

    data(eurotempforecast)
    rh <- Rankhist(ens, obs)

[Package SpecsVerification version 0.5-3 Index]