Dear @paulremo, dear all,
it seems that I observed a further well-hidden bug, i.e. a bug that only shows up if only one parameter is unknown and also the error distribution is given. In this situation, in the report generated by
uq_print_uq_inversion.m the value for the 0.05 Quantile is also printed as value for the 0.95 Quantile, see the output
SCR_BI_to_2261_to_send-out.zip (2.0 KB) of
SCR_BI_to_2261_to_send.m (1.8 KB) .
According to my investigations, it seems that this incorrect output is created from the line 166 in
uq_print_uq_inversion.m. There it is excepted that
percentiles(1,mm) is the value for for the 0.05 Quantile for the mm-th unkown parameter/Error and
percentiles(2,mm) is the value for the corresponding 0.95 Quantile., which seems to satisfied in most situations, but not in the simple one under considerations.
percentiles is copied from
uq_print_uq_inversion.m the corresponding field is printed in my test file, showing that in the special situation that there is only one unknown parameter and no unknown error size it holds that one can find the value for the 0.95 Quantile in
module.Results.PostProc.Percentiles.Values (1,2) and not in
module.Results.PostProc.Percentiles.Values (2,1) as expected somehow by
The content of
module.Results.PostProc.Percentiles.Values is copied
percentiles containing the output of
uq_inversion_percentiles(Sample3D, perc_probabilities), which is an output of the matlab funtion
So, I suggest that somewhere during one of this coping of data processes it should be to checked if one is in the spacial situation that the (copied) output of
quantile is a (1,2) field and to create then a copy that is the transposed (2,1) field, such that the report generated by
uq_print_uq_inversion.m is also correct in simple situations.
To decide where this should be done requires some more knowledge on the UQLab programming style then I have.