Allgene) <- Allgene$Geneidallgeneexpression <- Allgene[,3:9]geneLength <- Allgene$LengthtpmMatrix <- convertCounts
结果是一致的方法二:R包 DGEobj.utilsload("gfe.Rdata")#提取上面处理好的基因长度列effLen = gfe$lengthlibrary(DGEobj.utils)CC_res <- convertCounts 求高手指点~综合上述分析,暂时还是选择常规方法和DGEobj.utils R包中的convertCounts函数吧~参考资料:1、https://hbctraining.github.io/Training-modules /planning_successful_rnaseq/lessons/sample_level_QC.html2、https://rdrr.io/cran/DGEobj.utils/man/convertCounts.html3