w*******e 发帖数: 666 | 1 刚开始学习R,主要是想用survivalROC package。想请教下面两个问题,非常感谢!
第一个问题:
我做10年和1年的time dependt ROC得到的AUC 一样。我用的是survial ROC guide上的
code。做十年的时候只是把cutoff改成了3650还有year改成了10。不知道是不是还有什
么需要改的,为什么结果完全没有区别?
第二个问题是:我使用survivalROC 是因为我用传统的ROC 做出来我正在研究的那个
marker对outcome的AUC只有0.5几。做survivalROC 做出来也只有0.5几。有朋友提了下
面建议,说在其他factor加进去的时候看AUC是多少,然后再把我正在做的marker加进
去。看AUC能增加多少。我想请问survivalROC有这个功能吗?我在网上搜索到的资料就
只有这个package的guide。所以,不太清楚该怎么做,我是否找到了正确的package! |
Y******Y 发帖数: 8753 | 2 不清楚survivalROC有没有这个功能。不过就算没,你直接做个likelihood ratio test
也可以的
!
【在 w*******e 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 刚开始学习R,主要是想用survivalROC package。想请教下面两个问题,非常感谢! : 第一个问题: : 我做10年和1年的time dependt ROC得到的AUC 一样。我用的是survial ROC guide上的 : code。做十年的时候只是把cutoff改成了3650还有year改成了10。不知道是不是还有什 : 么需要改的,为什么结果完全没有区别? : 第二个问题是:我使用survivalROC 是因为我用传统的ROC 做出来我正在研究的那个 : marker对outcome的AUC只有0.5几。做survivalROC 做出来也只有0.5几。有朋友提了下 : 面建议,说在其他factor加进去的时候看AUC是多少,然后再把我正在做的marker加进 : 去。看AUC能增加多少。我想请问survivalROC有这个功能吗?我在网上搜索到的资料就 : 只有这个package的guide。所以,不太清楚该怎么做,我是否找到了正确的package!
|
w*******e 发帖数: 666 | 3 请问用R里的什么程序做呢?
算C statistic也可以。但是我也不知道是用哪个package。
test
【在 Y******Y 的大作中提到】![](/moin_static193/solenoid/img/up.png) : 不清楚survivalROC有没有这个功能。不过就算没,你直接做个likelihood ratio test : 也可以的 : : !
|
Y******Y 发帖数: 8753 | 4 lrtest(mod1,mod2)就行, from epicalc package
不过最简单的用手呀,把两个model的likelihood拿出来,减一下,算好df,查
chisqaure distribution的probability就行。 |
w*******e 发帖数: 666 | 5 非常感谢!
【在 Y******Y 的大作中提到】![](/moin_static193/solenoid/img/up.png) : lrtest(mod1,mod2)就行, from epicalc package : 不过最简单的用手呀,把两个model的likelihood拿出来,减一下,算好df,查 : chisqaure distribution的probability就行。
|