Hi Shekhar, It looks from your error that you have values outside the range 0 to 1. The beta distribution is only defined between 0 and 1. Can you please post your data set (or at least a portion of it)?
On Wed, May 4, 2011 at 3:24 PM, Shekhar <shekhar2...@gmail.com> wrote: > Hi Steven, > Thanks for the quick reply. i have tried but > its giving me error--->Error in optim(x = c(38.1815173696765, > -12.7988197976440, -3.88212459045077, : > initial value in 'vmmin' is not finite > > > i have tried something like this: > > library(MASS) > x<-rnorm(n=100,mean=10,sd=20); > fitdistr(x,dbeta,start=list(shape1=1,shape2=1) > > Please correct me if my understanding is wrong: > In the fitdistr fucntion we are providing the initial values of the > Beta distribution parameters as shape1=1 and shape2=1. This function > will try to fit the data and give us the new parameters of Beta > distribution that approximately fits this data. > > > I have tried the function with other distribution like Normal, Gamma, > Weibull...its working fine.. > > Regards, > Som Shekhar > > > > > On May 4, 1:25 am, Steven Kennedy <stevenkennedy2...@gmail.com> wrote: >> library(MASS) >> fitdistr(x,"beta",list(shape1=1,shape2=1)) >> >> >> >> On Tue, May 3, 2011 at 9:44 PM, Shekhar <shekhar2...@gmail.com> wrote: >> >> > Hi, >> > I have some random data and i want to find out the parameters of Beta >> > distribution ( a and b) such that this data approximately fits into >> > this distribution. I have tried by plot the histograms and graph, but >> > it requires lot of tuning and i am unable to do that. can anyone tell >> > me how to do it programmitically in R? >> >> > Regards, >> > Som Shekhar >> >> > ______________________________________________ >> > r-h...@r-project.org mailing list >> >https://stat.ethz.ch/mailman/listinfo/r-help >> > PLEASE do read the posting guidehttp://www.R-project.org/posting-guide.html >> > and provide commented, minimal, self-contained, reproducible code. >> >> ______________________________________________ >> r-h...@r-project.org mailing listhttps://stat.ethz.ch/mailman/listinfo/r-help >> PLEASE do read the posting guidehttp://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. > ______________________________________________ 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.
