User:Timothee Flutre/Notebook/Postdoc/2011/12/14

From OpenWetWare

Jump to: navigation, search
Project name Main project page
Previous entry      Next entry

Learn about mixture models and the EM algorithm

  • Motivation and examples: it's frequent to collect heterogeneous data, ie. for which we suspect that they come from several clusters. For instance, we measure the height of individuals without recording their gender, we measure the levels of expression of a gene in several individuals without recording which ones are healthy and which ones are sick, etc.
  • Data: we have N observations, noted X = (x1,x2,...,xN). For the moment, we suppose that each observation xi is univariate, ie. each corresponds to only one number.
  • Hypotheses and aim: let's assume that the data are heterogeneous and that they can be partitioned into K clusters (see examples above). This means that a subset of the observations come from cluster k = 1, another subset come from cluster k = 2, and so on.
  • Model: technically, we say that the observations were generated by a family of density functions. The density of all the observations is thus a mixture of densities, one per cluster. In our case, we will assume that each cluster k corresponds to a Normal distribution of mean μk and standard deviation σk. Moreover, as we don't know for sure from which cluster a given observation comes from, we define the mixture probability wk to be the probability that any given observation comes from cluster k. As a result, we have the following list of parameters: θ = (w1,...,wK1,...μK1,...,σK. Finally, for a given observation xi, we can write the model f(x_i/\theta) = \sum_{k=1}^{K} w_k g(x_i/\mu_k,\sigma_k) , wth g being the Normal distribution g(x_i/\mu_k,\sigma_k) = \frac{1}{\sqrt{2\pi} \sigma_k} \exp^{-\frac{1}{2}(\frac{x_i - \mu_k}{\sigma_k})^2}

Personal tools