# Bootstrapping the data to set up a prior

I am using a Gaussian model with a conjugate Normal-Inverse-Wishart (NIW) prior, as described here. The advantage of this approach is that the marginal likelihood $p(y)$, which is what I am interested in, is available in closed form.

My problem is that the results seem to be dependent of the NIW hyper-parameters (I have no prior information), with some of the dangers being described here.

As an alternative, I am considering bootstrapping my data in order to obtain $m$ estimate of the mean and covariance. Then I could calculate the marginal likelihood:

$$p(y) \approx \frac{1}{m} \sum_{i=1}^m p(y|\hat{\mu}_i, \hat{\Sigma}_i).$$

Would this prior be an approximation to an Empirical Bayes prior, something else or just nonsense?

Thank you