[R] How to fit a random data into Beta distribution?

[Shekhar]

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... at gmail.com> wrote:
> library(MASS)
> fitdistr(x,"beta",list(shape1=1,shape2=1))
>
>
>
> On Tue, May 3, 2011 at 9:44 PM, Shekhar <shekhar2... at 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... at 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... at 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.