Skip to contents

Minimal use of floodam.cba.sa

How to load data

path_input = system.file("extdata", package = "floodam.cba.sa")
input = load_data(case = "test", path_input = path_input)
#> Case 'test' successfully loaded.
#> Key informations are:
#>  - 'floodam.cba.sa' in version 1.1.0.0 is used with:
#>      - n=1000
#>      - rate=official
#>  - 'sensitivity' in version 1.31.0 is used.
#>  - input are taken from '/tmp/Rtmpgcx88f/temp_libpath1110f56104e92/floodam.cba.sa/extdata/test'
#>  - output are saved in '/tmp/RtmpCH7K31/test'

How to perform analysis

Compute nominal value

nominal = calculate_nominal(input)
#> Model EAD: 1.077
#> Model NPV: 26.97

Do uncertainty and sensitivity analysis

analysis = analyze_uncertainty(input)
#> Registered S3 method overwritten by 'sensitivity':
#>   method    from 
#>   print.src dplyr
#> Model EAD: 1.077
#> Model NPV: 26.97
#> Evaluating model_npv... elapsed time: 0.197 secs
#> sentivity analysis is perform using sensitity::sobolmartinez
#> Evaluating sobol_analysis... elapsed time: 2.4 secs
#> Evaluating analyze_graphic... elapsed time: 0.073 secs
#> Evaluating compress... elapsed time: 0.158 secs

Plot graphs from results

plot_uncertainty(input, analysis)

plot_sobol(input, analysis)

plot_sobol(input, analysis, "T")