logpdf or pdf in mhsample
2 次查看（过去 30 天）
I am using 'mhsample' function for Metropolis-Hastings sampling from a posterior distribution. I am not sure if I should use pdf or logpdf as input arguments. For example, if the prior is a Gaussian, can I give the 'logpdf' argument instead of 'pdf'. What are the situations in which I should use logpdf? And how does MatLab processes 'logpdf' differently than 'pdf', does it make any difference in my final solution.
the cyclist 2017-12-1
The functions that are input arguments to mhsample are probability density functions. If you are actually using the logarithm of the probability density functions as inputs, then specify that with 'logpdf'.
Since you say your prior is gaussian, I'm guessing that that is simply the density (not the log), so you should use 'pdf'.