Pedro, neste caso, vc precisa usar method="L-BFGS-B". Att., Rubem ________________________________ De: Pedro Rafael <pedro.rafael.marinho@gmail.com> Para: r-br@listas.c3sl.ufpr.br Enviadas: Quarta-feira, 21 de Março de 2012 20:48 Assunto: Re: [R-br] [Dúvida] Máxima Verossimilhança no R. Realmente Benilton. Existia um problema na minha função de log-verossimilhança. Correção: ######################################## dados = rbeta(1000000,1,3) # onde p = 1 e q = 3. hist(dados) # Cara dos dados. log.verossimilhanca <- function(teta,y){ p = teta[1] q = teta[2] logl <-sum((p-1)*log(y)+(q-1)*log(1-y)+log(gamma(p+q)) - log(gamma(q)) - log(gamma(p))) return(-logl) } optim(theta <- c(1,1),fn = log.verossimilhanca,y=dados,method="Nelder-Mead", lower = c(0.001,0.001), upper = c(40,40))