Calculates the likelihood of observing a set of discrete HI titres given a corresponding set of predicted titres
r_likelihood( expected, data, theta, expected_indices = NULL, measurement_shifts = NULL )
| expected | vector of expected HI titres |
|---|---|
| data | vector of observed discrete HI titres |
| expected_indices | the indices of the measurement_shifts vector that each predicted titre needs adding to it |
| measurement_shifts | the vector of measurement shifts for each cluster to add to the predicted titres |
a vector with the likelihood of making each observation given the predictions