data("obs")
data("sample_vars_sda")
obs$Cr_ppm <- log(obs$Cr_ppm)
krm <- rmvoutlier(obs$Cr_ppm)
y <- obs$Cr_ppm[-krm]
x <- list(sample_vars_sda[[1]][-krm, 1:11])
system.time({ # ~0.01s
sx <- selectvarsda(y, xlist = x)
})
Run the code above in your browser using DataLab