estimate_delay function

Estimate a delay distribution from a reporting triangle