similar to: Error in fitdist- mle failed to estimate parameters

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

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

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

Hi, Well, i need some help, practical and theoretical. I am wondering why the fitdistplus (mle function) is returning an error for this code: [code] x1 <- c(100,200,140,98,97,56,42,10,2,2,1,4,3,2,12,3,1,1,1,1,0,0); plotdist(x1); descdist(x1, boot =1000); y<- sum(x1); d= as.vector(length(x1)); for(i in 1:length(x1)){ d[i] = x1[i]/y; } fitdist(d, "beta") [/code] Error:

>From the documentation I have found, it seems that one of the functions from package plyr, or a combination of functions like split and lapply would allow me to have a really short R script to analyze all my data (I have reduced it to a couple hundred thousand records with about half a dozen records. I get the same result from ddply and split/lapply: >

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

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

The data.frame is constructed by one of the following functions: funweek <- function(df) if (length(df$elapsed_time) > 5) { rv = fitdist(df$elapsed_time,"exp") rv$year = df$sale_year[1] rv$sample = df$sale_week[1] rv$granularity = "week" rv } funmonth <- function(df) if (length(df$elapsed_time) > 5) { rv =

Dear all, I would like to use the following function fitdist(data, distr, method=c("mle", "mme", "qme", "mge"), start=NULL, fix.arg=NULL, ...) for many different distr values like distr=c("norm","lnorm","pois") (just a small example) and take back into a list the parameter name which is what is inside distr plus what the

I must have missed something simple, but still, i don't know what. I obtained my basic data as follows: x <- sprintf("SELECT m_id,sale_date,YEAR(sale_date) AS sale_year,WEEK(sale_date) AS sale_week,return_type,0.0001 + DATEDIFF(return_date,sale_date) AS elapsed_time FROM `merchants2`.`risk_input` WHERE DATEDIFF(return_date,sale_date) IS NOT NULL") moreinfo <- dbGetQuery(con,

Dear all, I need to plot an cumulative distribution plot of a variable and then to fit a distribution to that, probably a weibull or lognormal. I have plotted the ecdf as > plot(ecdf(x)) but I haven't managed to fit the distribution. I have as well attached the data. I would appreciate if you could help me on that. Thank you. Kind regards Maria -------------- next part --------------

Hello, I have detailed (with pictures and whatnot) my question on my weblog at http://www.cs-ed.org/blogs/mjadud/archives/2005/05/a_question_abou.html The short version of the question is this: When I ask 'fitdistr' to try and fit my distribution as a "weibull" distribution, it comes up with some rather wacky parameters. If I take the same distribution, and do something

Are there any functions in R for fitting distributions? In particular, I would like to fit the weibull and extreme-value distributions when: i) the data are given and the maximum likelihood estimates are required; ii) the sample moments are available but the complete data record is not (so moment estimates are required). Any suggestions will be gratefully received.

Hi R lovers Here is a numerical vector test > test [1] 206 53 124 112 92 77 118 75 48 176 90 74 107 126 99 84 114 147 99 114 99 84 99 99 99 99 99 104 1 159 100 53 [33] 132 82 85 106 136 99 110 82 99 99 89 107 99 68 130 99 99 110 99 95 153 93 136 51 103 95 99 72 99 50 110 37 [65] 102 104 92 90 94 99 76 81 109 91 98 96 104 104 93 99 125 89

Hi Folks! Just started using the gamlss package and I tried a simple code example (see below). Why the negative sigma? John > y <- rt(100, df=1)> m1<-fitDist(y, type="realline")Warning messages:1: In MLE(ll3, start = list(eta.mu = eta.mu, eta.sigma = eta.sigma, : possible convergence problem: optim gave code=1 false convergence (8)2: In MLE(ll4, start = list(eta.mu =