Services & Resources / Wolfram Forums / MathGroup Archive
-----

MathGroup Archive 2010

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: I look for a bi-modal distribution

  • To: mathgroup at smc.vnet.net
  • Subject: [mg114462] Re: I look for a bi-modal distribution
  • From: Andy <andyr at wolfram.com>
  • Date: Sun, 5 Dec 2010 21:50:15 -0500 (EST)

On 12/4/2010 5:11 AM, Alexei Boulbitch wrote:
> Dear community,
>
> I have an experimental data list having the form of lst={{x1, y1},...}.
> I make the analysis of statistics of the data. The distribution is
> non-Gaussian, as can be checked by various test criteria, and
> visualized by a histogram. I suspect however, that here there may be two
> superimposed Gaussian distributions in which both the means and the
> standard deviations are different.
>
> Among  the Mathematica built-in distributions in I did not succeed to
> find such a bi-modal distribution.  Do you know anything about such a
> distribution?
>
> Regards, Alexei
>
You can use MixtureDistribution to accomplish this. I'm assuming from 
the structure of your data that this is a bivariate distribution.

dist = MixtureDistribution[ { w1, w2 }, 
{MultinormalDistribution[{mu11,mu12},Sigma1], 
MultinormalDistribution[{mu21,mu22}, Sigma2]

Here the wi's are weights which determine how much each component 
contributes to the overall distribution. Note that the Sigmai's are 
covariance matrices. Also note that MixtureDistribution can take any 
distributions you want, not just gaussian distributions.

You can estimate the parameters in the component distributions and the 
weights with EstimatedDistribution.

Andy Ross
Wolfram Research


  • Prev by Date: Re: How to use "Apply" to do differentiation ?
  • Next by Date: Re: Changing Graph style
  • Previous by thread: I look for a bi-modal distribution
  • Next by thread: Re: I look for a bi-modal distribution