Suppose we know $\mu_1,\mu_2, \dots,\mu_n $ which are true values. Suppose we simulate data $x_{i1},\dots, x_{in}$ from some model or distribution involving, $\mu_i$. Using this we want to calculate the Mean Squared Error(MSE), thinking that the mean is a good estimate for $\mu_i$.
$$\hat{MSE}=\frac{1}{n}\sum^n_{i=1} (\mu_i-\bar{x_i})^2$$
My question is What is the variance of the estimate of this MSE and how we calculate it, analytically and computationally? I wonder whether the central limit theorem can be applied here.
Computationally, is it wise to just simulate the entire experiment $M$ times and compute it's variance. i.e. For simulation $M=1$, generate $x_{i1},\dots, x_{in}$ for each $i$,compute $\bar{x}_i$ for each $i$, compute $MSE$, repeat from start. We obtain $M$ different MSE estimates and so can compute the sample variance of that. This would be known as the variance or monte-carlo variance.