How can the Metropolis-Hastings algorithm help simulate a Normal Distribution?

  • Thread starter Thread starter Zaare
  • Start date Start date
  • Tags Tags
    Algorithm
Zaare
Messages
54
Reaction score
0
I'm having trouble understanding how to find an expression for \pi(x) and \pi(y) in the relation:
<br /> \alpha \left( {x,y} \right) = \min \left( {1,\frac{{\pi \left( y \right)q\left( {y,x} \right)}}{{\pi \left( x \right)q\left( {x,y} \right)}}} \right)<br />
For example, If I want to simulate Normal Distribution (Expectation value m and standard deviation s), how can I find expressions for \pi(x) and \pi(y)? Or are they equal: \pi(x)=\pi(y)?
 
Physics news on Phys.org
I've never used this algorithm, but the values of x and y are calculated/generated according to the algorithm and are different. One of the values, depending on how you are defining x and y, should come from the proposal distribution.

And pi(x) = the normal distribution probability of x, for your particular example.
 
I see. I was hopeing I could find a "short cut" which would simplify the expression for pi(x), but I suppose I can't.
Thanks for the help.
 
Namaste & G'day Postulate: A strongly-knit team wins on average over a less knit one Fundamentals: - Two teams face off with 4 players each - A polo team consists of players that each have assigned to them a measure of their ability (called a "Handicap" - 10 is highest, -2 lowest) I attempted to measure close-knitness of a team in terms of standard deviation (SD) of handicaps of the players. Failure: It turns out that, more often than, a team with a higher SD wins. In my language, that...
Hi all, I've been a roulette player for more than 10 years (although I took time off here and there) and it's only now that I'm trying to understand the physics of the game. Basically my strategy in roulette is to divide the wheel roughly into two halves (let's call them A and B). My theory is that in roulette there will invariably be variance. In other words, if A comes up 5 times in a row, B will be due to come up soon. However I have been proven wrong many times, and I have seen some...
Back
Top