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))