On Oct 19, 2013, at 2:26 AM, Charles Thuo wrote:
require(MASS)
n<- 100000
n
[1] 1e+05
w<- rgamma(n,rate=0.5,shape=5);head(w)
[1] 14.135679 18.839336 6.623502 17.165505 2.847572 7.511726
fit.w<- fitdistr(w,"weibull")
Warning messages:
1: In densfun(x, parm[1], parm[2], ...) : NaNs produced
2: In densfun(x, parm[1], parm[2], ...) : NaNs produced
3: In densfun(x, parm[1], parm[2], ...) : NaNs produced
4: In densfun(x, parm[1], parm[2], ...) : NaNs produced
5: In densfun(x, parm[1], parm[2], ...) : NaNs produced
6: In densfun(x, parm[1], parm[2], ...) : NaNs produced
7: In densfun(x, parm[1], parm[2], ...) : NaNs produced
8: In densfun(x, parm[1], parm[2], ...) : NaNs produced
summary(fit.w)
Length Class Mode
estimate 2 -none- numeric
sd 2 -none- numeric
vcov 4 -none- numeric
loglik 1 -none- numeric
n 1 -none- numeric
I do not understand why the summary does not have the fitted
parameters.
methods(summary)
There is no summary.fitdistr function so you are just seeing
summary.default's method of summarizing a list.
--
David Winsemius, MD
Alameda, CA, USA
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
