Compute the fitness of the given parameter set against observed data

compute_fitness(
  params,
  model,
  functions,
  observed,
  sequence_length,
  recombination_rate,
  mutation_rate = 0,
  engine = NULL,
  model_args = NULL,
  engine_args = NULL,
  statistics = FALSE
)