楼主: andy162639
13184 21

sas生存分析作图生存概率从0累积 [推广有奖]

  • 3关注
  • 6粉丝

教授

53%

还不是VIP/贵宾

-

威望
0
论坛币
5062 个
通用积分
229.8629
学术水平
10 点
热心指数
11 点
信用等级
0 点
经验
8851 点
帖子
849
精华
0
在线时间
1220 小时
注册时间
2010-5-5
最后登录
2024-4-14

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
关键词:生存分析 LIFETEST Censored censor Method method
有人说,统计归根结底就是个P
沙发
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
已有 3 人评分学术水平 热心指数 信用等级 收起 理由
wodewin + 1 + 1 + 1 精彩帖子
andy162639 + 3 + 3 + 2 精彩帖子
Tigflanker + 3 + 3 观点有启发

总评分: 学术水平 + 7  热心指数 + 7  信用等级 + 3   查看全部评分

使用道具

藤椅
jingju11 发表于 2015-7-9 06:49:13 |只看作者 |坛友微信交流群
无法用该过程直接做CDF.但是cdf=1-surv.把数据取出,然后用SGPLOT作出即可.好象我的拨可以前谈过这个问题.
其实生存做图,很多人都不需要利用过程的输出.使用SGPLOT可以得到满足个人需要的图示.
京剧
已有 1 人评分论坛币 学术水平 热心指数 信用等级 收起 理由
andy162639 + 1 + 1 + 1 + 1 精彩帖子

总评分: 论坛币 + 1  学术水平 + 1  热心指数 + 1  信用等级 + 1   查看全部评分

使用道具

板凳
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 |只看作者 |坛友微信交流群


SurvivalPlot8.png CIFPlot8.png
小试了一下:
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;

已有 1 人评分学术水平 热心指数 信用等级 收起 理由
andy162639 + 2 + 3 + 2 精彩帖子

总评分: 学术水平 + 2  热心指数 + 3  信用等级 + 2   查看全部评分

使用道具

地板
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;

使用道具

7
jingju11 发表于 2015-7-11 21:48:36 |只看作者 |坛友微信交流群
yangz98 发表于 2015-7-11 10:40
我猜可以使用
PROC UNIVARIATE PLOTS ;
  CDF Statement;
lifetest的要旨是处理生存分析。不能纯粹考虑事件,也要兼顾缺失。所以你的猜测是错误的。京剧

使用道具

8
yangz98 发表于 2015-7-11 23:37:51 |只看作者 |坛友微信交流群
jingju11 发表于 2015-7-11 21:48
lifetest的要旨是处理生存分析。不能纯粹考虑事件,也要兼顾缺失。所以你的猜测是错误的。京剧
哈哈 CDF(Cumulative density function)貌似是针对生存时间的拟合过程吧, CIF (Cumulative incidence function)貌似和CDF不太一样吧

使用道具

9
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。
不过。我们只是纯粹的讨论问题。只是争论而已。
京剧

使用道具

10
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)生存率是考虑了删失的。
已有 1 人评分论坛币 学术水平 热心指数 信用等级 收起 理由
andy162639 + 1 + 1 + 3 + 2 精彩帖子

总评分: 论坛币 + 1  学术水平 + 1  热心指数 + 3  信用等级 + 2   查看全部评分

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-19 16:41