Commentary and debuging of rectangle mean step

script.R

-
-###### A MODIFIER ######
+############ A MODIFIER ############
 
 
 # Path to the data
@@ -33,9 +32,9 @@ filename =
       # "A2250310_HYDRO_QJM.txt"
       # )
 
-    c("O3035210_HYDRO_QJM.txt",
-      "O3011010_HYDRO_QJM.txt",
-      "O1442910_HYDRO_QJM.txt")
+    c("S4214010_HYDRO_QJM.txt",
+      "Q6332510_HYDRO_QJM.txt",
+      "Q7002910_HYDRO_QJM.txt")
 
 
 
@@ -74,7 +73,7 @@ mean_period = list(period1, period2)
 p_thresold = 0.1 #c(0.01, 0.05, 0.1)
 
 
-########################
+####################################
 
 
 # FILE STRUCTURE #
@@ -206,13 +205,9 @@ res_VCN10trend = get_VCN10trend(df_data, df_meta,
 # df_break = get_break(res_VCN10trend$data, df_meta)
 
 # histogram(df_break$Date, df_meta,
-          # breaks=seq(min(df_break$Date),
-                     # max(df_break$Date), "years"),
           # figdir=figdir)
 
 # cumulative(df_break$Date, df_meta, dyear=8,
-          # breaks=seq(min(df_break$Date),
-                     # max(df_break$Date), "years"),
           # figdir=figdir)
 
 # TIME PANEL #
@@ -245,38 +240,6 @@ panels_layout(list(res_QAtrend$data, res_QMNAtrend$data,
               figdir=figdir,
               filename_opt='')
 
-# panels_layout(list(res_QAtrend$data, res_VCN10trend$data), 
-#               layout_matrix=c(1, 2),
-#               df_meta=df_meta, 
-#               df_trend=list(res_QAtrend$trend,
-#                     res_VCN10trend$trend), 
-#               type=list(bquote(Q[A]), bquote(V[CN10])),
-#               missRect=list(TRUE, TRUE),
-#               period=trend_period,
-#               info_header=TRUE,
-#               time_header=df_data,
-#               time_ratio=2, 
-#               var_ratio=5,
-#               figdir=figdir,
-#               filename_opt='')
-
-
-
-# panels_layout(list(res_QMNAtrend$data), 
-#               layout_matrix=c(1),
-#               df_meta=df_meta, 
-#               df_trend=list(res_QMNAtrend$trend), 
-#               type=list(bquote(Q[MNA])),
-#               missRect=list(TRUE),
-#               period=trend_period,
-#               info_header=TRUE,
-#               time_header=df_data,
-#               time_ratio=2, 
-              # var_ratio=5,
-              # figdir=figdir,
-#               filename_opt='')
-
-
 ### /!\ Removed 185 row(s) containing missing values (geom_path) -> remove NA ###