function ws = MultiBayes2(wc,L)
%
%%%%%%%%%%%%%%%It's my own function!!!!!
%
% MultiBayes2 -- Apply Shrinkage to Wavelet Coefficients 2-d
%  Usage 
%    ws = MultiBayes(wc,L)
%  Inputs 
%    wc    2-d Wavelet Transform of noisy sequence with N(0,1) noise
%    L     low-frequency cutoff for Wavelet Transform
%  Outputs 
%    ws    result of applying BayesThresh to each dyadic block
%
	[n,J] = dyadlength(wc);
	ws = wc;
    for j=(J-1):-1:L
        [t1,t2]=dyad2HH(j);
        ws(t1,t2)=BayesThresh2(wc(t1,t2));
        [t1,t2]=dyad2HL(j);
        ws(t1,t2)=BayesThresh2(wc(t1,t2));
        [t1,t2]=dyad2LH(j);
        ws(t1,t2)=BayesThresh2(wc(t1,t2));
    end