# Define settings
N <- 1000; params <- paste("X", 1:3, sep = ""); R <- 10
# Create sample matrix
mat <- sobol_matrices(N = N, params = params)
# Compute Ishigami function
Y <- ishigami_Fun(mat)
# Compute and bootstrap Sobol' indices
ind <- sobol_indices(Y = Y, N = N, params = params, boot = TRUE, R = R)
# Plot Sobol' indices
plot(ind)
Run the code above in your browser using DataLab