[R-br] [Dúvida] Máxima Verossimilhança no R.

[Rubem Kaipper Ceratti]

Pedro,


neste caso, vc precisa usar method="L-BFGS-B".


Att.,

Rubem



________________________________
 De: Pedro Rafael <pedro.rafael.marinho em gmail.com>
Para: r-br em 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))
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120321/678816d4/attachment.html>