Add a group to a simulation output

Add a group variable to a simulation from stics_eval().

add_group(stics_eval_out, y)

Arguments

stics_eval_out	The output of a stics_eval simulation
y	A vector of the group value, either of length 1 or nrow(stics_eval_out$outputs)

Details

This is usefull when comparing several simulations, to class outputs per simulations

See also

plot_output() to use this group variable

Examples

if (FALSE) {
library(sticRs)
library(magrittr)

# Make two simulations on two different years, and compare them in the same plot:
# Simulation from 2018 (note the add_group function call at the end)
Eval_stics_1=
   sticRs::stics_eval(dir.orig = "0-DATA/dummy/IC_Wheat_Pea_2018",
             dir.targ = "2-Simulations/Parameterization",
             stics = list(EquivDens= "0-DATA/stics_executable/EquDens/stics.exe"),
             obs_name =  c("6_IC_Wheat_N02.obs","6_IC_Pea_N02.obs"),
             Out_var = c("lai(n)", "masec(n)","hauteur"),
             Title = "Wheat-Pea Auzeville 2005-2006 N0", plot_it = T)%>%
             add_group("2017-2018")

# Simulation from 2017 (note the add_group function call at the end):
Eval_stics_2=
   sticRs::stics_eval(dir.orig = "0-DATA/dummy/IC_Wheat_Pea_2017",
             dir.targ = "2-Simulations/Parameterization",
             stics = list(EquivDens= "0-DATA/stics_executable/EquDens/stics.exe"),
             obs_name =  c("6_IC_Wheat_N02.obs","6_IC_Pea_N02.obs"),
             Out_var = c("lai(n)", "masec(n)","hauteur"),
             Title = "Wheat-Pea Auzeville 2005-2006 N0", plot_it = T)%>%
             add_group("2016-2017")

# Binding the two simulations:
Eval_stics= sticRs::rbind_sim(Eval_stics_1, Eval_stics_2)

# plotting the results:
plot_output(Eval_stics$outputs)

}