[R] Error when using fitdist function in R

[Paul Bernal]

Dear Dr. Brian,

Hope you are doing well. The thing is that my data does not have any 0s in
it. Furthermore, I also tried fitting other distributions (i.e. weibull)
and still did not work.

I will show you the results of str(k) so you have an idea of the tipe of
data I am dealing with:

 int [1:499] 363328 376216 367032 314826 311892 313340 367106 227742 386744
372002 ...
I will also attach, the .csv file that contains the data I am working with,
for your reference.

Any help you can provide me with will be greatly appreciated,

Best regards and thanks for your kind reply,

Paul



2013/4/10 Prof Brian Ripley <ripley at stats.ox.ac.uk>

> On 10/04/2013 08:29, PIKAL Petr wrote:
>
>> Hi
>>
>> where is function fitdist from? Do you by chance mean fitdistr from MASS
>> package?
>>
>
> No, from the output.  More likely
>
> fitdistrplus::fitdist   Fit of univariate distributions to non-censored
>                         data
>
>
> What is result of str(k)?
>>
>> Without these questions answered you hardly get any useful answer.
>>
>
> Also, the data seem to be integers: a gamma distribution is for real
> positive values only.  Having zero values would give the results you see.
>
>
>
>> Regards
>> Petr
>>
>>
>>
>> -----Original Message-----
>>> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-
>>> project.org] On Behalf Of Paul Bernal
>>> Sent: Tuesday, April 09, 2013 11:19 PM
>>> To: r-help at r-project.org
>>> Subject: [R] Error when using fitdist function in R
>>>
>>> Hello everyone,
>>>
>>> I  was trying to do some distribution fitting with a numerical field
>>> called Tolls. The sample size = 999 rows.
>>>
>>> Basically I assigned the Toll data to a new variable K by doing:
>>> k<-dtest$Toll
>>>
>>> After that, tried to fit a gamma distribution by doing: fitG<-
>>> fitdist(k,
>>> "gamma")
>>>
>>> Then the following messages showed (oh and I checked for empty rows
>>> before doing this):
>>>
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> Warning in dgamma(c(363328L, 376216L, 367032L, 314826L, 311892L,
>>> 313340L,  :
>>>    NaNs produced
>>> [1] "Error in optim(par = vstart, fn = fnobj, fix.arg = fix.arg, obs =
>>> data,  : \n  non-finite finite-difference value [2]\n"
>>> attr(,"class")
>>> [1] "try-error"
>>> attr(,"condition")
>>> <simpleError in optim(par = vstart, fn = fnobj, fix.arg = fix.arg, obs
>>> =
>>> data,     ddistnam = ddistname, hessian = TRUE, method = meth, lower =
>>> lower,     upper = upper, ...): non-finite finite-difference value [2]>
>>> Error in fitdist(k, "gamma") :
>>>    the function mle failed to estimate the parameters,
>>>                  with the error code 100
>>>
>>>>
>>>>
>>>         [[alternative HTML version deleted]]
>>>
>>> ______________________________**________________
>>> R-help at r-project.org mailing list
>>> https://stat.ethz.ch/mailman/**listinfo/r-help<https://stat.ethz.ch/mailman/listinfo/r-help>
>>> PLEASE do read the posting guide http://www.R-project.org/**posting-<http://www.r-project.org/posting->
>>> guide.html
>>> and provide commented, minimal, self-contained, reproducible code.
>>>
>>
>> ______________________________**________________
>> R-help at r-project.org mailing list
>> https://stat.ethz.ch/mailman/**listinfo/r-help<https://stat.ethz.ch/mailman/listinfo/r-help>
>> PLEASE do read the posting guide http://www.R-project.org/**
>> posting-guide.html <http://www.r-project.org/posting-guide.html>
>> and provide commented, minimal, self-contained, reproducible code.
>>
>>
>
> --
> Brian D. Ripley,                  ripley at stats.ox.ac.uk
> Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~**ripley/<http://www.stats.ox.ac.uk/~ripley/>
> University of Oxford,             Tel:  +44 1865 272861 (self)
> 1 South Parks Road,                     +44 1865 272866 (PA)
> Oxford OX1 3TG, UK                Fax:  +44 1865 272595
>