clear all
close all


 load 'C:\STAT\Survival\Survivaldat\limphoma.mat'
 tA = limphoma.survive(limphoma.stage == 1);
            cA = limphoma.censored(limphoma.stage == 1);
 tB = limphoma.survive(limphoma.stage == 2);
            cB = limphoma.censored(limphoma.stage == 2);
  
  ple(tA, cA); 
  hold on
  ple(tB, cB); 
  legend off
  str1 = 'Stage III'; str2 = 'Stage IV';
       text(200, 0.75, str1, 'fontsize',16);
       text(100, 0.20, str2, 'fontsize',16);    
   hold off
  
  %print -depsc 'C:\STAT\Survival\Survivaleps\limphomaLR.eps'
  
  [p chi2] = manteltwo(tA,  cA,  tB,  cB)

  
