Maximum likelihood estimation in R
[This article was first published on Quantitative Finance Collector, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
Maximum likelihood estimation can be implemented like Quasi-maximum likelihood in Matlab, You can also write an R function which computes out the likelihood function. As always in R, this can be done in several different ways.Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
One issue is that of restrictions upon parameters. When the search algorithm is running, it may stumble upon nonsensical values – such as a sigma below 0 – and you do need to think about this. One traditional way to deal with this is to “transform the parameter space”. As an example, for all positive values of sigma, log(sigma) ranges from -infinity to +infinity. So it’s safe to do an unconstrained search using log(sigma) as the free parameter.
For detail about methodology and sample codes see http://www.mayin.org/ajayshah/KB/R/documents/mle/mle.html.
Tags – mle
Read the full post at Maximum likelihood estimation in R.
To leave a comment for the author, please follow the link and comment on their blog: Quantitative Finance Collector.
R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.