How is it possible to calculate the variance $\sigma^2$ for the Normal distribution if only $\alpha$ and $\beta$ (based on data) from the Inverse-gamma distribution are available? I followed the definitions on the Wiki and found out that for a Normal distribution in case of unknown variance and known mean the cojugate distribution is the Inverse-gamma with the updates $\alpha = \alpha + \frac{1}{2}$ and $\beta = \beta + \frac{(x - \mu)^2}{2}$.
I tried $\sigma^2 = \frac{\beta^2}{(\alpha - 1)^2 (\alpha - 2)} \text{for } \alpha > 2$ but $\alpha$ becomes less than $2$ after an update, and the variance returns NaN
.
Thanks for any suggestions.