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

From: <mmiller_at_mail.usp.edu>
Date: Mon, 18 Feb 2008 16:03:10 +0000

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


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)

Where B=background, H=hits, N=noise, and D=my observed distribution

I have seen analyses like this in the past, but I can't remember what it is called. If somebody out there can point me towards an R function, or even the cannonical name for this kind of model, I think I can write the necessary code.

Thanks in advance,

R-help_at_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. Received on Mon 18 Feb 2008 - 16:53:55 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Tue 19 Feb 2008 - 02:30:15 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive