diff --git a/detectRUNS/DESCRIPTION b/detectRUNS/DESCRIPTION
index a4e7595..e30e09c 100644
--- a/detectRUNS/DESCRIPTION
+++ b/detectRUNS/DESCRIPTION
@@ -28,7 +28,7 @@ Imports:
     Rcpp,
     gridExtra,
     data.table
-RoxygenNote: 7.1.1
+RoxygenNote: 7.2.3
 Suggests:
     testthat,
     knitr,
diff --git a/detectRUNS/R/plots.R b/detectRUNS/R/plots.R
index 0ac2db9..66b648c 100644
--- a/detectRUNS/R/plots.R
+++ b/detectRUNS/R/plots.R
@@ -783,7 +783,7 @@ plot_InbreedingChr<- function(runs, mapFile , groupSplit=TRUE, style=c("ChrBarPl
     g1 <- g1 +  scale_x_discrete(labels=list_chr)
     g1 <- g1 +  xlab("Inbreeding by Chromosome") + ylab("Froh")
     if(!is.null(plotTitle)) { g1 <- g1 +  ggtitle(mainTitle1) + theme(plot.title = element_text(hjust = 0.5)) }
-    if (groupSplit) { g1 <- g1 + facet_grid(group ~. ) + guides(fill=FALSE) }     # if you want split or not!
+    if (groupSplit) { g1 <- g1 + facet_grid(group ~. ) + guides(fill="none") }     # if you want split or not!
     if (savePlots){ ggsave(filename = fileNameOutput1 , plot = g1, device = "pdf") } else { print(g1) }
   }
 
@@ -795,7 +795,7 @@ plot_InbreedingChr<- function(runs, mapFile , groupSplit=TRUE, style=c("ChrBarPl
     g2 <- g2 + scale_x_discrete(labels=list_chr)
     g2 <- g2 + xlab("Inbreeding by Chromosome") + ylab("Froh")
     if(!is.null(plotTitle)) { g2 <- g2 +  ggtitle(mainTitle2) + theme(plot.title = element_text(hjust = 0.5)) }
-    if (groupSplit) { g2 <- g2 + facet_grid(group ~. ) + guides(fill=FALSE) }    # if you want split or not!
+    if (groupSplit) { g2 <- g2 + facet_grid(group ~. ) + guides(fill="none") }    # if you want split or not!
     if (savePlots){ ggsave(filename = fileNameOutput2 , plot = g2, device = "pdf") } else { print(g2) }
   }
 
@@ -949,7 +949,7 @@ plot_DistributionRuns <- function(runs, mapFile , groupSplit=TRUE, style=c("Mean
     g1 <- g1 + geom_bar(stat="identity", position=position_dodge())
     g1 <- g1 + xlab("Class Length Category") + ylab("Mean (Mb)") + scale_x_discrete(limits=unique(long_DF$chrom))
     g1 <- g1 +  ggtitle(mainTitle1) + theme(plot.title = element_text(hjust = 0.5))
-    if (groupSplit) { g1 <- g1 + facet_grid(group ~. ) + guides(fill=FALSE) }    # if you want split or not!
+    if (groupSplit) { g1 <- g1 + facet_grid(group ~. ) + guides(fill="none") }    # if you want split or not!
     if (savePlots){ ggsave(filename = fileNameOutput1 , plot = g1, device = "pdf") } else { print(g1) }
   }
 
@@ -962,7 +962,7 @@ plot_DistributionRuns <- function(runs, mapFile , groupSplit=TRUE, style=c("Mean
     g2 <- g2 + geom_bar(stat="identity", position=position_dodge())
     g2 <- g2 + xlab("Chromosome") + ylab("Mean (Mb)") + scale_x_discrete(limits=unique(long_DF$chrom))
     g2 <- g2 +  ggtitle(mainTitle2) + theme(plot.title = element_text(hjust = 0.5))
-    if (groupSplit) { g2 <- g2 + facet_grid(group ~. ) + guides(fill=FALSE) }    # if you want split or not!
+    if (groupSplit) { g2 <- g2 + facet_grid(group ~. ) + guides(fill="none") }    # if you want split or not!
     if (savePlots){ ggsave(filename = fileNameOutput2 , plot = g2, device = "pdf") } else { print(g2) }
   }
 
@@ -974,7 +974,7 @@ plot_DistributionRuns <- function(runs, mapFile , groupSplit=TRUE, style=c("Mean
     g3 <- g3 + geom_bar(stat="identity", position=position_dodge()) + scale_x_discrete(limits=unique(long_DF$CLASS))
     g3 <- g3 + xlab("Class Length Category") + ylab("Frequency")
     g3 <- g3 +  ggtitle(mainTitle3) + theme(plot.title = element_text(hjust = 0.5))
-    if (groupSplit) { g3 <- g3 + facet_grid(group ~. ) + guides(fill=FALSE) }    # if you want split or not!
+    if (groupSplit) { g3 <- g3 + facet_grid(group ~. ) + guides(fill="none") }    # if you want split or not!
     if (savePlots){ ggsave(filename = fileNameOutput3 , plot = g3, device = "pdf") } else { print(g3) }
   }
 
@@ -987,7 +987,7 @@ plot_DistributionRuns <- function(runs, mapFile , groupSplit=TRUE, style=c("Mean
     g4 <- g4 + geom_bar(stat="identity", position=position_dodge()) + scale_x_discrete(limits=unique(long_DF$chrom))
     g4 <- g4 + xlab("Chromosome") + ylab("Frequency")
     g4 <- g4 +  ggtitle(mainTitle4) + theme(plot.title = element_text(hjust = 0.5))
-    if (groupSplit) { g4 <- g4 + facet_grid(group ~. ) + guides(fill=FALSE) }    # if you want split or not!
+    if (groupSplit) { g4 <- g4 + facet_grid(group ~. ) + guides(fill="none") }    # if you want split or not!
     if (savePlots){ ggsave(filename = fileNameOutput4 , plot = g4, device = "pdf") } else { print(g4) }
   }
 }
diff --git a/detectRUNS/man/tableRuns.Rd b/detectRUNS/man/tableRuns.Rd
index 83c2da8..1e2bf7d 100644
--- a/detectRUNS/man/tableRuns.Rd
+++ b/detectRUNS/man/tableRuns.Rd
@@ -27,8 +27,9 @@ proportion of individuals carrying that run (e.g. 70\%)}
 }
 \value{
 A dataframe with the most common runs detected in the sampled individuals
-(the group/population, start and end position of the run, chromosome and number of SNP
-included in the run are reported in the output dataframe)
+(the group/population, start and end position of the run, chromosome, number of SNP
+included in the run and average percentage of SNPs in run
+are reported in the output dataframe)
 }
 \description{
 This function takes in input either the run results or the output from
diff --git a/detectRUNS/src/RcppExports.cpp b/detectRUNS/src/RcppExports.cpp
index c2b5579..e4721e6 100644
--- a/detectRUNS/src/RcppExports.cpp
+++ b/detectRUNS/src/RcppExports.cpp
@@ -5,6 +5,11 @@
 
 using namespace Rcpp;
 
+#ifdef RCPP_USE_GLOBAL_ROSTREAM
+Rcpp::Rostream<true>&  Rcpp::Rcout = Rcpp::Rcpp_cout_get();
+Rcpp::Rostream<false>& Rcpp::Rcerr = Rcpp::Rcpp_cerr_get();
+#endif
+
 // fast_factor
 SEXP fast_factor(SEXP x);
 RcppExport SEXP _detectRUNS_fast_factor(SEXP xSEXP) {