%%
  close all
  clear all
  disp('Cox Oakes leukemia Remission Data Set')
  lw = 2.5; 
  set(0, 'DefaultAxesFontSize', 16);
  fs = 15;
  msize = 10;
 
 
 times =  [6 6 6 6 7 9 10 10  11 13 16 17 19 20 22 23 25 32 32 34 35];
 censor =[1 0 0 0 0 1  1   0    1   0    0    1   1   1   0    0  1    1    1   1   1];
[f,t] = ecdf(times,'censoring',censor);
[t 1-f]
figure()
[f,x,flo,fup] = ecdf(times,'censoring', censor);
figure(2)
stairs(x,1-f,'LineWidth',2)
table =[x 1-f];
hold on
stairs(x,1-flo,'r:','LineWidth',2)
stairs(x,1-fup,'r:','LineWidth',2)
legend('Empirical','LCB','UCB',...
       'Location','SW')
hold off

  
 %%
 
 times =  [6 6 6 6 7 9 10 10  11 13 16 17 19 20 22 23 25 32 32 34 35];
 censor =[0 0 0 1 0 1  1   0    1   0    0    1   1   1   0    0  1    1    1   1   1];
%plorkm(times, censor)
 table = ple(times, censor)