disp('Running MATLAB script file SiZerDist1Fig11.m') ;
%
%    For Making Figure 11
%    in the Advanced Distribution Theory for SiZer paper
%
%
%    Copied from: SiZerSimul2power.m via SiZerDist1Fig9.m
%
%    For investigating power (using standard test bed examples)
%    of improved SiCon methods
%
%    Copied from sizerSMtest.m
%
%    FOR DEVELOPMENT AND TESTING OF MATLAB FUNCTION sizerSM,
%    Default version of SiCon,


%vitest = [2] ;
                %  14 ER Wavelet Examples (regression)


%iout = 1 ;    %  1 - 4 panel SiCon plot, Regular, Hsing-R, H-Bonf, H-Global

addpath 'D:\Research\SiZer\simul';
%path

fighand = figure(1) ;
clf ;

matfilestr = '\Research\GeneralData\incomes' ;
vrange = [0;3] ;

load(matfilestr) ;


flag = data <= 3 ;
data = data(flag) ;

disp('        Working on Hsing Global SiCon') ;
paramstruct = struct('titlestr',['Figure 11'],...
                     'ibdryadj',0, ...
                     'simflag',3) ;
paramstruct = setfield(paramstruct,'vxgp',[vrange(1);vrange(2);401]) ;

  %  Note 'ibdryadj',0 is DIFFERENT from the 2 used in sims
  %  simflag is:
  %       3  -  Use global adjustment (sum independent blocks calculated as in 2)

sc2SMsim(data,paramstruct) ;





pstr = 'SiZerDist1Fig11' ;

  orient portrait ;
  set(gcf,'PaperSize',[6.0, 2.0]) ; 
  set(gcf,'PaperPosition',[0.25, 0.25, 5.5, 1.5]) ; 
print('-dpsc', [pstr '.ps']) ;

