code: def underlying_likelihood(self, binary_outcomes, modelparams, expparams):
          
          original_mps = modelparams[..., self._orig_mps_slice]
          return self.underlying_model.likelihood(binary_outcomes, original_mps, expparams)
docs: Given outcomes hypothesized for the underlying model, returns the likelihood
          which which those outcomes occur.