@@ -15,12 +15,16 @@ authors:
1515 affiliation : 1
1616 corresponding : true
1717 - name : Woodrow L. Fields
18+ orcid : 0009-0008-7454-5552
1819 affiliation : 1
1920 - name : Julian Gonzalez
21+ orcid : 0009-0009-9058-7653
2022 affiliation : 1
2123 - name : Sadie Niblett
24+ orcid : 0009-0008-8588-4816
2225 affiliation : 1
2326 - name : Brennan Beam
27+ orcid : 0009-0003-0515-3727
2428 affiliation : 2
2529 - name : Brian Skahill
2630 orcid : 0000-0002-2164-0301
@@ -92,26 +96,31 @@ The following example illustrates a bootstrap uncertainty analysis [@efron1993]
9296``` csharp
9397using Numerics .Distributions ;
9498
95- // Generate a synthetic flood record from a Log-Pearson Type III distribution
99+ // Generate a synthetic flood record
96100int n = 80 ;
97101int B = 10000 ;
98- var lp3 = new LogPearsonTypeIII (meanOfLog : 3 , standardDeviationOfLog : 0 . 5 , skewOfLog : 0 . 25 );
99- double [] sample = lp3 .GenerateRandomValues (n , seed : 12345 );
102+ var lp3 = new LogPearsonTypeIII (
103+ meanOfLog : 3 ,
104+ standardDeviationOfLog : 0 . 5 ,
105+ skewOfLog : 0 . 25 );
106+ double [] sample = lp3 .GenerateRandomValues (n , 12345 );
100107
101- // Estimate distribution parameters from the sample using L-moments
102- lp3 .Estimate (sample , ParameterEstimationMethod .MethodOfLinearMoments );
108+ // Estimate parameters using L-moments
109+ lp3 .Estimate (sample ,ParameterEstimationMethod .MethodOfLinearMoments );
103110
104111// Configure the bootstrap analysis
105- var boot = new BootstrapAnalysis (distribution : lp3 ,
106- estimationMethod : ParameterEstimationMethod .MethodOfLinearMoments ,
107- sampleSize : n ,
108- replications : B ,
109- seed : 12345 );
110-
111- // Define non-exceedance probabilities for quantile estimation
112- var probabilities = new double [] { 0 . 999 , 0 . 99 , 0 . 9 , 0 . 8 , 0 . 7 , 0 . 5 , 0 . 3 , 0 . 2 , 0 . 1 , 0 . 05 , 0 . 02 , 0 . 01 };
113-
114- // Compute 90% BCa confidence intervals for each quantile
112+ var boot = new BootstrapAnalysis (
113+ distribution : lp3 ,
114+ estimationMethod : ParameterEstimationMethod .MethodOfLinearMoments ,
115+ sampleSize : n ,
116+ replications : B ,
117+ seed : 12345 );
118+
119+ // Non-exceedance probabilities for quantile estimation
120+ var probabilities = new double [] { 0 . 999 , 0 . 99 , 0 . 9 , 0 . 8 , 0 . 7 , 0 . 5 , 0 . 3 ,
121+ 0 . 2 , 0 . 1 , 0 . 05 , 0 . 02 , 0 . 01 };
122+
123+ // Compute 90% BCa confidence intervals
115124var CIs = boot .BCaQuantileCI (sample , probabilities , alpha : 0 . 1 );
116125```
117126
0 commit comments