I have some data and several experimental models with 3 free parameters. I would like to use aicbic to compare the models, and to do this, I need maximum likelihood estimates. I am trying to use the mle function, but it only seems to generate parameter estimates and not the actual likelihood values. Is there a fast way I can compute maximum likelihood values?
I am using a custom probability density function, by the way.
Other solutions not using mle are also appreciated! Thanks!