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