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@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-Gaussia... _____________________________________ Sent from http://r-br.2285057.n4.nabble.com