Skip to contents

Flow Frequency Sampler

Source: R/flow_frequency_sampler.R
flow_frequency_sampler.Rd

Generates a stratified matrix of flow values from a single set of frequency distribution parameters using stratified Monte Carlo sampling. Used internally by rfa_simulate() for median-only and full uncertainty modes.

Usage

flow_frequency_sampler(
  bestfit_params,
  freq_dist = "LP3",
  strat_dist = "ev1",
  Nbin = NULL,
  Mevent = NULL
)

Arguments

bestfit_params

Numeric vector of length 3 containing distribution parameters. For LP3: c(mean_log, sd_log, skew_log). For GEV: c(location, scale, shape).

freq_dist

Character. Distribution type. Either "LP3" (default) or "GEV".

strat_dist

Character. Probability space for stratification bins. Passed to stratified_sampler(). One of "ev1" (default), "normal", or "uniform". See stratified_sampler() for details.

Nbin

Integer. Number of stratified bins. Default is 50.

Mevent

Integer. Number of events per bin. Default is 200.

Value

A list containing:

flow

Matrix of sampled flow values [Mevent x Nbin]

nbins

Number of stratified bins

mevents

Number of events per bin

weights

Probability weights for each bin from stratified_sampler()

See also

stratified_sampler(), qp3(), rfa_simulate()

Examples

# Single parameter set (posterior mode)
params <- c(4.85, 0.39, -0.15)
result <- flow_frequency_sampler(params, freq_dist = "LP3",
                                 Nbin = 20, Mevent = 500)
dim(result$flow)  # 500 x 20
#> [1] 500  20
hist(result$flow)