om pmtk3.googlecode.com


%PMTKauthor Emtiyaz Khan
%PMTKneedsStatsToolbox iwishrnd

data = [59  77
43  39
34  46
32  26
42  38
38  43
55  68
67  86
64  77
45  60
49  50
72  59
34  38
70  48
34  55
50  58
41  54
52  60
60  75
34  47
28  48
35  33];

params.mean0 = [50 50]';
params.covMat0 = [625 312.5; 312.5 625];
params.nu0 = 4;
params.S0 = params.covMat0;

setSeed(0);
nSamples = 1000;
% data should be D*N
samples = gibbsSampler(data', params, nSamples);
burnin = 250;
samples.mu = samples.mu(:,burnin:end);
