When I try to compute the Viterbi_EM_fit() in an HMM model using DiagonalGaussian prior I get the error:
n, muhat, sumsq = self._get_statistics(data) ValueError: too many values to unpack
the parameters seem to be inverted in the max_likelihood function.