The paper addresses a global sensitivity analysis of complex models. The work presents a generalization of the hierarchical statistical models where uncertain parameters determine the distribution of statistical models. The double randomization method is applied to increase the efficiency of the Monte Carlo estimation of Sobol indices. Numerical computations are provided to study the accuracy and efficiency of the proposed technique. The issue of optimization of the suggested approach is considered.
