  Let $\alpha_n(\cdot)=P\bigl(X_{n+1}\in\cdot\mid X_1,\ldots,X_n\bigr)$ be the
predictive distributions of a sequence $(X_1,X_2,\ldots)$ of $p$-variate random
variables. Suppose $$\alpha_n=\mathcal{N}_p(M_n,Q_n)$$ where
$M_n=\frac{1}{n}\sum_{i=1}^nX_i$ and
$Q_n=\frac{1}{n}\sum_{i=1}^n(X_i-M_n)(X_i-M_n)^t$. Then, there is a random
probability measure $\alpha$ on $\mathbb{R}^p$ such that
$\alpha_n\rightarrow\alpha$ weakly a.s. If $p\in\{1,2\}$, one also obtains
$\lVert\alpha_n-\alpha\rVert\overset{a.s.}\longrightarrow 0$ where
$\lVert\cdot\rVert$ is total variation distance. Moreover, the convergence rate
of $\lVert\alpha_n-\alpha\rVert$ is arbitrarily close to $n^{-1/2}$. These
results (apart from the one regarding the convergence rate) still apply even if
$\alpha_n=\mathcal{L}_p(M_n,Q_n)$, where $\mathcal{L}_p$ belongs to a class of
distributions much larger than the normal. Finally, the asymptotic behavior of
copula-based predictive distributions (introduced in [13]) is investigated and
a numerical experiment is performed.
