search for: fitdistrplus

Dear Charles, Please, when you have questions about fitdistrplus, contact directly the authors of the package and not R-help. When fitting non ? standard ? distributions with fitdistrplus, you should define by yourself the density and the cumulative distribution functions, or load a package which define them. See FAQ for a general example : https://cran.r-proj...

Dear useRs, We are pleased to announce you a new version of "fitdistrplus" on CRAN : https://cran.r-project.org/web/packages/fitdistrplus/index.html "fitdistrplus" is a package dedicated to help the fit of a parametric distribution to non-censored or censored data. The main new features in this release are few new topics in the FAQ vignette (https://c...

Dear useRs, We are pleased to announce you a new version of "fitdistrplus" on CRAN : https://cran.r-project.org/web/packages/fitdistrplus/index.html "fitdistrplus" is a package dedicated to help the fit of a parametric distribution to non-censored or censored data. The main new features in this release are few new topics in the FAQ vignette (https://c...

We are pleased to announce your a new version of fitdistrplus ( https://cran.r-project.org/package=fitdistrplus). Among the new features of the package (https://cran.r-project.org/web/packages/fitdistrplus/NEWS), a FAQ vignette is now available ( https://cran.r-project.org/web/packages/fitdistrplus/vignettes/FAQ.html). We will be delighted to update it with...

We are pleased to announce your a new version of fitdistrplus ( https://cran.r-project.org/package=fitdistrplus). Among the new features of the package (https://cran.r-project.org/web/packages/fitdistrplus/NEWS), a FAQ vignette is now available ( https://cran.r-project.org/web/packages/fitdistrplus/vignettes/FAQ.html). We will be delighted to update it with...

Dear Sirs, The below listed code fits a gamma and a pareto distribution to a data set danishuni. However the distributions are not appropriate to fit both tails of the data set hence a mixed distribution is required which has ben defined as "mixgampar" as shown below. library(fitdistrplus) x<- danishuni$Loss fgam<- fitdist(x,"gamma",lower=0) fpar<- fitdist(x,"pareto",start = list(shape=2,scale=2),lower=0) fmixgampar<- fitdist(x,"mixgampar",start = list(prob=1/2,nu=1,lambda=1,alpha=1,theta=1),lower=0) Error in fitdist(x, "mixgampar&quo...

Dear All, Apologies for not providing a reproducible example, but if I could, then I would be able to answer myself my question. Essentially, I am trying to fit a very complicated custom probability distribution to some data. Fitdistrplus does in principle everything which I need, but if require me to specify not only the density function d, but also the cumulative p and and inverse cumulative function q (see for instance http://www.stat.umn.edu/geyer/old/5101/rlook.html to understand what these quantities are in the case of a nor...

...<lorenzo.isella at gmail.com> wrote: > Dear All, > Apologies for not providing a reproducible example, but if I could, then I > would be able to answer myself my question. > Essentially, I am trying to fit a very complicated custom probability > distribution to some data. > Fitdistrplus does in principle everything which I need, but if require me > to specify not only the density function d, but also the cumulative p and > and inverse cumulative function q (see for instance > > http://www.stat.umn.edu/geyer/old/5101/rlook.html > > to understand what these quantit...

Hello, I like to fit data against a negative binominal distribution x2<-c(rep(10,14),rep(9,8),rep(8,13),rep(7,11),rep(6,6),rep(5,18),rep(4,7),re p(3,21),rep(2,33),rep(1,55),rep(0,225)) f2<-fitdist(x2,"nbinom",method="mle") plot(f2) summary(f2) gofstat(f2) I receive the following result: Fitting of the distribution ' nbinom ' by maximum

...a variety of chi-squared fitting techniques to fit a linear background and Gaussian directly to the binned counts to obtain a position and area under a peak. I understand that this is poor practice, especially given my rather low counting rate, so I have been attempting to fit the raw data with the fitdistrplus package (linked to from http://www.r-bloggers.com/fitting-distributions-with-r). I can fit an isolated Gaussian peak successfully using these methods (the central peak in my attached image). However, I cannot find a suitable distribution to account for the background in the peaks to the left of th...

Dear all, I have been trying to fit my data (only right censored) with gumbel distribution using fitdistrplus. I am getting very high standard error. I have been wondering why. The followings are the outputs: fit1=fitdistcens(dr0, "gumbel", start=list(a=99, b=0.6), optim.method= "L-BFGS-B", lower = 0.0, upper = Inf) > summary(fit1) FITTING OF THE DISTRIBUTION ' gumbel ' BY M...

Hi, Let say I have data by two columns A and B, and I have fit each column using the gamma distribution by 'fitdist' . I just want the result show only the shape and rate only. Eg: library(fitdistrplus) A <-c(1,2,3,4,5) B<-c(6,7,8,9,10) C <-cbind(A,B) apply(C, 2, fitdist, "gamma") Output show like this: $A Fitting of the distribution ' gamma ' by maximum likelihood Parameters: estimate Std. Error shape 3.702253 2.2440052 rate 1.234126 0.8011369 $B Fitting...

Hi, sorry for repeating the question but this is kind of important to me and i don't know whom should i ask. So as noted before when I do a parameter fit to the beta distr i get: fitdist(vectNorm,"beta"); Fitting of the distribution ' beta ' by maximum likelihood Parameters: estimate Std. Error shape1 2.148779 0.1458042 shape2 810.067515 61.8608126 Warning

Hi there, I have fitted a sample (with size 20) to a normal and/or logistic distribution using fitdistr() in MASS or fitdist() in fitdistrplus package. It's easy to get the parameter estimates. Now, I hope to report the confidence interval for those parameter estimates. However, I don't find a function that could give the confidence interval in R. I hope to write a function, however, I don't find some detailed information...

Hi, I am trying to estiamte parameters for gamma distribution using mle for below data using fitdist & fitdistr functions which are from "fitdistrplus" & "MASS"packages . I am getting errors for both functions. Can someone please let me know how to overcome this issue?? data y1<- c(256656, 76376, 6467673, 46446, 3400, 3100, 5760, 4562, 8000, 512, 4545, 4562, 4645, 4521, 4624, 5110, 4678, 4782, 4587, 491100,48426, 457940, 45...

Hi all, I am trying to fit a distribution to some data about survival times. I am interested only in a specific interval, e.g., while the data lies in the interval (0,...., 600), I want the best for the interval (0,..., 24). I have tried both fitdistr (MASS package) and fitdist (from the fitdistrplus package), but I could not get them working, e.g. fitdistr(left, "weibull", upper=24) Error in optim(x = c(529L, 528L, 527L, 526L, 525L, 524L, 523L, 522L, 521L, : L-BFGS-B needs finite values of 'fn' In addition: Warning message: In dweibull(x, shape, scale, log) : NaNs produc...

Hi all, I'm currently working on the fitdistrplus package (that basically fit distributions). There is something I do not understand about the generic function summary. In the current version on CRAN, there is no NAMESPACE saying S3method(summary, fitdist) . However if we use summary on an object send by fitdist function it works fine......

Hi, I am trying to fit a Johnson SB distribution using fitdist function in fitdistrplus Library. I have defined the Johnson SB distribution from ( http://www.ntrand.com/johnson-sb-distribution/) . But it gives me the follwing errors. Any help would be appreciated #xi = xi #lambda =l #delta =d #gamma = g djohn = function(x,xi,l,d,g) (d/(l*sqrt(2*pi)*((x-xi)/l)*(1-((x-xi)/l))))*exp[-0...

Hi, I am trying fit certain data into Beta distribution. I get the error saying "Error in fitdist(discrete_random_variable_c, "beta", start = NULL, fix.arg = NULL) : the function mle failed to estimate the parameters, with the error code 100" Below is the sorted data that I am trying to fit. Where am I going wrong. Thanks a lot for any help. Vinod

...ibution. I have also had a look at the 'Mixtools' and 'Flexmix' packages. With both of these it is possible to fit a mixture of distributions, but, as far as I can tell, these are pre-defined and do not have a function for a mixture of two Weibull distributions. I have used the 'fitdistrplus' package to fit the single distribution. Does anyone have any suggestions either for getting the log-likelihood for models fitted with the mix() function, or other packages / functions which allow me to fit a double Weibull and get the log-likelihood for the model? Any suggestions are much ap...