Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Home
- /
- Analytics
- /
- Stat Procs
- /
- Proc lIfetest : Different results from product-limit-estimate table an...

Options

- RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

🔒 This topic is **solved** and **locked**.
Need further help from the community? Please
sign in and ask a **new** question.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Posted 05-12-2020 03:40 AM
(810 views)

Hi,

My code -

```
ods output productlimitestimates = sds ;
proc lifetest data=graph2_B_G timelist= 0 100 200 300 400 500 600 700 800 900 1000 plots=(survival( atrisk=0 to 1000 by 100 ) )maxtime = 1000;
time time*status(0);
strata drg_class_detail;
run;
```

The numbers that i get from the product-limit-estimate table (number surviving) are not the same as the numbers i get on the graph with the atrisk number. Can anyone please help explain why that is the case?

What is the difference in the calculation?

even the numbers for the time-0 are different

Eg;

Atrisk | Product-lim-estimate | |||

time | group1 | group2 | group1 | group2 |

0 | 2050 | 2050 | 2043 | 2048 |

100 | 1454 | 1481 | 1451 | 1480 |

200 | 1186 | 1217 | 1182 | 1213 |

300 | 894 | 848 | 892 | 845 |

1 ACCEPTED SOLUTION

Accepted Solutions

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Hi @riyaaora275,

Your survival times are discrete. The numbers of subjects at risk in the graph (denoted by *Y*_{i} in section Breslow, Fleming-Harrington, and Kaplan-Meier Methods of the PROC LIFETEST documentation) are counted "*just prior to*" the respective time points, whereas the columns "Number Left" and "Number Failed" in the table "Product-Limit Survival Estimates" (i.e., variables Left and Failed in your dataset sds) take the events (failures or censoring) *at* the respective time points already into account (event by event) -- not at time 0, though, when using the TIMELIST= option. So the small differences are the numbers of subjects whose event times *equal* 0, 100, 200, etc., respectively.

1 REPLY 1

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Hi @riyaaora275,

Your survival times are discrete. The numbers of subjects at risk in the graph (denoted by *Y*_{i} in section Breslow, Fleming-Harrington, and Kaplan-Meier Methods of the PROC LIFETEST documentation) are counted "*just prior to*" the respective time points, whereas the columns "Number Left" and "Number Failed" in the table "Product-Limit Survival Estimates" (i.e., variables Left and Failed in your dataset sds) take the events (failures or censoring) *at* the respective time points already into account (event by event) -- not at time 0, though, when using the TIMELIST= option. So the small differences are the numbers of subjects whose event times *equal* 0, 100, 200, etc., respectively.

Are you ready for the spotlight? We're accepting content ideas for **SAS Innovate 2025** to be held May 6-9 in Orlando, FL. The call is **open **until September 16. Read more here about **why** you should contribute and **what is in it** for you!

What is ANOVA?

ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.

Find more tutorials on the SAS Users YouTube channel.