[R] Error when using fitdist function in R
Paul Bernal
paulbernal07 at gmail.com
Wed Apr 10 15:09:35 CEST 2013
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
>
More information about the R-help
mailing list