reltest_simulate

Random training data from one model

Generates predictive distributions based on calibrating priors
for various commonly used statistical models, including models with
predictors. Routines for densities, probabilities, quantiles, random
deviates and the parameter posterior are provided. The predictions are
generated from the Bayesian prediction integral, with priors chosen to
give good reliability (also known as calibration). For homogeneous models,
the prior is set to the right Haar prior, giving predictions which are
exactly reliable. As a result, in repeated testing, the frequencies of
out-of-sample outcomes and the probabilities from the predictions agree.
For other models, the prior is chosen to give good reliability. Where
possible, the Bayesian prediction integral is solved exactly. Where exact
solutions are not possible, the Bayesian prediction integral is solved
using the Datta-Mukerjee-Ghosh-Sweeting (DMGS) asymptotic expansion.
Optionally, the prediction integral can also be solved using posterior
samples generated using Paul Northrop's ratio of uniforms sampling package
('rust'). Results are also generated based on maximum likelihood, for
comparison purposes. Various model selection diagnostics and testing
routines are included. Based on "Reducing reliability bias in assessments
of extreme weather risk using calibrating priors", Jewson, S., Sweeting, T.
and Jewson, L. (2024); <doi:10.5194/ascmo-11-1-2025>.

Stephen Jewson

fitdistcp

Distribution Fitting with Calibrating Priors for Commonly Used
Distributions

reltest_simulate function

<dl><dt>model</dt>
<dd>which distribution to test. Possibles values are
"<code>exp</code>",
"<code>pareto_k2</code>",
"<code>halfnorm</code>",
"<code>unif</code>",
"<code>norm</code>",
"<code>norm_dmgs</code>",
"<code>gnorm_k3</code>",
"<code>lnorm</code>",
"<code>lnorm_dmgs</code>",
"<code>logis</code>",
"<code>lst_k3</code>",
"<code>cauchy</code>",
"<code>gumbel</code>",
"<code>frechet_k1</code>",
"<code>weibull</code>",
"<code>gev_k3</code>",
"<code>exp_p1</code>",
"<code>pareto_p1k2</code>",
"<code>norm_p1</code>",
"<code>lnorm_p1</code>",
"<code>logis_p1</code>",
"<code>lst_p1k3</code>",
"<code>cauchy_p1</code>",
"<code>gumbel_p1</code>",
"<code>frechet_p2k1</code>",
"<code>weibull_p2</code>",
"<code>gev_p1k3</code>",
"<code>norm_p12</code>",
"<code>lst_p12k3</code>",
"<code>gamma</code>",
"<code>invgamma</code>",
"<code>invgauss</code>",
"<code>gev</code>",
"<code>gpd_k1</code>",
"<code>gev_p1</code>".
"<code>gev_p12</code>".
"<code>gev_p123</code>".</dd>
<dt>nx</dt>
<dd>the length of the training data to use.</dd>
<dt>tt</dt>
<dd>predictor vector</dd>
<dt>tt1</dt>
<dd>predictor vector 1</dd>
<dt>tt2</dt>
<dd>predictor vector 2</dd>
<dt>tt3</dt>
<dd>predictor vector 2</dd>
<dt>params</dt>
<dd>values for the parameters for the specified distribution</dd>
<dt>minxi</dt>
<dd>minimum value for EVT shape parameter</dd>
<dt>maxxi</dt>
<dd>maximum value for EVT shape parameter</dd></dl>

Arguments

Random training data from one model — reltest_simulate

<dl>

<dt>model</dt>
<dd>which distribution to test. Possibles values are
"<code>exp</code>",
"<code>pareto_k2</code>",
"<code>halfnorm</code>",
"<code>unif</code>",
"<code>norm</code>",
"<code>norm_dmgs</code>",
"<code>gnorm_k3</code>",
"<code>lnorm</code>",
"<code>lnorm_dmgs</code>",
"<code>logis</code>",
"<code>lst_k3</code>",
"<code>cauchy</code>",
"<code>gumbel</code>",
"<code>frechet_k1</code>",
"<code>weibull</code>",
"<code>gev_k3</code>",
"<code>exp_p1</code>",
"<code>pareto_p1k2</code>",
"<code>norm_p1</code>",
"<code>lnorm_p1</code>",
"<code>logis_p1</code>",
"<code>lst_p1k3</code>",
"<code>cauchy_p1</code>",
"<code>gumbel_p1</code>",
"<code>frechet_p2k1</code>",
"<code>weibull_p2</code>",
"<code>gev_p1k3</code>",
"<code>norm_p12</code>",
"<code>lst_p12k3</code>",
"<code>gamma</code>",
"<code>invgamma</code>",
"<code>invgauss</code>",
"<code>gev</code>",
"<code>gpd_k1</code>",
"<code>gev_p1</code>".
"<code>gev_p12</code>".
"<code>gev_p123</code>".</dd>


<dt>nx</dt>
<dd>the length of the training data to use.</dd>


<dt>tt</dt>
<dd>predictor vector</dd>


<dt>tt1</dt>
<dd>predictor vector 1</dd>


<dt>tt2</dt>
<dd>predictor vector 2</dd>


<dt>tt3</dt>
<dd>predictor vector 2</dd>


<dt>params</dt>
<dd>values for the parameters for the specified distribution</dd>


<dt>minxi</dt>
<dd>minimum value for EVT shape parameter</dd>


<dt>maxxi</dt>
<dd>maximum value for EVT shape parameter</dd>

</dl>

reltest_simulate: Random training data from one model

Description

Usage

Value

Arguments