[R] newbie (me) needs to model distribution as two overlapping gaussians

David Winsemius dwinsemius at comcast.net
Mon Feb 18 21:36:56 CET 2008


<mmiller at mail.usp.edu> wrote in
news:BAY105-W146A57D6AFC49F8003C73E8F210 at phx.gbl: 

> 
> Recently, I have been working with some data that look like two
> overlapping gaussian distributions.  I would like to either 
> 
> 1) determine the mean and SD for each of the two distributions
> 
> OR
> 
> 2) get some (bayesian ?) statistic that estimates how likely an
> observation is to belong to the left-hand or right-hand distribution 
> 
> In case I'm using the wrong language, my data looks something like
> this: 
> 
> B <- rnorm(500,40,10)
> H <- rnorm(500,80,5 )
> N <- runif(200,0,99)
> D <- c(B,H,N)

First hit on <http://www.rseek.org/> with: functions mixture gaussians

<http://research.stowers-institute.org/efg/R/Statistics/MixturesOfDistributions/index.htm>

Task Views:
<http://cran.us.r-project.org/web/views/Cluster.html>
<http://cran.r-project.org/web/views/Multivariate.html>
<http://cran.r-project.org/web/views/Bayesian.html>

nor1mix: Normal (1-d) Mixture Models (S3 Classes and Methods)
<http://lib.stat.cmu.edu/R/CRAN/>  

In the nor1mix documentation, Mächler recommends the mclust package 
for estimation.

-- 
David Winsemius



More information about the R-help mailing list