Skip to content
GitLab
Select Git revision
  • 191-speed-up-runmodel_lag-when-it-is-used-alone
  • dev default
  • master
  • 171-check-error-direct-call-of-as-data-frame-posixct-is-deprecated-use-as-data-frame-vector-or-as
  • 161-runmodel_lag-proposal-for-a-better-transformation-and-screening-of-the-parameters
  • 165-add-missing-isintstore-argument-in-createcaliboptions
  • 33-merge-airgr-and-airgrplus-functions
  • 115-add-pcict-date-format-managment
  • rmOptions
  • 136-decreased-performance-of-calibration-execution-time
  • 118-plot-outputsmodel-remove-case-sensitivity-to-the-which-argument
  • 96-parallelize-the-grid-screening-step-during-calibration
  • 60-migrate-transfoparam-functions-and-createcaliboptions-to-s3-methods
  • airGRplus
  • CRAN_v1.7.6
  • CRAN_v1.7.5
  • CRAN_v1.7.4
  • CRAN_v1.7.3
  • CRAN_v1.7.2
  • CRAN_v1.7.1
  • CRAN_v1.7.0
  • CRAN_v1.6.12
  • CRAN_v1.6.11
  • CRAN_v1.6.10.4
  • CRAN_v1.6.9.27
  • CRAN_v1.6.9.24
  • CRAN_v1.6.9.23
  • CRAN_v1.6.9.21
  • CRAN_v1.4.3.65
  • CRAN_v1.4.3.60
  • CRAN_v1.4.3.52
  • CRAN_v1.3.2.42
  • CRAN_v1.3.2.41
  • CRAN_v1.3.2.40
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
7172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140
# entry of the same function reprepare = function(df_XEx, df_Xlist, colnamegroup=NULL) { # Changes the column name of the results of the # 'extract.Var' function colnames(df_XEx) = c('Date', 'group', 'Qm3s') # Converts Date column as character df_XEx$Date = as.character(df_XEx$Date) # Takes the first date as example exDate = df_XEx$Date[1] # Finds the number of dash in the date nbt = lengths(regmatches(exDate, gregexpr('-', exDate))) # If there is only one dash if (nbt == 1) { # Converts it to date from a year and a month df_XEx$Date = paste(df_XEx$Date, '01', sep='-') # If there is no dash } else if (nbt == 0) { # Converts it to date from only a year df_XEx$Date = paste(df_XEx$Date, '01', '01', sep='-') # If there is more than 2 dashes } else if (nbt != 2) { # This is not a classical date stop('erreur of date format') } # Recreates the outing of the 'extract.Var' function nicer df_XEx = bind_cols(Date=as.Date(df_XEx$Date, format="%Y-%m-%d"), df_XEx[-1], df_Xlist$info[df_XEx$group, 2:ncol(df_Xlist$info)]) # Prepares the nicer outing df_XlistEx = prepare(df_XEx, colnamegroup=colnamegroup) return (df_XlistEx) } ## 2. OUTPUT # Cleans the trend results of the function 'Estimate.stats' in the # 'StatsAnalysisTrend' package. It adds the station code and the # intercept of the trend to the trend results. Also makes the data # more presentable. clean = function (df_Xtrend, df_XEx, df_Xlist) { # Reprepares the list of data and info in order to be presentable df_Xlist = reprepare(df_XEx, df_Xlist, colnamegroup=c('code')) # Adds a column of station code df_Xlist$data$code = NA # For all the group for (g in df_Xlist$info$group) { # Adds the station code corresponding to each group info df_Xlist$data$code[which(df_Xlist$data$group == g)] = df_Xlist$info$code[df_Xlist$info$group == g] } # Adds the info to trend tibble df_Xtrend = bind_cols(df_Xtrend, df_Xlist$info[df_Xtrend$group1, 2:ncol(df_Xlist$info)]) # Renames the column of group of trend results colnames(df_Xtrend)[1] = 'group' # Adds the intercept value of trend df_Xtrend = get_intercept(df_Xtrend, df_Xlist, unit2day=365.25) # Changes the position of the intercept column df_Xtrend = relocate(df_Xtrend, intercept, .after=trend) # Creates a list of results to return res = list(trend=df_Xtrend, data=df_Xlist$data, info=df_Xlist$info)