Compute Sobol indices from a finished PEcAn run

Loads model outputs from a Sobol ensemble, calculates summary statistics for a chosen variable, feeds them to sensitivity::tell(), and returns the updated Sobol object.

Usage

compute_sobol_indices(outdir, sobol_obj, var = "GPP", stat_fun = mean)

Arguments

outdir

PEcAn run output directory that contains runs.txt

sobol_obj

object produced by PEcAn.uncertainty::generate_joint_ensemble_design()

var

Variable name to summarise (default "GPP").

stat_fun

Summary statistic applied to var default mean .

Value

sobol_obj .