sas生存分析作图生存概率从0累积

andy162639

14864

收藏 2015-07-08

悬赏 100 个论坛币已解决

sas生存分析作图，生存概率是从1下降的，是否有参数可以设置，如果没有，怎么调整图让生存概率从0开始向上呢谢谢
proc lifetest data=a method=km PLOT=(S) censoredSYMBOL='*'
timelist= 0 12 24 36 48 60 72;
time p1004*p1003(0);
strata group;
run;

最佳答案

jingju11 查看完整内容

Sorry, my bad. I asked my senior coworker and she said we could do the plot by: plots=(survival(failure)) in PROC LIFETEST; JingJu

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

全部回复

jingju11

2015-7-8 15:23:23

andy162639 发表于 2015-7-9 14:16
有没有示例供参考呢？谢谢

Sorry, my bad. I asked my senior coworker and she said we could do the plot by:
plots=(survival(failure)) in PROC LIFETEST;

JingJu

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jingju11

2015-7-9 06:49:13

无法用该过程直接做CDF.但是cdf=1-surv.把数据取出,然后用SGPLOT作出即可.好象我的拨可以前谈过这个问题.
其实生存做图,很多人都不需要利用过程的输出.使用SGPLOT可以得到满足个人需要的图示.
京剧

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

andy162639

2015-7-9 14:16:49

jingju11 发表于 2015-7-9 06:49
无法用该过程直接做CDF.但是cdf=1-surv.把数据取出,然后用SGPLOT作出即可.好象我的拨可以前谈过这个问题.
...

有没有示例供参考呢？谢谢

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

yangz98

2015-7-10 13:49:42

小试了一下：
Modeling the cumulative incidence.

Data BMT; Set SASHelp.BMT;Run;
Data Cov;
  Length Group $13.;
  Group = "ALL"; Output;
  Group = "AML-Low Risk"; Output;
  Group = "AML-High Risk"; Output;
Run;
Ods Graphics On/ Width = 7In Height = 5In;
Proc PHReg Data = BMT Plots(Overlay) = CIF;
  Ods Select  CIFPlot;
  Class Group(Ref = First);
  Model T*Status(0) = Group / Eventcode = 1;
  Baseline Covariates = Cov Out = Res CIF = _ALL_ / Seed = 1000;
Run;
Proc Lifetest Data = BMT Plots = Survival(Outside Atrisk(Maxlen = 13 Outside(0.15) ) );
  Ods Select SurvivalPlot;
  Time T*Status(0);
  Strata Group / Order = Internal;
Run;
Ods Graphics Off;

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

yangz98

2015-7-11 10:40:53

jingju11 发表于 2015-7-9 06:49
无法用该过程直接做CDF.但是cdf=1-surv.把数据取出,然后用SGPLOT作出即可.好象我的拨可以前谈过这个问题.
...

我猜可以使用
PROC UNIVARIATE PLOTS ;
CDF Statement;
WHERE STATUS = 1;
Run;

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

点击查看更多内容…

jingju11

2015-7-11 21:48:36

yangz98 发表于 2015-7-11 10:40
我猜可以使用
PROC UNIVARIATE PLOTS ;
CDF Statement;

lifetest的要旨是处理生存分析。不能纯粹考虑事件，也要兼顾缺失。所以你的猜测是错误的。京剧

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

yangz98

2015-7-11 23:37:51

jingju11 发表于 2015-7-11 21:48
lifetest的要旨是处理生存分析。不能纯粹考虑事件，也要兼顾缺失。所以你的猜测是错误的。京剧

哈哈 CDF(Cumulative density function)貌似是针对生存时间的拟合过程吧, CIF (Cumulative incidence function)貌似和CDF不太一样吧

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

jingju11

2015-7-12 10:29:09

没有理解你的意思。他要求的1-survival function.我想一些人成其为CDF（according to wiki: the survival function is the complementary cumulative distribution function.）。然后你说CDF
可以用univariate作。我认为不是一回事。然后你又开始说CIF。
不过。我们只是纯粹的讨论问题。只是争论而已。
京剧

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

yangz98

2015-7-12 20:15:41

jingju11 发表于 2015-7-12 10:29
没有理解你的意思。他要求的1-survival function.我想一些人成其为CDF（according to wiki: the survival f ...

谢谢大神点拨撒刚刚做了功课自己理解错误了生存分析中的CDF和Univariate众的CDF还是不一样的累积分布函数（CDF, Cumulative distribution function）F(t) = 1 - S(t); S(t) = Pr(T > t), 这里的S(t)生存率是考虑了删失的。

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

yangz98

2015-7-12 20:27:19

jingju11 发表于 2015-7-12 10:29
没有理解你的意思。他要求的1-survival function.我想一些人成其为CDF（according to wiki: the survival f ...

http://www.ats.ucla.edu/stat/sas ... urvival/default.htm 没有理解透彻再好好研究研究吧不过还是非常感谢jingju大神

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

yangz98

2015-7-13 12:32:00

jingju11 发表于 2015-7-11 21:48
lifetest的要旨是处理生存分析。不能纯粹考虑事件，也要兼顾缺失。所以你的猜测是错误的。京剧

哈哈生存分析处理删失时，我们假设全是随机删失；而针对右删失类型的生存资料，univariate过程步的结果应该没错

扫码加我拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

littlecan2006

2015-8-3 17:31:07

yangz98 发表于 2015-7-10 13:49
小试了一下：
Modeling the cumulative incidence.

您好，请问这种类型的生存曲线是怎么作出来的？（具体的程序过程是什么？）。我买了冯国双的《医学实验分析与SAS实现》及《SAS统计分析实用宝典》多种书籍，学了好久都没学会，尤其是下面那几行NO. at risk的数据是怎么做出来的！

另外，下面这两种图我也不会做！如果您会做的话，能否教我一下呢（最好是具体的程序过程）？不胜感激！另外，有什么好的SAS书籍推荐么？