Hi. I think there may be one or more zeros in your data set, causing the problem: x <- rgamma(100) fitdistr(x, "weibull") fitdistr(c(x,0), "weibull") Maybe you should omit the zeros. Christian