We present a multivariate GARCH (MGARCH) model to forecast model averaged (co-)variances, correlations, and means in behavioral data for single individuals. We consider four MGARCH model parameterizations; three classic models (CCC, DCC, andBEKK) as well as a recently introduced model (pdBEKK) optimized for behavioral research. To obtain the averaged forecasts across the four models, we will first need to compute the model weights, obtained from a model stacking method. To do so, we need to approximate the expected log predictive density via a fully Bayesian leave-future-out cross validation technique. This approach has not been described so far in the literature. We provide an illustrative implementation using real data on two individuals from the longitudinal Intelligent Systems for Assessing Aging Change (ISAAC) study covering up to 4 years of daily measurements on individual computer use and walking speed. The individual participants show distinct patterns in the model weights, suggesting that individuals differ in the parameterizations that best capture their behavior. We generate weighted forecasts for up to 5 consecutive two-week periods for both individuals. These foreceast are compared to forecasts from the single best model. The resulting predictions are shown to be superior or at least equivalent to the forecasts from the single best model. We close with a discussion on limitations and an outlook. We provide an R package