diff --git a/R/evaluation_identity.R b/R/evaluation_identity.R
index 50538280046b642af3ef1492ed5bf6d1cd9b0559..76d9efac6fa51f625a90bb4de27bf05da4055ba9 100644
--- a/R/evaluation_identity.R
+++ b/R/evaluation_identity.R
@@ -61,6 +61,7 @@ get_rsquare_2plot <- function(data_rsquare){
 #'
 #' @param data_identity A data frame containing comparison results with "actual" and "estimate" columns.
 #' @param palette_color dict-like palette default: palette_color = c(DESeq2 = "#500472", HTRfit ="#79cbb8")
+#' @param palette_shape Optional parameter that sets the point shape for plots.Default : c(DESeq2 = 17, HTRfit = 19).
 #' @param ... additional parameters to pass geom_point aes 
 #' @return A ggplot2 identity plot and R2 metric associated
 #'
diff --git a/R/mock_rnaseq.R b/R/mock_rnaseq.R
index 8f60f98461b81909516bbc6fdb4098843fbae19c..a2227cdc1c1e2bf1c92e245f15a752a4ec5e6b97 100644
--- a/R/mock_rnaseq.R
+++ b/R/mock_rnaseq.R
@@ -96,8 +96,7 @@ mock_rnaseq <- function(list_var, n_genes, min_replicates, max_replicates, seque
   
   if (!is.null(sequencing_depth)) {
     message("Scaling count table according to sequencing depth.")
-    message("Scaling count table according to sequencing depth.")
-
+    message("Scaling counts by sequencing depth may exhibit some randomness due to certain parameter combinations, resulting in erratic behavior. This can be minimized by simulating more genes. We advise verifying the simulated sequencing depth to avoid drawing incorrect conclusions.")
     mu_ij_dtf_rep <- scaleCountsTable(mu_ij_matx_rep, sequencing_depth)
     mu_ij_matx_rep <- as.matrix(mu_ij_dtf_rep)
   }
diff --git a/R/prepare_data2fit.R b/R/prepare_data2fit.R
index e2cdd246261e2254e7e1cef30dd5483cd4cd8c82..482c5b3a004e1521d454259a68a963a57f652c9c 100644
--- a/R/prepare_data2fit.R
+++ b/R/prepare_data2fit.R
@@ -34,7 +34,7 @@ countMatrix_2longDtf <- function(countMatrix, value_name = "kij", id_vars = "gen
 #' @export
 #' @examples
 #' list_var <- init_variable()
-#' mock_data <- mock_rnaseq(list_var, n_genes = 3, 2,2, 2)
+#' mock_data <- mock_rnaseq(list_var, n_genes = 3, 2,2)
 #' dtf_countLong <- countMatrix_2longDtf(mock_data$counts)
 #' getColumnWithSampleID(dtf_countLong, mock_data$metadata)
 getColumnWithSampleID <- function(dtf_countsLong, metadata) {
@@ -67,7 +67,7 @@ getColumnWithSampleID <- function(dtf_countsLong, metadata) {
 #' @export
 #' @examples
 #'  list_var <- init_variable()
-#'  mock_data <- mock_rnaseq(list_var, n_genes = 3, 2,2, 2)
+#'  mock_data <- mock_rnaseq(list_var, n_genes = 3, 2,2)
 #'  data2fit <- prepareData2fit(mock_data$counts, mock_data$metadata)
 prepareData2fit <- function(countMatrix, metadata, normalization = TRUE , response_name = "kij", groupID = "geneID" ) {
   
diff --git a/man/eval_identityTerm.Rd b/man/eval_identityTerm.Rd
index 07a47adcaa398b7f86671c361a339922c26e7d96..faa44035d54e810d6e6f4a044f520c943b7965b4 100644
--- a/man/eval_identityTerm.Rd
+++ b/man/eval_identityTerm.Rd
@@ -16,6 +16,8 @@ eval_identityTerm(
 
 \item{palette_color}{dict-like palette default: palette_color = c(DESeq2 = "#500472", HTRfit ="#79cbb8")}
 
+\item{palette_shape}{Optional parameter that sets the point shape for plots.Default : c(DESeq2 = 17, HTRfit = 19).}
+
 \item{...}{additional parameters to pass geom_point aes}
 }
 \value{
diff --git a/man/getColumnWithSampleID.Rd b/man/getColumnWithSampleID.Rd
index 288ba8714815f8667ac4755f51762cc2064fa49a..9439edb5e3e992d9763709ef4c8a16e2fd187028 100644
--- a/man/getColumnWithSampleID.Rd
+++ b/man/getColumnWithSampleID.Rd
@@ -19,7 +19,7 @@ Returns the column name in the metadata data frame that corresponds to the given
 }
 \examples{
 list_var <- init_variable()
-mock_data <- mock_rnaseq(list_var, n_genes = 3, 2,2, 2)
+mock_data <- mock_rnaseq(list_var, n_genes = 3, 2,2)
 dtf_countLong <- countMatrix_2longDtf(mock_data$counts)
 getColumnWithSampleID(dtf_countLong, mock_data$metadata)
 }
diff --git a/man/prepareData2fit.Rd b/man/prepareData2fit.Rd
index 4158249ff833b4560270f36c1ad8838acc285488..94555f4ade08d256a6491ed05ca62d1e62f847de 100644
--- a/man/prepareData2fit.Rd
+++ b/man/prepareData2fit.Rd
@@ -34,6 +34,6 @@ Prepares the countMatrix and metadata for fitting by converting the countMatrix
 }
 \examples{
  list_var <- init_variable()
- mock_data <- mock_rnaseq(list_var, n_genes = 3, 2,2, 2)
+ mock_data <- mock_rnaseq(list_var, n_genes = 3, 2,2)
  data2fit <- prepareData2fit(mock_data$counts, mock_data$metadata)
 }
diff --git a/pkgdown/_pkgdown.yml b/pkgdown/_pkgdown.yml
index 21b84db3f2a1d02cfd7031b59b4d77f75477a503..58de1b13fd33b42dd1f0dff69292013d6a636886 100644
--- a/pkgdown/_pkgdown.yml
+++ b/pkgdown/_pkgdown.yml
@@ -3,7 +3,7 @@ repo:
     home: https://gitbio.ens-lyon.fr/LBMC/yvertlab/vortex/plasticity_mutation/HTRfit
 navbar:
   structure:
-    left:  [home, intro, reference, articles, additionalPArt ,news]
+    left:  [home, intro, reference, articles ,news]
     right: [search , gitlab]
   components:
     gitlab:
@@ -18,10 +18,6 @@ template:
   includes:
     in_header:
       <script defer data-domain="pkgdown.r-lib.org,all.tidyverse.org" src="https://plausible.io/js/plausible.js"></script>
-
-additionalPArt:
-- title: Testing
-  
 news:
   releases:
   - text: "Version 2.0.0"
diff --git a/tests/testthat/test-prepare_data2fit.R b/tests/testthat/test-prepare_data2fit.R
index b288e616b40823007a427bef587ed4245e6c0719..beb3720b6cba2b6e3fd11cb8b6d0cad45126b1d8 100644
--- a/tests/testthat/test-prepare_data2fit.R
+++ b/tests/testthat/test-prepare_data2fit.R
@@ -21,7 +21,7 @@ test_that("countMatrix_2longDtf converts count matrix to long data frame", {
 test_that("getColumnWithSampleID returns column name with sampleID", {
   # dummy data
   list_var <- init_variable()
-  mock_data <- mock_rnaseq(list_var, n_genes = 3, 2,2, 2)
+  mock_data <- mock_rnaseq(list_var, n_genes = 3, 2,2)
   dtf_countLong <- countMatrix_2longDtf(mock_data$counts)
   
   # Expected output
@@ -38,7 +38,7 @@ test_that("getColumnWithSampleID returns column name with sampleID", {
 test_that("prepareData2fit prepares data for fitting", {
   # dummy data
   list_var <- init_variable()
-  mock_data <- mock_rnaseq(list_var, n_genes = 3, 2,2, 2)
+  mock_data <- mock_rnaseq(list_var, n_genes = 3, 2,2)
   
   # Prepare data for fitting
   data2fit <- prepareData2fit(mock_data$counts, mock_data$metadata)