function fkmt = kme_all_x(dt) % this function performs Kaplan-Meier estimation with given parameter % and gives estimated F() for all data points [kmest sortdat] = KMcdfSM(dt(:,1), dt(:,2), 0); 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]; temp_val = []; %calculate each CDF F() value for all data points for i=1:length(data) if sum(data(i) >= sortdat) > 0 temp_val = [temp_val kmest(sum(data(i) >= sortdat), 1)]; else % when there is no observation, CDF F() value is simply 0 temp_val = [temp_val 0]; end end fkmt = temp_val;