UQWorld

Morris sensitivity analysis with vectors as output

Hi there.
I’m trying to do a sensitivity analysis with 11-dimensional inputs, and the output is a 196-dimensional binary vector. The output is obtained by taking the 11-d inputs and calling external programs to run simulation.
What I want to do is to combine the Morris method with Hamming Distance calculation. That is, the Hamming Distance between two output vectors (again, 196-d binary vector) as a result of the input change. I don’t know if it’s possible.
If not, is there a way to run sensitivity analysis with UQLab when the output is a 196-dimensional binary vector?
Thank you very much.
Tianjiao

Hi Tianjiao,
I can’t tell if UQLab can perform SA with your configuration. But your idea sounds very interesting. My understanding is that with the Hamming Distance you will obtain one single value per input of the Morris statistic and not 196 ones. I’ve never seen such an example before.

By the way, there is another way to conduct the analysis which will provide a set of 11x192 sensitivity indices (which might be cumbersome but at least you can guess which variable is important for each element of the output vector, if this can be of interest). Compute the delta-importance measure by comparing x_i to x_i\vert y_t=0 or to x_i\vert y_t=1 for t=1,...,196 and i=1,...,11.
Regards,
Thierry

Dear Tianjiao

First of all, from a practical point of view, the vast majority of the methods available in UQLab work seamlessly when the model has a vector of outputs. You have nothing to do!
So if you apply the method of Morris to a model with 196 outputs, you’ll get 196 different (\mu,\,\sigma) diagrams. Same if you use Sobol’ indices: you’ll get a matrix of first-order indices, a matrix of total indices, etc.
Now this usually applies to continuous outputs. And I guess it does not necessarily make sense to perform a sensitivity analysis with each 0/1 output in your case. If so, the approach proposed by Thierry would work. Otherwise, you can compute the Hamming distance with the null vector (Hamming weight) and use it as a single scalar output.
Hope it is useful!
Best regards
Bruno