Skip to content

Commit

Permalink
Update DiffBind_v2_EdgeR_block.Rmd
Browse files Browse the repository at this point in the history
Removing up.bed and down.bed
  • Loading branch information
tovahmarkowitz authored Sep 12, 2024
1 parent 1398faa commit e0db6d0
Showing 1 changed file with 1 addition and 13 deletions.
14 changes: 1 addition & 13 deletions bin/DiffBind_v2_EdgeR_block.Rmd
Original file line number Diff line number Diff line change
Expand Up @@ -35,8 +35,6 @@ csvfile <- params$csvfile
contrasts <- params$contrasts
peakcaller <- params$peakcaller
peak_counts <- params$counts
down_file <- params$down_file
up_file <- params$up_file
list_file <- params$list_file
```

Expand Down Expand Up @@ -178,16 +176,6 @@ dba.plotHeatmap(DBAnalysisEdgeR, contrast=1, method=DBA_EDGER_BLOCK, correlation
## Top 500 or less differentially bound peaks

```{r EdgeRReport}
UpPeaks <- DBReportEdgeR[which(DBReportEdgeR$Fold > 0)]
print("up ")
print(dim(UpPeaks))
rtracklayer::export(UpPeaks, up_file)
DownPeaks <- DBReportEdgeR[which(DBReportEdgeR$Fold < 0)]
print("down ")
print(dim(DownPeaks))
rtracklayer::export(DownPeaks, down_file)
D2i <- length(DBReportEdgeR)
i <- as.integer(min(c(500, as.integer(max(c(D2i, 1))))))
DT::datatable(data.frame(DBReportEdgeR)[1:i,], rownames=F)
Expand All @@ -209,4 +197,4 @@ write.table(report2, list_file, quote=F, sep="\t", row.names=F)
sessionInfo()
```

<div class="tocify-extend-page" data-unique="tocify-extend-page" style="height: 0;"></div>
<div class="tocify-extend-page" data-unique="tocify-extend-page" style="height: 0;"></div>

0 comments on commit e0db6d0

Please sign in to comment.