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