注意:我建议使用我的新GWMCMC采样器，也可以从文件交换下载://www.tatmou.com/matlabcentral/fileexchange/49820-the-mcmc-hammer--gwmcmc
马尔可夫链蒙特卡罗抽样后验分布

大都市采样器
[mmc, logP] =密度(initialm、loglikelihood logmodelprior, stepfunction, mccount,跳过)
---------
Initialm:随机游走的起始点
loglikelihood:函数句柄到似然函数:logL(m)
logprior:对数模型的函数句柄先验概率:logPapriori(m)
Stepfunction:没有输入的函数句柄，返回一个随机的
在随机游走中。(注意stepfunction也可以是a
描述正态分布大小的矩阵
步骤。)
mccount:马尔可夫链应该有多长?
跳过:通过只存储第N步来精简链[默认=10]



示例用法:拟合数据的正态分布
-------------------------------------------
data = randn (100 1) * 2 + 3;
logmodelprior = @ (m) 0;%使用平坦优先级。
loglike = @ (m)和(日志(normpdf(数据、m(1)、(2))));
minit = [0 1];
m =密度(minit loglike logmodelprior,(。2。5),10000);
米(:1:10 0)= [];%作物漂移
plotmatrix (m);


——Aslak Grinsted 2010