similar to: Fitting Mixed Distributions in the fitdistrplus package

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

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

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

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 new questions sent by users. Do not

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 new questions sent by users. Do not

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

Created corresponding unit test, and the optimization is bit exact with C function. --- silk/SigProc_FIX.h | 7 ++- silk/arm/arm_silk_map.c | 12 ++++ silk/arm/inner_prod_aligned_arm.h | 58 +++++++++++++++++++ silk/arm/inner_prod_aligned_neon_intr.c | 66 ++++++++++++++++++++++ silk/enc_API.c

Why not define your own functions based on d? e.g. myCumDist <- function(x) { integrate(d, lower=-Inf, upper=x)$value } myQuantile <- function(x) { uniroot(f=function(y) { h(y) - x }, interval=c(-5,5)) } # limits -5,5 should be replaced by your own which might require some fiddling e.g. d <- function(x) { exp(-x^2/2)/(sqrt(2*pi)) } # just an example for you to test with; use your own

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

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,

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... According to R-lang, we have " The most

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.5*(g +

Capture the results of the apply command into an object and then work with that. Here is one way to do it: > res <- apply(C, 2, fitdist, "gamma") > out <- c( res$A$estimate["shape"], res$B$estimate["shape"], res$A$estimate["rate"], res$B$estimate["rate"]) > names(out) <- c("A shape","B shape","A

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

Dear R-helpers, I am trying to fit my data to a 4-parameter lognormal distribution (aka Johnson Sb dist) with fitdist function from the library(fitdistrplus). So far, I have learnt that with "mle" method it's not always possible to estimate the gamma and delta parameters even if the bounding estimates are "known"/"guessed". Therefore, I tried to fit it with the

I would have thought this to be relatively elementary, but I can't find it mentioned in any of my stats texts. Please consider the following: library(fitdistrplus) fp = fitdist(y,"exp"); rate = fp$estimate; sd = fp$sd fOneWeek = exp(-rate*7); #fraction that happens within a week - y is measured in days fr = exp(-rate*dt); #fraction remaining - dt = elapsed time from

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,

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)

Previously, by mistake I added request in wish list for slide like MSN. Actully it was Blind like msn. However it was my mistake so I tried to created script by my slef - and able to do it. But I am new to this, so my script is not in good condition. Effect.Scale2 = Class.create(); Object.extend(Object.extend(Effect.Scale2.prototype, Effect.Base.prototype), { initialize: function(element,

I'm trying to estimate two parameter for a markov chain (exponential). I used maximum likelihood to do it. (fitdistr{MASS},fitdist{fitdistrplus}), but I get just one parameter. Do you have any suggestion to help me? THANKS!! -- View this message in context: http://r.789695.n4.nabble.com/Estimate-two-parameter-for-exponential-dist-tp4648917.html Sent from the R help mailing list archive at