全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
12831 6
2015-06-13
使用 survival package

>km<-survdiff(formula = Surv(overallsurvival, status) ~ antt[,54],data=antt)
>km
Call:survdiff(formula = Surv(overallsurvival, status) ~ antt[, 54],data = antt)
N Observed Expected (O-E)^2/E (O-E)^2/V
antt[, 54]=C 77 36 41.43 0.711 5.46
antt[, 54]=mut 21 12 6.57 4.484 5.46
Chisq= 5.5 on 1 degrees of freedom, p= 0.0195

在最后面有个p=0.0195 如何把这个P值提取出来?
> summary(km)
Length Class Mode
n 2 table numeric
obs 2 -none- numeric
exp 2 -none- numeric
var 4 -none- numeric
chisq 1 -none- numeric
call 3 -none- call  
> str(km)
List of 6
$ n : 'table' int [1:2(1d)] 77 21
..- attr(*, "dimnames")=List of 1
.. ..$ groups: chr [1:2] "antt[, 54]=C" "antt[, 54]=mut"
$ obs : num [1:2] 36 12
$ exp : num [1:2] 41.43 6.57
$ var : num [1:2, 1:2] 5.4 -5.4 -5.4 5.4
$ chisq: num 5.46$ call : language survdiff(formula = Surv(overallsurvival, status) ~ antt[, 54], data = antt)- attr(*, "class")= chr "survdiff"

km 为包含6个成分的list变量
使用km$n,km$obs,km$exp,km$var,km$chisq, km$call 均得到不到p 值 。

求大神赐教!
二维码

扫码加我 拉你入群

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

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

全部回复
2015-6-13 17:05:26
学习了!
二维码

扫码加我 拉你入群

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

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

2015-6-13 17:11:33
提示: 作者被禁止或删除 内容自动屏蔽
二维码

扫码加我 拉你入群

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

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

2015-6-13 18:42:09
ryoeng 发表于 2015-6-13 17:11
https://stat.ethz.ch/pipermail/r-help/2007-April/130676.html
果然厉害!!佩服佩服!
二维码

扫码加我 拉你入群

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

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

2015-6-13 18:59:08
查了下ryong提供的链接,

实现办法:
p.value <- 1 - pchisq(km$chisq, length(km$n) - 1)

就是利用卡方值跟自由度算概率值啊。。。[lol]
二维码

扫码加我 拉你入群

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

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

2016-1-21 20:22:10
MARK  
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群