docs: def rand(self, unique_pstates, random_state=None):
          """
          Randomize the POHMM parameters
          """
          self._init_pstates(unique_pstates)
          self._init_random(random_state=random_state)
          self._compute_marginals()
          return self