mcmc模拟代码

mcmc模拟有详细的描述，对mcmc需要了解的人有很大帮助

% BETALPR - Beta Distribution - Log Probability Ratio
% Copyright （c） 1998 Harvard University. Full copyright in the file Copyright
% 
%   [ lpr ] = betalpr（p1p2alphabeta）
%
% returns the log of the p（p1） / p（p2） when both
% are distributed Beta（alphabeta）.
%
% See also: METROP *LPR

function [ lpr ] = betalpr（p1p2alphabeta）
lpr = （alpha-1） * log（p1/p2） + （beta-1） * log（（1-p1）/（1-p2）） ;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  1999-02-22 15:27  mcmc\
     文件        1479  1999-02-22 15:11  mcmc\About_MCMC.html
     文件         399  1999-02-22 15:11  mcmc\betalpr.m
     文件        3085  1999-02-22 15:11  mcmc\Contents.m
     文件         487  1999-02-22 15:11  mcmc\Copyright
     文件         360  1999-02-22 15:11  mcmc\gamlpr.m
     文件        1479  1999-02-22 15:11  mcmc\index.html
     文件         851  1999-02-22 15:11  mcmc\invwishirnd.m
     文件         842  1999-02-22 15:11  mcmc\invwishlpr.m
     文件        1251  1999-02-22 15:11  mcmc\invwishpdf.m
     文件         842  1999-02-22 15:11  mcmc\invwishrnd.m
     文件         689  1999-02-22 15:11  mcmc\ltindex.m
     文件         554  1999-02-22 15:11  mcmc\ltvec.m
     文件         346  1999-02-22 15:11  mcmc\mcmcacf.m
     文件        1426  1999-02-22 15:11  mcmc\mcmcdemo.m
     文件         731  1999-02-22 15:11  mcmc\mcmcgr.m
     文件         786  1999-02-22 15:11  mcmc\mcmclt.m
     文件        2722  1999-02-22 15:11  mcmc\mcmcsumm.m
     文件         777  1999-02-22 15:11  mcmc\mcmctrace.m
     文件         860  1999-02-22 15:11  mcmc\metrop.m
     文件         641  1999-02-22 15:11  mcmc\mvnormlpr.m
     文件         898  1999-02-22 15:11  mcmc\mvnormrnd.m
     文件         504  1999-02-22 15:11  mcmc\randrand.m
     文件         601  1999-02-22 15:11  mcmc\veclt.m
     文件          29  1999-02-22 15:11  mcmc\Version
     文件         771  1999-02-22 15:11  mcmc\wishirnd.m
     文件        1454  1999-02-22 15:11  mcmc\wishrnd.m