# NOT RUN {
data("copals_camera")
Y <- scale(Y[, 1:50])
X <- model.matrix(~ group + 0)
residuals <- lm(as.matrix(Y) ~ X - 1)$residuals
S12_inv <- whitening(residuals, "AR1", pAR = 1, qMA = 0)
Frequencies <- variable_selection(
Y = Y, X = X,
square_root_inv_hat_Sigma = S12_inv,
nb_repli = 10, nb.cores = 1, parallel = FALSE
)
# }
# NOT RUN {
# Parallel computing
require(doMC)
registerDoMC(cores=4)
Freqs <- variable_selection(Y,X,square_root_inv_hat_Sigma,
nb_repli=10,parallel=TRUE,nb.cores=4)
# }
Run the code above in your browser using DataLab