Sampling from normalized and un-normalized posterior

  • Context: Graduate 
  • Thread starter Thread starter michaelmas
  • Start date Start date
  • Tags Tags
    Sampling
Join the discussion
Registration is free. Ask a follow-up in this thread, or start your own.
5 replies · 2K views
michaelmas
Messages
3
Reaction score
0
Help me understand something.

I get that the posterior ##p(\theta|y) \propto p(y|\theta)p(\theta)## should be normalized by ##\frac{1}{p(y)}## for the probability to sum to 1, but what about the mean and variance?

Am I not right understanding that sampling from the un-normalized posterior gives the same mean and variance as sampling from the normalized posterior?

Can I prove it mathematically?

Can't find it and can't figure it out.
 
Physics news on Phys.org
michaelmas said:
Am I not right understanding that sampling from the un-normalized posterior gives the same mean and variance as sampling from the normalized posterior?

.

How do you define "sampling" from a distribution that doesn't integrate to 1.0 ?
 
I don't know.

What exactly is the point with MCMC and not normalizing?

Isn't that what the Bayesians are doing?
 
You should explain where in the Markov Chain Monte Carlo method that you think sampling is done from a non-normalized distribution. As far as I know, that technique is never defined or used.
 
Ok, let me rephrase the question.

if ##p(\theta|y)## is the distribution of interest, then what good is ##p(y|\theta)p(\theta)## if the mean and variance aren't the same?
 
michaelmas said:
if ##p(\theta|y)## is the distribution of interest, then what good is ##p(y|\theta)p(\theta)## if the mean and variance aren't the same?

You are speaking as if you've read that one distribution is of some "good" in answering questions about the other, but until you say exactly what the "good" is, it isn't possible to understand what you are asking.

In general, if [itex]W[/itex] is a random variable and [itex]X = k\ W[/itex] for some constant [itex]k[/itex] then you can figure out the mean and variance of [itex]X[/itex] if you know the mean and variance of [itex]W[/itex]. If we treat [itex]y[/itex] a fixed event then [itex]k = \frac{1}{p(y)}[/itex].