Dear rehena, You can use optim() with method="Brent", or optimize() directly, for your case.
But for you case, a popular approach is Newton-Raphson root finding algorithm. You can use newton function under my package OOmisc. You can check whether the global maximum or local maximum is found by replacing the solution in the second derivative of the log-likelihood function and if its value at the solution is less than 0, you can be sure that global maximum is secured. Hope this helps Ozgur -- View this message in context: http://r.789695.n4.nabble.com/Need-help-in-using-OPTIM-OPTIMIZATION-function-tp4633796p4633921.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.