diff --git a/DESCRIPTION b/DESCRIPTION
index a3f1322..cf47618 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -1,7 +1,7 @@
 Package: tbea
-Title: Tools for Pre- and Post-processing in Bayesian Evolutionary Analyses
-Version: 1.4.1
-Authors@R: c(person("Gustavo A.", "Ballen", email = "gaballench@gmail.com", role = c("aut", "cre")), 
+Title: Tools for Pre- and Post-Processing in Bayesian Evolutionary Analyses
+Version: 1.4.2
+Authors@R: c(person("Gustavo A.", "Ballen", email = "gustavo.a.ballen@gmail.com", role = c("aut", "cre")), 
     person("Sandra", "Reinales", email = "spreinalesl@gmail.com", role = c("aut")))
 Description: Package for bayesian inference in phylogenetics and evolution.
     It provides functions for prior specification in divergence time
diff --git a/R/xintercept.R b/R/xintercept.R
index b80d46a..a3bc7dd 100644
--- a/R/xintercept.R
+++ b/R/xintercept.R
@@ -39,16 +39,15 @@
 #' ages <- ages[complete.cases(ages)] # remove NAs
 #' ages <- ages[which(ages < 10)] # remove outliers
 #' 
+#' \dontrun{
 #' # Draper-Smith, OLS
 #' draperSmithNormalX0 <- xintercept(x = ages, method = "Draper-Smith", alpha = 0.05, robust = FALSE)
 #' # Draper-Smith, Robust fit
-#' draperSmithRobustX0 <-xintercept(x = ages, method = "Draper-Smith", alpha = 0.05, robust = TRUE)
+#' draperSmithRobustX0 <- xintercept(x = ages, method = "Draper-Smith", alpha = 0.05, robust = TRUE)
 #' # Bootstrap, OLS
 #' bootstrapNormalX0 <- xintercept(x = ages, method = "Bootstrap", p = c(0.025, 0.975), robust = FALSE)
 #' # Bootstrap, Robust fit
 #' bootstrapRobustX0 <- xintercept(x = ages, method = "Bootstrap", p = c(0.025, 0.975), robust = TRUE)
-#'
-#' \donttest{
 #' # plot the estimations
 #' hist(ages, probability = TRUE, col = rgb(red = 0, green = 0, blue = 1, alpha = 0.3),
 #'      xlim = c(0, 10), main = "CDF-based on confidence intervals", xlab = "Age (Ma)")
diff --git a/man/xintercept.Rd b/man/xintercept.Rd
index 9fcba21..359eacb 100644
--- a/man/xintercept.Rd
+++ b/man/xintercept.Rd
@@ -49,16 +49,15 @@ ages <- andes$ages
 ages <- ages[complete.cases(ages)] # remove NAs
 ages <- ages[which(ages < 10)] # remove outliers
 
+\dontrun{
 # Draper-Smith, OLS
 draperSmithNormalX0 <- xintercept(x = ages, method = "Draper-Smith", alpha = 0.05, robust = FALSE)
 # Draper-Smith, Robust fit
-draperSmithRobustX0 <-xintercept(x = ages, method = "Draper-Smith", alpha = 0.05, robust = TRUE)
+draperSmithRobustX0 <- xintercept(x = ages, method = "Draper-Smith", alpha = 0.05, robust = TRUE)
 # Bootstrap, OLS
 bootstrapNormalX0 <- xintercept(x = ages, method = "Bootstrap", p = c(0.025, 0.975), robust = FALSE)
 # Bootstrap, Robust fit
 bootstrapRobustX0 <- xintercept(x = ages, method = "Bootstrap", p = c(0.025, 0.975), robust = TRUE)
-
-\donttest{
 # plot the estimations
 hist(ages, probability = TRUE, col = rgb(red = 0, green = 0, blue = 1, alpha = 0.3),
      xlim = c(0, 10), main = "CDF-based on confidence intervals", xlab = "Age (Ma)")