[R] Error code 100 when using the function “fitdist” from the fitdistrplus package

Nelly Reduan nell.redu at hotmail.fr
Thu Aug 4 01:42:05 CEST 2016


Nelly Reduan has shared OneDrive?files with you. To view them, click the link or image below.


<https://1drv.ms/u/s!Apkg2VlgfYyDgQo8BnoB_Ds4KXgR>
[https://bzmvxw.by3301.livefilestore.com/y3mor2T_TYssPck9iMngzQsuiM_z140uCxN_MOvDhcRAILsrvwtWQ8cMMtzUvEuWjFYFytobNNvH8TJGzIrV7tjUHcKnVG_E76ru3RJEpGNPM3v-gKVJgYLhZPb9gLcQST6h-N6UtMumzrF9xVxCaWwJjLJ7amvLtDs816OaSSttsA?width=200&height=150&cropmode=center]<https://1drv.ms/u/s!Apkg2VlgfYyDgQo8BnoB_Ds4KXgR>


<https://1drv.ms/u/s!Apkg2VlgfYyDgQmEe6HOjP3d_4e2>
[https://r1.res.office365.com/owa/prem/images/dc-jpg_20.png]<https://1drv.ms/u/s!Apkg2VlgfYyDgQmEe6HOjP3d_4e2>

Figure_1.jpeg<https://1drv.ms/u/s!Apkg2VlgfYyDgQmEe6HOjP3d_4e2>




Hello,

I’m trying to fit distributions to data. To do this, I used the function “fitdist” from the fitdistrplus package and I drew a Cullen and Frey graph (attached Figure 1). From this graph, I am attempting to fit different distributions: Beta, Gamma and Weibull. The function “fitdist” works with Gamma distribution from this code:

Fit.dist <- fitdist(x[x!=0], distr="gamma", method="mle",lower=c(0, 0),start=list(scale=1,shape=1))

However, with Beta and Weibull distributions, I obtain this error message:

  the function mle failed to estimate the parameters,

                with the error code 100



Here is my code to fit Beta and Weibull distributions:

fit.dist <- fitdist(x_scaled, distr="beta", method = "mle")

fit.dist <- fitdist(x, distr="weibull", method="mle", lower=c(0, 0))

For the Beta distribution, I transformed the variable to have values between 0 and 1 as follows:

x_scaled  <- (x-min(x))/max(x)

Here are some information about data (summary() and attached Figure 2)

> summary(x)

   Min. 1st Qu.  Median    Mean 3rd Qu.    Max.

    0.0    90.4   244.8   437.4   549.4  4904.0

 How can I fit Gamma and Weibull distributions to my data without having the error message ?

Thank you very much for your time.

Nell


	[[alternative HTML version deleted]]



More information about the R-help mailing list