function update_multivariate_normal_inverse_gamma_gaussian(x:Real[_], ν:Real[_], Λ:LLT, γ:Real, α:Real, β:Real) -> (Real[_], LLT, Real, Real, Real)

Update the parameters of a normal inverse-gamma distribution with a multivariate Gaussian likelihood.

  • x: The variate.
  • ν: Prior precision times mean.
  • Λ: Prior precision.
  • γ: Prior scale accumulator.
  • α: Prior shape of the inverse-gamma.
  • β: Prior scale of the inverse-gamma.

Returns: the posterior hyperparameters μ', Λ', α' and β'.