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

Benilton Carvalho beniltoncarvalho em gmail.com
Quarta Março 21 20:03:30 BRT 2012


Informe o intervalo desejado usando os argumentos lower e upper da
funcao optim. b

2012/3/21 Pedro Rafael <pedro.rafael.marinho em gmail.com>:
> Pessoal estou com um problema para estimar os parâmetros de dados com
> distribuição Beta com parâmetros p e q. Gerei aleatoriamente um
> conjuntos de dados com distribuição densidade de probabilidade beta no
> intervalo (0,1). Encontrei a distribuição de log verossimilança e como
> não existe solução fechada para os parâmetros p e q, temos que estimar
> por métodos numéricos como o de Newton Raphson, BFGS, etc.
>
> Código:
>
> ########################################
> 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((1-p)*log(y)+(1-q)*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="BFGS")
>
> Como informar o espaço paramétrico de cada um dos parâmetros da
> distribuição, no caso particular da distribuição beta ambos sendo
> maiores que zero?
>
> --
> Saudações,
> Pedro Rafael Diniz Marinho.
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.


Mais detalhes sobre a lista de discussão R-br