Tell me more ×
Cross Validated is a question and answer site for statisticians, data analysts, data miners and data visualization experts. It's 100% free, no registration required.

I have a data set that contains two types of points. The first type of points come from an N(0,1) distribution. The second type of points come from an N(m,v) distribution for some real m and some positive, real v. The objective is to classify each point as type 1 or type 2, and to identify m & v. We have no apriori information about m & v. Any ideas?

share|improve this question
You can find more ideas and code in "Separating Two Populations in a Sample," stats.stackexchange.com/questions/899/…. (Is there perhaps some better way of cross-referencing threads than I've done here?) – whuber Aug 18 '10 at 21:49

1 Answer

up vote 1 down vote accepted

You could use a mixture model to separate out the components. The data generating process can be represented as follows:

Let:

$z_i$: be the type (1 or 2) for the $i^{th}$ observation,

$y_i$ be the $i^{th}$ observation.

Then you have:

$f(y_i|z_i=1) \sim N(0,1)$

$f(y_i|z_i=2) \sim N(m,v)$

$P(z_i=1) = \pi$ and

$P(z_i=2) = 1-\pi$.

The likelihood function is given by:

$L(m,v,\pi|-) = \sum_{y_i} \pi f(y_i|z_i=1) + (1-\pi) f(y_i|z_i=2)$

You can then use either EM or MCMC to estimate for the model parameters.

share|improve this answer
Have no idea why latex looks so bad. Can someone please help and fix this? – user28 Aug 9 '10 at 18:05

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.