would be great if we could get CIs for the estimated parameters in fit_disc_gamma, using Hessian = TRUE in optim should do the job
would be great if we could get CIs for the estimated parameters in fit_disc_gamma, using Hessian = TRUE in optim should do the job