R. python新手日记R语言从零开始

R语言日常笔记(5)一些小问题的集合

2019-07-27  本文已影响0人  柳叶刀与小鼠标

(1)生存分析的KM曲线绘制问题
在绘制之前,我们会

group <- ifelse(gene > = median(gene), 'high', 'low')

一般常用中位数将样本分为高低表达组,这样便于绘制,但是假如说某个基因表达量为0的样本数目超过了半数,这样的结果就是所有该基因的所有表达量被修改成‘high’,这样会导致,生存曲线绘制错误。


image.png

error in ggsurvplot_df(d, fun = fun, color = color, palette = palette, :
The length of legend.labs should be 1

修改的代码是

group <- ifelse(gene > median(gene), 'high', 'low') 取消等号

==============================================

对表达量矩阵进行反log2+1


a <- data.frame(a = c(1, 1, 5, 6, 8, 9),
                b = c(2, 7, 5, 2, 6, 9),
                c = c(1, 4, 9, 0, 8, 4))


a_new <-   data.frame(matrix(nrow=nrow(a),ncol=ncol(a)))


for (i in 1:nrow(a)) {
  for (j in 1:ncol(a)) {
   
    num <- a[i,j]
    a_new[i,j] <- 2^num - 1
    
  }
  
}

rownames(a_new) <- rownames(a)
colnames(a_new) <- colnames(a)

a
a_new
上一篇 下一篇

猜你喜欢

热点阅读