diff --git a/workflow/rules/annotations.smk b/workflow/rules/annotations.smk
index 6b9f040..966e0b9 100644
--- a/workflow/rules/annotations.smk
+++ b/workflow/rules/annotations.smk
@@ -103,7 +103,7 @@ rule combine_homer:
         rscript=join(SCRIPTSDIR,"_combine_macs2_homer.R")
     shell:
         """
-        Rscript {params.rscript} {input.peaks_file} {input.annotation} {output.combined}
+        Rscript {params.rscript} {input.peaks_file} {input.annotation} {output.combined_tsv} {output.combined_xlsx}
         """
 
 rule rose:
diff --git a/workflow/scripts/_combine_macs2_homer.R b/workflow/scripts/_combine_macs2_homer.R
index c36804e..88b8d1b 100644
--- a/workflow/scripts/_combine_macs2_homer.R
+++ b/workflow/scripts/_combine_macs2_homer.R
@@ -1,5 +1,6 @@
 # Add peak statistics to HOMER annotations file
 
+library(openxlsx)
 library(readr)
 library(dplyr)
 library(tidyr)