Estava tentando rodar vosso exemplo e não consegui, não reconhece a função fitdistr. Existe alguma bibliotéca para ela, que não foi eventualmente mencionada? Muito obrigado Leonard pelo esclarecimento. Em 26 de abril de 2011 18:16, Leonard Assis <assis.leonard@gmail.com> escreveu: Acredito q seja o erro padrão, repare este exemplo, direto do help da função fitdistr set.seed(123) x <- rgamma(100, shape = 5, rate = 0.1) fitdistr(x, "gamma") ## now do this directly with more control. fitdistr(x, dgamma, list(shape = 1, rate = 0.1), lower = 0.001) ainda no help, observe o que diz: " An object of class "fitdistr", a list with four components, estimate the parameter estimates, sd the estimated standard errors, vcov the estimated variance-covariance matrix, and loglik the log-likelihood. " executando zzz <-fitdistr(x, dgamma, list(shape = 1, rate = 0.1), lower = 0.001) zzz$sd repare que retorna os valores que estão entre parêntesis Leonard