void PoissonModel::mle() {
double n = suf()->n();
double s = suf()->sum();
if (n > 0)
set_lam(s / n);
else
set_lam(1.0);
}
void ExponentialModel::mle() {
double number_of_observations = suf()->n();
double sum_of_durations = suf()->sum();
set_lam(number_of_observations / sum_of_durations);
}
