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
)

Arguments

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

Value

a vector with the likelihood of making each observation given the predictions