前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >R语言之生信(9)R语言多个生存分析曲线比较

R语言之生信(9)R语言多个生存分析曲线比较

作者头像
用户1359560
发布2019-03-22 11:12:59
5.7K0
发布2019-03-22 11:12:59
举报
文章被收录于专栏:生信小驿站

====================================== 这篇博客的目的主要是计算当需要计算多个不同组之间的成对比较,并计算P值。

主要使用的函数是:pairwise_survdiff(formula,data,p.adjust.method =“BH”,na.action,rho = 0)

  • 参数formula:类似其他生存模型的公式表达式,形式为Surv(time,status)〜variable。
  • 参数data:一个数据框,用于做生存分析的数据。
  • 参数p.adjust.method:p值矫正方法(参见p.adjust)。 允许的参数包含(“holm”,“hochberg”,“hommel”,“bonferroni”,“BH”,“BY”,“fdr”,“none”)。 如果不想对p值矫正(不推荐),请使用p.adjust.method =“none”。
  • 参数na.action:缺失数据过滤功能。
  • 参数RHO参数:用于控制测试类型。 允许值包括0(Log-Rank检验)和1(peto和peto检验)。 值
  • 函数计算返回对象是包含p值的列表。
代码语言:javascript
复制
> library(survival)
> library(survminer)
> data(myeloma)
> 
> # Pairwise survdiff
> res <- pairwise_survdiff(Surv(time, event) ~ molecular_group,
+                          data = myeloma)
> res

查看输出结果

代码语言:javascript
复制
    Pairwise comparisons using Log-Rank test 

data:  myeloma and molecular_group 

                 Cyclin D-1 Cyclin D-2 Hyperdiploid Low bone disease MAF   MMSET
Cyclin D-2       0.723      -          -            -                -     -    
Hyperdiploid     0.943      0.723      -            -                -     -    
Low bone disease 0.723      0.988      0.644        -                -     -    
MAF              0.644      0.447      0.523        0.485            -     -    
MMSET            0.328      0.103      0.103        0.103            0.723 -    
Proliferation    0.103      0.038      0.038        0.062            0.485 0.527

P value adjustment method: BH 

通过这个矩阵我们可以看到molecular_group 水平上两两之间的生存分析对比的P值。 下一步我们可以将数字转化成符号,代表我们不同的P值水平。

代码语言:javascript
复制
> symnum(res$p.value, cutpoints = c(0, 0.0001, 0.001, 0.01, 0.05, 0.1, 1),
+        symbols = c("****", "***", "**", "*", "+", " "),
+        abbr.colnames = FALSE, na = "")
                 Cyclin D-1 Cyclin D-2 Hyperdiploid Low bone disease MAF MMSET
Cyclin D-2                                                                    
Hyperdiploid                                                                  
Low bone disease                                                              
MAF                                                                           
MMSET                                                                         
Proliferation               *          *            +                         
attr(,"legend")
[1] 0 ‘****’ 1e-04 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘+’ 0.1 ‘ ’ 1 \t    ## NA: ‘’
本文参与 腾讯云自媒体同步曝光计划,分享自作者个人站点/博客。
原始发表:2019.03.22 ,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 作者个人站点/博客 前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档