%WinBUGS wavelets
clear all
close all


  disp('Wavelet Shrinkage in WinBUGS')
  lw = 2; 
  set(0, 'DefaultAxesFontSize', 16);
  fs = 15;
  msize = 5;
  

% 
% This mat file has three columns 
% 1. input -> wavelet coefs d
 ds =[ -50.0000  -10.0000   -5.0000   -4.0000   -3.0000   -2.0000   -1.0000 ...
   -0.5000   -0.1000         0    0.1000    0.5000    1.0000    2.0000...
    3.0000    4.0000    5.0000   10.0000   50.0000 ];
 % 2. WinBUGS output  E(p)=0.1
 dst1 = [ -49.7300   -6.2970   -1.6700   -1.2000   -0.8320   -0.5193   -0.2434...
   -0.1175   -0.0171   -0.0055    0.0201    0.1162    0.2417    0.5246...
    0.8288    1.1970    1.6600    6.2790   49.7200];
 % 3. WinBUGS output  E(p)=0.05  
 dst2 =[ -49.6500   -4.6460   -0.9458   -0.6825   -0.4808   -0.2960   -0.1466...
   -0.0669   -0.0159   -0.0029    0.0112    0.0739    0.1406    0.3016 ...
    0.4796    0.6803    0.9403    4.6770   49.6700];
figure(1)
plot(ds, ds, ':','linewidth',lw)
hold on
plot(ds, dst2, 'k-','linewidth',lw)
plot(ds, dst2, 'ko','markersize',msize)
print -deps 'C:\NPBook\NPBayes\wavbugs.eps'