%Script to save a set of gibbs samples for the lab
figure(3)
for(datatype = 1:7)
   data=getdata(datatype,20);
   tmp=gibbs_sampler_reml_like(data', 1000,50000);
   saved_gibbs_reml{datatype}=tmp;
   subplot(3,3,datatype)
   hist(log(tmp.s2p));
   title(datatype);
end
save saved_gibbs_reml saved_gibbs_reml;