[R-br] Distribuições Gamma e Inversa Gaussiana

[mauricio.soares em globo.com]

Felipe, obrigado!
Eu preciso rodar algum pacote antes de executar essa função?

<quote author='R-br mailing list'>
Pode usar a função optim:

### Gamma
G <- function(par, x)
{
   sum(dgamma(x, shape = par[1], rate = par[2], log = TRUE))
}


emv.G   <- optim(par = c(1, 10), fn = G, x = dados,
                   ,control=list("fnscale"=-1), hessian=TRUE)
emv.G$par


##### standard error
sqrt(diag(solve(-emv.G$hessian)))



-- 
Atenciosamente,


==================================
Felipe E. Barletta Mendes
Estatístico(UFPR) - Conre3 9766-A
Mestrando em Bioestatística(UEM)
+55 (41)-92077191
+55 (41)-33287216
===================================

_______________________________________________
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 fornea cdigo
mnimo reproduzvel.
</quote>
Quoted from: 
http://r-br.2285057.n4.nabble.com/R-br-Distribuicoes-Gamma-e-Inversa-Gaussiana-tp4667391p4667392.html


_____________________________________
Sent from http://r-br.2285057.n4.nabble.com