Giters

lcolladotor / derfinder

Annotation-agnostic differential expression analysis of RNA-seq data via expressed regions-level or single base-level approaches

Home Page:http://lcolladotor.github.io/derfinder

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

rename "foldChange" => "log2FoldChange" for the column in the output regions

mikelove opened this issue · comments

commented

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]])
                        }
                }
commented

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.