hi how are you?
i use sobol2002 for sensitivity analysis which have two mesures a first-order sensitivity index and a total sensitivity index. but i need anly inorder to measure a sensitivity of four rainfall runoff models namely, IHACRES_CWI ,IHACRES_CMD GR4J and AWBM. when i use
```{r'}
# a gain of samples from the parameter sets
n <- 1000
# Find the parameter ranges from your model
# Make sure you define the right model!
X1 <- hydromad::parameterSets(getFreeParsRanges(dmody),n)
X2 <- hydromad::parameterSets(getFreeParsRanges(dmody),n)
# run the sensitivity
cm <- sobol2002(model = evalPars,
# using these parameter ranges
X1 = X1, X2 = X2,
# relates to the boot strapping in the analysis
nboot = 100,
# The model to use
object=dmody,
# The objective function to use
objective=~hmadstat("r.squared")(Q,X))
```
both measures are desplay
```{r}
print(cm)
```
Call:
sobol2002(model = evalPars, X1 = X1, X2 = X2, nboot = 100, object = wmod, objective = ~hmadstat("r.squared")(Q, X))
Model runs: 7000
First order indices:
original bias std. error min. c.i. max. c.i.
tw 0.28315064 0.0039998922 0.10168815 0.05889264 0.47335592
f 0.03454740 0.0017048688 0.02792947 -0.03010949 0.08094929
tau_s 0.03240779 0.0095650946 0.06490069 -0.11337796 0.15394195
v_s 0.36057461 -0.0010769543 0.08716641 0.20435177 0.55149404
tau_q 0.06478449 0.0006725721 0.03769932 -0.01995087 0.13772147
Total indices:
original bias std. error min. c.i. max. c.i.
tw 0.343931601 -0.007191064 0.09325466 0.174822883 0.53162276
f -0.007130158 -0.001698604 0.02951521 -0.059275505 0.06529394
tau_s 0.192021305 -0.005224342 0.06186857 0.063911156 0.31227424
v_s 0.529243827 0.011284257 0.08655102 0.317865031 0.69202420
tau_q 0.082840863 -0.001366429 0.04424592 0.001479365 0.18021432
the plot also display both main and total effects
but i need only total effect how to correct?