Error code 100 when using the function “fitdist” from the fitdistrplus package
Hi Nelly, The message David suggested was about scaling the values, not adjusting the parameters. It is quite possible that the empirical distribution is nothing like beta or Weibull. Have you tried plotting the values with "density"? Jim
On Fri, Aug 5, 2016 at 6:56 AM, Nelly Reduan <nell.redu at hotmail.fr> wrote:
Hello, Thanks David for your answer. I tested different parameters for Beta and Weibull distributions based on forums but I didn't find solutions. Thank you very much for your time. Have a nice day Nell
________________________________ De : David Winsemius <dwinsemius at comcast.net> Envoy? : jeudi 4 ao?t 2016 08:45:01 ? : Nelly Reduan Cc : r-help at r-project.org Objet : Re: [R] Error code 100 when using the function ?fitdist? from the fitdistrplus package On Aug 3, 2016, at 4:42 PM, Nelly Reduan <nell.redu at hotmail.fr> wrote: 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> Cannot get anything useful from that URL. <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 A bit of searching produces this: http://markmail.org/search/?q=list%3Aorg.r-project.r-help+fitdistrplus+error+code+100#query:list%3Aorg.r-project.r-help%20fitdistrplus%20error%20code%20100+page:1+mid:esp2okcorrdichtj+state:results -- David. 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]] ______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see 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. David Winsemius Alameda, CA, USA [[alternative HTML version deleted]] ______________________________________________ R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.