which includes handling of `numpy.nan` for missing data in observations
and predictions. This is documented in the `evalp` docstring and new
unit tests are added to check this new functionality.