Closed mikelove closed 10 years ago
for example, limma and edgeR give a column "logFC", DESeq gives "log2FoldChange"
## Calculate fold coverage vs group 1 if(length(regionGroupMean) > 1){ foldChange <- vector("list", length(regionGroupMean) - 1) names(foldChange) <- names(regionGroupMean)[-1] for(group in names(foldChange)) { foldChange[[group]] <- log2(regionGroupMean[[group]] / regionGroupMean[[1]]) } }
Good idea. We would need to change downstream uses of foldChange with log2FoldChange. Could you issue a pull request with these changes made? That would be huge.
Fixed in 78b70658b35fb30dc12b79237d68ef1710b5b712
for example, limma and edgeR give a column "logFC", DESeq gives "log2FoldChange"