diff --git a/doc/Documentation.pdf b/doc/Documentation.pdf
index e412a80d343be1e8e06f9ca1362af6d9d3058c22..58638f6bb17390dccc273a25b8cab26826c6182c 100644
Binary files a/doc/Documentation.pdf and b/doc/Documentation.pdf differ
diff --git a/src/R/4.summary1.R b/src/R/4.summary1.R
index 8327b92ed6ac675c5109ff8ead199db709e13b72..a5f946d8f6de645021818e9ac67d72c5d7171bfb 100755
--- a/src/R/4.summary1.R
+++ b/src/R/4.summary1.R
@@ -154,9 +154,9 @@ mode_peaks = mlv(round(peaks.df$D2_l2FC, 2), method = "mfv", na.rm = TRUE)
 
 summary.df = summary.df %>%
               dplyr::mutate(
-                  adj_pvalue = p.adjust(Ttest_pval, method = "fdr"),
+                  pvalue_adj = p.adjust(pvalue_raw, method = "fdr"),
                   Diff_mean  = Mean_l2FC   - mean  (peaks.df$D2_l2FC, na.rm = TRUE), 
-                  DiffMedian = Median_l2FC - median(peaks.df$D2_l2FC, na.rm = TRUE),
+                  Diff_median = Median_l2FC - median(peaks.df$D2_l2FC, na.rm = TRUE),
                   Diff_mode  = Mode_l2FC - mode_peaks[[1]],    
                   Diff_skew  = Modeskewness - mode_peaks[[2]])  %>%
               na.omit(summary.df)