[R] newbie (me) needs to model distribution as two overlapping gaussians
Greg Snow
Greg.Snow at imail.org
Tue Feb 19 18:31:50 CET 2008
The mixreg package may help. Just fit an intercept only model to your
data (may need to have x be a column of 1's and tell it not to include
an intercept) and see if it finds your mixtures.
Hope this helps,
--
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 tcornulier
> Sent: Monday, February 18, 2008 12:48 PM
> To: mmiller at mail.usp.edu
> Cc: r-help at r-project.org
> Subject: Re: [R] newbie (me) needs to model distribution as
> two overlapping gaussians
>
> Is a finite mixture of 2 gaussians the name you are looking for?
> This specific model will not deal with your N component however.
>
> you can find some functions here:
> http://cran.r-project.org/web/views/Cluster.html
>
> Thomas
>
> mmiller at mail.usp.edu wrote:
> > 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)
> >
> > 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,
> > Mark
> >
> > ______________________________________________
> > 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.
> >
>
> --
>
>
> **************************
> Thomas Cornulier
>
> Mammal Research Institute
> Polish Academy of Sciences
> ul. Waszkiewicza 1c
> 17-230 Bialowieza
> http://www.zbs.bialowieza.pl/?en
>
> tel (0048) 85 682 77 88
>
> ______________________________________________
> 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.
>
More information about the R-help
mailing list