表达量太低的话取log 可能出现负值的;
医学癌症TCGA-生存分析课程中基因表达数据处理时通过cmp对表达数据进行标准化处理
我分析的是mRNA的数据
运行下面代码,得到的数据如下图
norExpr <- DGEList(counts=datExpr)
norExpr <- calcNormFactors(norExpr)
运行下面的代码后,得到数据如下图
norExpr_cpm <- cpm(norExpr,prior.count=2, log=TRUE)
cpm标准化后,数据中有很多是负值,这是不是有问题呀?因为课程里提到cpm计算后数值比较大,所以选择进行log进行数据标准化。
而且后续进行单因素cox回归分析时,运行到
# 将单因素的结果汇总,形成一个dataframe
Univar <- ldply(Univar, data.frame)
得不到结果。
请各位老师和同学指点迷津,十分感谢!