Re: [R] Optimizing nested function with nlminb()

On 22-11-2012, at 22:55, Дмитрий Островский wrote: > I am trying to optimize custom likelyhood with nlminb() > Arguments h and f are meant to be fixed. > > example.R: > > compute.hyper.log.likelyhood <- function(a, h, f) { > a1 <- a[1] > a2 <- a[2] > l <- 0.0 > for (j in 1:length(f)) { >

[R] Optimizing nested function with nlminb()

I am trying to optimize custom likelyhood with nlminb() Arguments h and f are meant to be fixed. example.R: compute.hyper.log.likelyhood <- function(a, h, f) { a1 <- a[1] a2 <- a[2] l <- 0.0 for (j in 1:length(f)) { l <- l + lbeta(a1 + f[j], a2 + h - f[j]) - lbeta(a1, a2) } return