function val=tlike2(par,data) % log posterior of t-error problem % on darwin's data mu=par(:,1); lsigma=par(:,2); sigma=exp(lsigma); df=4; n=length(data); val=0*mu; c1=(df+1)/2; c2=df*sigma.^2; for i=1:n val=val-c1*log(1+(data(i)-mu).^2./c2)-lsigma; end