Odtnhj

I'm trying to use the sensitivity package to calculate sobol sensitivity indices from my parameters and model outputs. This function sobolSmthSpl in the package is supposed to do this but I just keep getting the error

 Error: $ operator is invalid for atomic vectors
 In addition: Warning message:
 In parallel::mclapply(1:ncol(Y), function(i) smooth.spline(X[, i], :
 all scheduled cores encountered errors in user code*

I've made the classes of the variables the same to the example, but just can't get rid of the error. Has anyone used this function before or have met a similar problem?

Here is the code of the function

 function (Y, X) 
 
 ANS = list()
 ANS[["call"]] = match.call()
 ANS[["X"]] = X
 ANS[["Y"]] = Y
 par.names = colnames(X)
 if (is.null(colnames(X))) 
 par.names = paste0("X", 1:ncol(X))
 X = normalize(X)
 Y = Y - mean(Y)
 Y = sapply(1:ncol(X), function(i) return(Y[order(X[, i])]))
 X = sapply(1:ncol(X), function(i) return(X[order(X[, i]), 
 i]))
 SMTH = optSmooth(Y, X, c(-2, 2))
 SA.tab = t(sapply(SMTH, est.Si))
 colnames(SA.tab) = c("Si", "se", "q0.05")
 rownames(SA.tab) = par.names
 ANS[["S"]] = SA.tab
 class(ANS) = "sobolSmthSpl"
 return(ANS)