diff --git a/src/sobol_sensitivity.jl b/src/sobol_sensitivity.jl index 154be8da..235781b4 100644 --- a/src/sobol_sensitivity.jl +++ b/src/sobol_sensitivity.jl @@ -219,20 +219,20 @@ function gsa_sobol_all_y_analysis( Eᵢs, [ ( - sum(fA .^ 2 + fAⁱ[k] .^ 2) ./ (2n) .- + sum(fA .^ 2 + fAⁱ[k] .^ 2) ./ (2n) .- (sum(fA + fAⁱ[k]) ./ (2n)) .^ 2 - ) * ( - 1.0 .- + ) * ( + 1.0 .- ( - 1 / n .* sum(fA .* fAⁱ[k]) + 1 / n .* sum(fA .* fAⁱ[k]) .- (1 / n .* sum((fA .+ fAⁱ[k]) ./ 2)) .^ 2 - ) ./ + ) ./ ( - 1 / n .* sum((fA .^ 2 .+ fAⁱ[k] .^ 2) ./ 2) - + 1 / n .* sum((fA .^ 2 .+ fAⁱ[k] .^ 2) ./ 2) - (1 / n .* sum((fA .+ fAⁱ[k]) ./ 2)) .^ 2 - ) ) + ) for k in 1:d ] ) @@ -299,19 +299,19 @@ function gsa_sobol_all_y_analysis( hcat, [ ( - sum(fA .^ 2 + fAⁱ[k] .^ 2, dims = 2) ./ (2n) .- + sum(fA .^ 2 + fAⁱ[k] .^ 2, dims = 2) ./ (2n) .- (sum(fA + fAⁱ[k], dims = 2) ./ (2n)) .^ 2 - ) .* ( - 1.0 .- + ) .* ( + 1.0 .- ( - 1 / n .* sum(fA .* fAⁱ[k], dims = 2) .- + 1 / n .* sum(fA .* fAⁱ[k], dims = 2) .- (1 / n * sum((fA .+ fAⁱ[k]) ./ 2, dims = 2)) .^ 2 - ) ./ + ) ./ ( - 1 / n .* sum((fA .^ 2 .+ fAⁱ[k] .^ 2) ./ 2, dims = 2) .- + 1 / n .* sum((fA .^ 2 .+ fAⁱ[k] .^ 2) ./ 2, dims = 2) .- (1 / n * sum((fA .+ fAⁱ[k]) ./ 2, dims = 2)) .^ 2 - ) - ) for k in 1:d + ) + ) for k in 1:d ] ) )