# Define settings
N <- 2^9; params <- paste("X", 1:8, sep = "")
# Create sample matrix
mat <- sobol_matrices(N = N, params = params, matrices = "A")
# Compute the Sobol' G function
Y <- sobol_Fun(mat)
# Compute the S-ersatz discrepancy values
ind <- discrepancy_ersatz(mat = mat, Y = Y, params = params)
Run the code above in your browser using DataLab