DESeq2 normalized
2020-07-03 本文已影响0人
一路向前_莫问前程_前程似锦
conds <- factor( rep("samples", times=ncol(countsTable)) )
cds <- newCountDataSet( countsTable[setdiff(rownames(countsTable), tailrows),], conds )
cds <- estimateSizeFactors( cds )
countsNorm <- counts( cds, normalized=TRUE )
countsNorm <- round(countsNorm, digits=1)
write.table(countsNorm, file=file.path(res.path, paste(tumorname, ".normalized.count.all.features.txt", sep="")), row.names=T, col.names=NA, sep="\t", quote=F)
