diff --git a/workflow/main.nf b/workflow/main.nf index 40a9990263e1245c54c72c7c225ae1cc0bae803a..577903a78c3fda3a20a47a1ad104e27397e6f10b 100644 --- a/workflow/main.nf +++ b/workflow/main.nf @@ -417,10 +417,10 @@ process diffPeaks { output: - file "design_diffpeak_annotatePeaks.csv" into diffPeaksDesignAnnotatePeaks - file "design_diffpeak_annotatePeaks.csv" into diffPeaksDesignMeme - file "*_diffbind.bed" into diffpeaks_meme - file "*_diffbind.bed" into diffpeaks_chipseeker + file "design_diffpeak_annotatePeaks.tsv" into diffPeaksDesignAnnotatePeaks + file "design_diffpeak_annotatePeaks.tsv" into diffPeaksDesignMeme + file "*_diffbind.bed" into diffPeaks + file "*_diffbind.csv" into diffPeaksCounts file '*.pdf' into diffPeaksStats file 'normcount_peaksets.txt' into normCountPeaks diff --git a/workflow/scripts/dba.R b/workflow/scripts/dba.R index 3081db20f73b0e22313c9c2926e72faf3c35a3e1..42986d23afa1b5cc0c699ad27ec68468be8ea9df 100644 --- a/workflow/scripts/dba.R +++ b/workflow/scripts/dba.R @@ -42,7 +42,7 @@ for (i in c(1:length(data$contrasts))) { contrast_bed_name = paste(data$contrasts[[i]]$name1,"vs", data$contrasts[[i]]$name2,"diffbind.bed",sep="_") contrast_name = paste(data$contrasts[[i]]$name1,"vs", - data$contrasts[[i]]$name2,"diffbind.xls",sep="_") + data$contrasts[[i]]$name2,"diffbind.csv",sep="_") new_SampleID = c(new_SampleID,paste(data$contrasts[[i]]$name1,"vs",data$contrasts[[i]]$name2,sep="_")) new_Peaks = c(new_Peaks, contrast_bed_name) report <- dba.report(data, contrast=i, th=1, bCount=TRUE) @@ -55,4 +55,4 @@ for (i in c(1:length(data$contrasts))) { # Write new design file newdesign = data.frame(Condition=new_SampleID, Peaks=new_Peaks) -write.csv(newdesign,"design_diffpeak_annotatePeaks.csv",row.names=F, quote=F) +write.csv(newdesign,"design_diffpeak_annotatePeaks.csv",row.names=F, quote=F, sep='\t')