diff --git a/workflow/main.nf b/workflow/main.nf
index 658835264a92f08f1dfecc433e7ea7fc944345fa..b12867bf600de39425aac9bb2c4e81877e55da13 100644
--- a/workflow/main.nf
+++ b/workflow/main.nf
@@ -100,10 +100,9 @@ process diffbind {
    file "*"
    """
    module load R/3.2.1-intel
-   cp design.txt design.shiny.txt
-   cp geneset.gmt geneset.shiny.gmt
-   Rscript  $baseDir/scripts/dea.R
-   perl $baseDir/scripts/concat_edgeR.pl *.edgeR.txt
+   #cp design.txt design.shiny.txt
+   #cp geneset.gmt geneset.shiny.gmt
+   Rscript  $baseDir/scripts/runDiffBind.R
  """
 }
 
diff --git a/workflow/scripts/runDiffBind.R b/workflow/scripts/runDiffBind.R
index e95d45421cead97070c3c93cf60d29725c9763fc..2eb624ba86c7b252c8ed7108eb2cac1e37338929 100644
--- a/workflow/scripts/runDiffBind.R
+++ b/workflow/scripts/runDiffBind.R
@@ -16,7 +16,7 @@ dba.plotPCA(data, DBA_TISSUE, label=DBA_CONDITION)
 dev.off()
 
 #Save peak reads count
-normcount <- dba.peakset(dta, bRetrieve=T)
+normcount <- dba.peakset(data, bRetrieve=T)
 write.table(as.data.frame(normcount),"diffbind.normcount.txt",sep="\t",quote=F,row.names=F)
 
 #Retriving the differentially bound sites