2
$\begingroup$

i have just started out in MCMC and I am not sure if I fully understand the concepts of MCMC with respect to the above terms. Let me try to explain that in my own words and include some thoughts / questions that I still have on it.

My understanding for MCMC (in Bayesian inference) is that it help us identify the final distribution of the samples (also known as the posterior distribution). And we can then use it to determine the sample's mean, range etc.

Example: We wanted to find out the mean test score of students. We are given 2 samples (i.e. the test scores of 2 students only and let's say the score is 65 and 75). We conduct the MCMC.

So to implement MCMC, we need the following "ingredients" - Target distribution: This is actually the posterior distribution that we actually want the MCMC to sample from. In other words, we need to determine a distribution of the mean score.

Q1: If our goal is to determine the final distribution of the samples and yet we are required to specify the target distribution that the MCMC is working on, isnt that like forcing the solution to a distribution that we already have in mind?

Q2: What if we have no idea what this distribution is supposed to be, including the mean and standard deviation. What should we do to proceed from here?

  • Proposal distribution: I see this as a distribution to determine the step / noise to be added on the parameter when we sample from the target distribution. Let say we had a starting test score at 68. We use the proposal distribution to generate a random noise (e.g. 3) and this is added to the test score for testing. So the next state is 78.

Q3: Am i right in my interpretation here?

  • Likelihood function: As far I understand so far, this is used in the acceptance testing and is also where we incorporate our observed data in (ie. 2 students test score we have on hand).

Q4: Did i interpret this wrongly? So should the incorporation of our observed data be when we are determining the target distribution?

I think I may be already confused. I really hope you can help me out with this. Thank you in advance.

$\endgroup$

1 Answer 1

1
$\begingroup$

The question reflects on a misunderstanding of the respective roles of MCMC algorithms and Bayesian inference principles, and calls for some prior investment in the basics or core of what constitutes Bayesian inference. The answer below cannot address this missing background and I cannot but recommend investing some time in an introductory Bayesian textbook like Albert or Kruschke.

MCMC is a tool to simulate from a posterior distribution, not a method of (Bayesian) inference per se. This means that it requires two inputs:

  • the distribution of the data, written as a density $f(x;\theta)$, which also writes as the likelihood $\ell(\theta|x)$. This function has to be specified completely, possibly as a marginal across latent variables, i.e., as an integral over these variables. In which cases the latent variables are added to the parameter to define an augmented target, of which the true target is a marginal. (Pseudo-marginal Metropolis-Hastings techniques partake from the same principle.)
  • the prior distribution on the parameter $\theta$, written as a density $\pi(\theta)$. Again, this function has to be specified completely. (Or again be expressed as a marginal.)

Given those two elements, MCMC aims at simulating from the distribution with density $\pi(\theta|x)$ which is proportional to $\pi(\theta)\times \ell(\theta|x)$. It can operate without the normalising constant. And without these elements, MCMC makes no sense.

Concerning the proposal distribution, this is a technical component in some of the MCMC algorithms, namely a new value for the Markov chain is proposed from a certain distribution and accepted or rejected according to an acceptance probability, computed in such a way that the stationary distribution of the Markov chain is the posterior.

$\endgroup$
7
  • $\begingroup$ Thanks Xi'an. I think my misunderstanding in concepts comes more from MCMC than Bayesian. In Bayesian, if we can have a simple prior (e.g. Gaussian) and simple prior (e.g. Gaussian), it would be easy to calculate the posterior distribution. In that case, we wont need MCMC. However, sometimes, the likelihood and/or the prior distribution can look horrendous and calculating the posterior will be difficult. This is when MCMC becomes useful. Yet, from my reading on MCMC so far, we needed target distribution which is the est. posterior distribution. So i think I may miss some link in MCMC. $\endgroup$
    – LearnW
    Commented Mar 27, 2019 at 13:55
  • $\begingroup$ Therefore this lead to my first question: "If our goal is to determine the final distribution of the samples and yet we are required to specify the target distribution that the MCMC is working on, isnt that like forcing the solution to a distribution that we already have in mind?" $\endgroup$
    – LearnW
    Commented Mar 27, 2019 at 13:59
  • $\begingroup$ I do not understand what you mean by (a) estimated posterior distribution or by (b) final distribution of the sample. The likelihood function has to be numerically available or expressed as the marginal of a numerically available joint distribution. Or estimable by an unbiased estimator as in pseudo-marginal techniques. Otherwise MCMC does not apply and more approximate techniques like VB or ABC enter. $\endgroup$
    – Xi'an
    Commented Mar 27, 2019 at 16:08
  • $\begingroup$ Thanks Xi'an again for your help. Maybe I should clarify first if I am right that in implementing the MCMC we need a target distribution to be expressed? If yes, this target distribution that put in the MCMC system needs to be as close as the posterior distribution that we are trying to calculate right? If yes again, then i am thinking the reason we need MCMC is because we cannot really calculate the posterior but yet we are required to input in a target distribution. This is the part that is confusing me...thank you again so much for the help. $\endgroup$
    – LearnW
    Commented Mar 27, 2019 at 16:20
  • 1
    $\begingroup$ Great, i think i got it now. I misunderstand that we need a target distribution to be expressed in MCMC (i think the article that i have read have expressed a target distribution to illustrate/show how close is the result of MCMC with the target distribution or the final posterior distribution). Thank you again Xi'an :) $\endgroup$
    – LearnW
    Commented Mar 27, 2019 at 16:48

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.