% Demo of wavelet-based function estimation
clear all
close all
        %--- set plotting parameters -------
        lw = 2;  
        set(0, 'DefaultAxesFontSize', 15);
        fs = 15;
        msize = 10;

data = [36.3, 41.7, 43.9, 49.9, 50.1, 50.8, 51.9, 52.1, 52.3, 52.3, 52.4, ...
52.6, 52.7, 53.1, 53.6, 53.6, 53.9, 53.9, 54.1, 54.6, 54.8, 54.8, 55.1, ...
55.4, 55.9,56.0, 56.1, 56.5, 56.9, 57.1, 57.1, 57.3, 57.7, 57.8, 58.1, ...
58.9, 59.0, 59.1, 59.6, 60.4, 60.7, 26.8, 29.6, 33.4, 35.0, 40.0, 41.9, 42.5];
censor=[ones(1,41), zeros(1,7)];
ci = ciboot([data' censor'], 'kme_all_x', 5, .95, 1000);
   figure;	
   plot(data', ci(:,2)', '.');
   hold on;
   plot(data', ci(:,1)', '+');
   plot(data', ci(:,3)', '*');	
print -deps 'C:\NPBook\Bootstrap\sky.eps'