Biology Stack Exchange is a question and answer site for biology researchers, academics, and students. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I wonder why gene expression data are very frequently modeled by multivariate normal distributions. What is the reason for those strong assumptions that the genes follow multivariate gaussian? Are there any reasons specific for genetics other than the reasons for general gaussian assumptions (the ease of calculation, etc.)?

share|improve this question

Usually if something is not expected to behave according to some scheme then the measured values for such a parameter are assumed to be normal. It is not just with gene expression but with all types of measurements like dimensions of an object, luminosity of an electric bulb, range of a bullet etc. In any measurement, the random error is modeled using normal distribution. I don't have a very intuitive explanation for why random errors follow normal distribution but mathematically it comes from the central limit theorem.

Now each gene is a variable and measurement of each gene suffers some random error; so a multivariate normal distribution is used.

When we dismiss a null hypothesis in a t-test or z-test what we are actually doing is dismissing our parsimonious notion that a sample is drawn from a given normal distribution. This means two things:

  1. The sample belongs to some other normal distribution (different $\mu$ and $\sigma$)
  2. The sample follows some other distribution

But a t-test will never be able to point out the exact reason. All it tells you is that the sample is not from some given normal distribution.

share|improve this answer

Your Answer


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.