I am currently fitting the following distributions using JMP and looking for ways to fit the same distributions in R: Zero Inflated Lognormal Zero Inflated Loglogistic Zero Inflated Frechet Zero Inflated Weibull Threshold Frechet Threshold Loglogistic Threshold Lognormal Log Generalized Gamma Threshold Weibull LEV Logistic Normal SEV Are there any packages that contain these distributions??? I am specifically interested in the zero inflated distributions since the data I have contains quite a bit of zeros. Thank you. Ravi -- View this message in context: http://r.789695.n4.nabble.com/Zero-Inflated-Distributions-tp3334861p3334861.html Sent from the R help mailing list archive at Nabble.com.
library(sos) findFn("Zero Inflated Lognormal") ---------------------------------------------------------------------------- ir. Thierry Onkelinx Instituut voor natuur- en bosonderzoek team Biometrie & Kwaliteitszorg Gaverstraat 4 9500 Geraardsbergen Belgium Research Institute for Nature and Forest team Biometrics & Quality Assurance Gaverstraat 4 9500 Geraardsbergen Belgium tel. + 32 54/436 185 Thierry.Onkelinx at inbo.be www.inbo.be To call in the statistician after the experiment is done may be no more than asking him to perform a post-mortem examination: he may be able to say what the experiment died of. ~ Sir Ronald Aylmer Fisher The plural of anecdote is not data. ~ Roger Brinner The combination of some data and an aching desire for an answer does not ensure that a reasonable answer can be extracted from a given body of data. ~ John Tukey> -----Oorspronkelijk bericht----- > Van: r-help-bounces at r-project.org > [mailto:r-help-bounces at r-project.org] Namens vioravis > Verzonden: vrijdag 4 maart 2011 9:40 > Aan: r-help at r-project.org > Onderwerp: [R] Zero Inflated Distributions > > I am currently fitting the following distributions using JMP > and looking for ways to fit the same distributions in R: > > Zero Inflated Lognormal > Zero Inflated Loglogistic > Zero Inflated Frechet > Zero Inflated Weibull > Threshold Frechet > Threshold Loglogistic > Threshold Lognormal > Log Generalized Gamma > Threshold Weibull > LEV > Logistic > Normal > SEV > > Are there any packages that contain these distributions??? I > am specifically interested in the zero inflated distributions > since the data I have contains quite a bit of zeros. > > Thank you. > > Ravi > > -- > View this message in context: > http://r.789695.n4.nabble.com/Zero-Inflated-Distributions-tp33 > 34861p3334861.html > Sent from the R help mailing list archive at Nabble.com. > > ______________________________________________ > 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. >
Thanks, Thierry. Has anyone used the "bayescount" for estimating zero inflated distributions? It states that it is a "crude function". Does that mean the estimates are only approximate??? The example they have given seems to work only with Gamma Poisson. data <- rpois(100, rgamma(100, shape=1, scale=8)) data[1:15] <- 0 maximise.likelihood(data, "ZIGP") However, when I tried fitting Gamma/LogNormal/Weibull (assuming that data is continuous), it throws out the following error: shape scale zi 9.532 4 21 Error in optim(c(shape, scale, zi), f6, control = list(fnscale = -1)) : function cannot be evaluated at initial parameters What is this error about??? Moreover, the function seems extremely slow. For the 100 data point example considered, it takes around 8 seconds for the estimation. Please let me know your opinions on this package and alternative packages, if any. Thank you. Ravi -- View this message in context: http://r.789695.n4.nabble.com/Zero-Inflated-Distributions-tp3334861p3335122.html Sent from the R help mailing list archive at Nabble.com.