Random number from a unknown distribution

  • Thread starter aihaike
  • Start date
  • #1
55
0

Main Question or Discussion Point

Dear all,

I apologize if it is the wrong place, I don't know where I had to post this question since I'm not a mathematician.

Well, suppose you have a set of numbers which can be describe by a unknown distribution.
I just like to know whether we can use those numbers to generate a set of random number so that they will have the same probability distribution.

More precisely, I have a set of atomic velocities which do not obey to the Maxwell distribution. I'd like to generate a new (random) set of velocity. I thought about multiplying each initial value by a uniform random number between -1 and 1 but it does not seem to work.

Thank you,

Éric.
 

Answers and Replies

  • #2
chiro
Science Advisor
4,790
132
Hey aihaike.

In terms of generating a random number process that generates something from an arbitrary distribution, there are techniques for doing this, but you need to be able to describe your distribution at least by some sort of histogram. It doesn't need to be an analytic one, but it needs to have at least some kind of bin-representation (i.e. histogram). You could fit a curve to the histogram but this is not necessary.

Once you have enough samples to build a good histogram, you can then use something like the Metropolis-Hastings algorithm to generate a random number from your histogram distribution.

http://en.wikipedia.org/wiki/Metropolis–Hastings_algorithm

You will need to implement some code to get the maximum and I recommend you use the uniform distribution where the constant * the uniform PDF "boxes in" your histogram distribution.

Also you will need to have your histogram be defined on a finite-interval. If it's defined on the entire real line or equivalent in multiple dimensions, you need to create a finite-interval that captures the majority of the distribution that's good enough. For example getting five standard deviations from the mean of a normal is going to be more than enough for most purposes.

Once you use MH and use the right test distribution taking into account your histogram and making sure the assumptions are good for using MH (there are other techniques you can look into), then this is your random number generator for your histogram distribution.

Just make sure you get a lot of samples to construct the histogram.
 
  • #3
55
0
Hey Chiro,

Thank you for your quick and useful reply!
I looks like exactly what I was looking for.
I'll try it out.
Thanks again.

Éric.
 

Related Threads on Random number from a unknown distribution

  • Last Post
Replies
15
Views
4K
  • Last Post
Replies
3
Views
2K
Replies
4
Views
616
  • Last Post
Replies
13
Views
3K
  • Last Post
Replies
3
Views
1K
Replies
1
Views
10K
  • Last Post
Replies
7
Views
6K
Replies
12
Views
5K
Replies
8
Views
771
Replies
8
Views
737
Top