(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') 取消等号