Skip to content
GitLab
Select Git revision
  • master default
Find file BlameHistoryPermalink
format.R 7.99 KiB
1 
# \\\
2 
# Copyright 2021-2022 Louis Héraut*1
3 
#
4 
# *1   INRAE, France
5 
#      louis.heraut@inrae.fr
6 
#
7 
# This file is part of ash R toolbox.
8 
#
9 
# ash R toolbox is free software: you can redistribute it and/or modify
10 
# it under the terms of the GNU General Public License as published by
11 
# the Free Software Foundation, either version 3 of the License, or (at
12 
# your option) any later version.
13 
#
14 
# ash R toolbox is distributed in the hope that it will be useful, but
15 
# WITHOUT ANY WARRANTY; without even the implied warranty of
16 
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
17 
# General Public License for more details.
18 
#
19 
# You should have received a copy of the GNU General Public License
20 
# along with ash R toolbox.  If not, see <https://www.gnu.org/licenses/>.
21 
# ///
22 
#
23 
#
24 
# processing/format.R
25 
#
26 
# Manages all the format problem of data and info. Mainly problem of
27 
# input and output of the 'StatsAnalysisTrend' package. It also allows
28 
# to join different selections of station and to gets exact period of
29 
# trend analysis.
30
31
32 
# Usefull library
33 
library(dplyr)
34
35
36 
## 1. INPUT
37 
### 1.1. Preparation
38 
# Prepares the data in order to have a list of a data tibble with
39 
# date, group and flow column and a info tibble with the station code
40 
# and group column to fit the entry of the 'extract.Var' function in
41 
# the 'StatsAnalysisTrend' package
42 
prepare = function(df_data, colnamegroup=NULL) {
43
44 
    # Forces the column name to group to be a vector
45 
    colnamegroup = c(colnamegroup)
46 
    # Converts it to index of the column to group
47 
    colindgroup = which(colnames(df_data) == colnamegroup)
48 
    # Groups the data by those indexes
49 
    df_data = group_by_at(df_data, colindgroup)
50
51 
    # Creates a new tibble of data with a group column
52 
    data = tibble(Date=df_data$Date,
53 
                  group=group_indices(df_data),
54 
                  Qm3s=df_data$Qm3s)
55
56 
    # Gets the different value of the group
57 
    Gkey = group_keys(df_data)
58 
    # Creates a new tibble of info of the group
59 
    info = bind_cols(group=seq(1:nrow(Gkey)),
60 
                     Gkey)
61
62 
    # Stores data and info tibble as a list that match the entry of
63 
    # the 'extract.Var' function
64 
    res = list(data=data, info=info)
65 
    return (res)
66 
}
67
68 
### 1.2. Re-preparation
69 
# Re-prepares the data in outing of the 'extract.Var' function in
70 
# the 'StatsAnalysisTrend' package in order to fit again to the