Vfrdtyky

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)

    }