# \donttest{
# Data simulation
set.seed(1)
simul <- SimulateRegression(n = 1000, pk = 50, family = "gaussian")
# Training/test split
ids <- Split(data = simul$ydata, tau = c(0.8, 0.2))
stab <- VariableSelection(
xdata = simul$xdata[ids[[1]], ],
ydata = simul$ydata[ids[[1]], ]
)
# Constructing the ensemble model
ensemble <- Ensemble(
stability = stab,
xdata = simul$xdata[ids[[1]], ],
ydata = simul$ydata[ids[[1]], ]
)
# Making predictions
yhat <- EnsemblePredictions(
ensemble = ensemble,
xdata = simul$xdata[ids[[2]], ]
)
# Calculating Q-squared
cor(simul$ydata[ids[[2]], ], yhat)^2
# }
Run the code above in your browser using DataLab