Binomical vs poisson distribution in Simulations

Click For Summary
SUMMARY

This discussion focuses on the application of Binomial and Poisson distributions in simulating DNA mutations using C++. The user, Hagen, seeks to determine whether a binomially distributed random variable can replace a Poisson distributed random variable for simulating mutations in DNA sequences. The conversation highlights the statistical relationship between the two distributions, specifically that the Poisson distribution can be approximated by the binomial distribution under certain conditions, particularly when the number of trials is large and the probability of success is small.

PREREQUISITES
  • Understanding of Binomial Distribution
  • Understanding of Poisson Distribution
  • Basic knowledge of C++ programming
  • Familiarity with statistical concepts such as mean and variance
NEXT STEPS
  • Research the mathematical relationship between Binomial and Poisson distributions
  • Learn how to implement random number generation in C++ using libraries like random
  • Study the implications of using Poisson approximation for large sample sizes
  • Explore statistical simulation techniques in C++ for modeling biological processes
USEFUL FOR

This discussion is beneficial for statisticians, data scientists, and software developers interested in computational biology, particularly those working on simulations involving genetic mutations and statistical modeling.

hagen
Hey, I want to write a Computer Simulation in C++, which simulates the development of a DNA sequence with a probability to mutate x in one "generation". I do have a variable number (=n) of copies of this DNA. Now one might think, to simulate the mutation by simply:
sum(n*Poisson distributed random variable(x) )​
to get the number of mutated DNA copies. But this would be too slow.
So my question is, could I also just create a
binomically distributed random variable and multiply it by n * x​
to get the number of mutated DNA's? Or is this statistically incorrect?
If not, how might I set the Params for the Bin. dis.? Can I take 1 as a mean and multiply the result x or has the mean to be x? And how do I set / transform the variance of the distribution in a ratio to number of copies.
As you might probably have guessed, I'm a beginner in statistics, but i would be really grateful for any help. Thanks in advance,
Hagen
 
Physics news on Phys.org
Hm, if I remember correct, there is a link between the binomial and Poisson distribution.. Poisson's probability function is given with [tex]f(x)=\frac{\lambda^x}{x!}e^{-\lambda}[/tex]. Now, I think you can put [tex]\lambda=mp[/tex], where the number of repetitions of a Bernoulli scheme experiment [tex]m \rightarrow \infty[/tex] and it's probability [tex]p \rightarrow 0[/tex].
 

Similar threads

  • · Replies 12 ·
Replies
12
Views
5K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 13 ·
Replies
13
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
5K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 15 ·
Replies
15
Views
3K