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