Hey Guillaume, I have a question, @ 10:40 , I see you use the "Method of Moments" to get the estimator for the initial value (iv = c(a.mm, I.mm) ), did you have any article to support this method? Because I am searching around to get the support for my study, and I can't find it... You are the only one have the same approach as me to be honest...