Error when using fitdist function in R

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130409/005dabbc/attachment.pl>
An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130409/6969b077/attachment.pl>
Hi

where is function fitdist from? Do you by chance mean fitdistr from MASS package?
What is result of str(k)?

Without these questions answered you hardly get any useful answer.

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
PLEASE do read the posting guide http://www.R-project.org/posting-
guide.html
and provide commented, minimal, self-contained, reproducible code.
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
PLEASE do read the posting guide 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
PLEASE do read the posting guide 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/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595
Dear PIKAL,

Thank you for your kind reply. The function fitdist comes from package
fitdistrplus. I am attaching a .csv document that contains the data I want
to fit a distribution to, so you can have a good idea of what kind of data
I am dealing with.

Basically what str(k) contains is the information that appears in the
column Toll in the .csv document I am attaching. I will also show you the
result of str(k), for your reference:

 int [1:499] 363328 376216 367032 314826 311892 313340 367106 227742 386744
372002 ...

Hope this helps explain my situation better,

Best regards,

Paul

2013/4/10 PIKAL Petr <petr.pikal at precheza.cz>
Hi

where is function fitdist from? Do you by chance mean fitdistr from MASS
package?
What is result of str(k)?

Without these questions answered you hardly get any useful answer.

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
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.

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

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20130411/0f744e7a/attachment.pl>