If I was to do a 99% VaR calculation on a portfolio with normally distributed returns $\mathcal{N} (\mu,\sigma)$, the 99% VaR would be $\mu - 2.33\sigma$.
Instead of having a constant volatility, let's say volatility is gamma distributed, i.e. $\sigma \sim \Gamma(k, \theta)$.
Explain in as many details as possible (either derive a formula or explain a numerical solution using a computer program) how to compute the VaR of the portfolio when returns have a normal distribution conditional on sigma, and sigma is distributed according to a gamma distribution.