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.pi...@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-boun...@r-project.org [mailto:r-help-bounces@r- > > project.org] On Behalf Of Paul Bernal > > Sent: Tuesday, April 09, 2013 11:19 PM > > To: r-help@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@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. >
______________________________________________ R-help@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.
