Hello,
Is there any test for bimodality in R that
x <- c(rnorm(1000,0,1),rnorm(1000,3,1))
hist(x,nclass=100)
Thank you in advance for any help.
Regards,
Samor
[[alternative HTML version deleted]]
Samor, A somewhat indirect answer: you could fit a mixture distribution to your data and test how many components are needed to best describe your data. hth, Ingmar On Wed, Feb 24, 2010 at 1:22 PM, Samor Gandhi <samorgandhi@yahoo.com> wrote:> Hello, > > Is there any test for bimodality in R that > > x <- c(rnorm(1000,0,1),rnorm(1000,3,1)) > hist(x,nclass=100) > > Thank you in advance for any help. > > Regards, > Samor > > > > [[alternative HTML version deleted]] > > > ______________________________________________ > 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. > >[[alternative HTML version deleted]]
You could use: library(TeachingDemos) x <- c(rnorm(1000,0,1),rnorm(1000,3,1)) TeachingDemos:::vis.test( x, TeachingDemos:::vt.normhist ) Then click on the plot that looks most different. Sorry about needing TeachingDemos::: for now, the functions were accidentally left out of the NAMESPACE file, this has been fixed and will not be needed for version 2.6 on. -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.snow at imail.org 801.408.8111> -----Original Message----- > From: r-help-bounces at r-project.org [mailto:r-help-bounces at r- > project.org] On Behalf Of Samor Gandhi > Sent: Wednesday, February 24, 2010 5:23 AM > To: r-help at r-project.org > Subject: [R] Bimodal distribution > > Hello, > > Is there any test? for bimodality in R that > > x <- c(rnorm(1000,0,1),rnorm(1000,3,1)) > hist(x,nclass=100) > > Thank you in advance for any help. > > Regards, > Samor > > > > [[alternative HTML version deleted]]