Compute Stage-Frequency Curve

Converts a matrix of routed peak stages into a stage-frequency curve using the law of total probability across stratified bins. For each candidate stage, the conditional exceedance probability within each bin is weighted by the bin probability and summed to produce the unconditional annual exceedance probability (AEP).

Usage

stage_frequency_curve(peakStage, weights, stage_bins = 1000)

Arguments

peakStage

Matrix of peak stages [Mevents x Nbins] from routing.

weights

Numeric vector of bin weights from stratified_sampler(). Must have length equal to ncol(peakStage).

stage_bins

Numeric value of exceedance stages used to compute stage-frequency in each bin. Default is 1000.

Value

A data frame with columns:

stage

Evaluated stage values

AEP

Annual exceedance probability at each stage

See also

stratified_sampler(), flow_frequency_sampler(), rfa_simulate()