Module2.md
Environmental Characterization Module
- Environmental characterization is a fundamental step to understand how the environment regulates the phenotypic expression and adaptation of the genotypes under different growing conditions. For this reason, based on envirotyping (environmental + typing) studies alredy published (see references), we provide mechanisms that enable the typing of environmental factors in terms of frequency of occurrence. We have also developed functions for collecting environmental factors and organizing them as covariates to be used in reaction norm studies.
### Environmental Covariables (ECs, W matrix)
> * Environmental variables can be used as indicators of the quality of an environment (experiment, location). A double entry table (*q* environments x *k* environmental factors) can be built as suggested by Jarquin et al (2014). Hereafther we will refer to this matrix as **W**, and therefore, it will be obtained by the function *W_matrix*:
```{r}
W_matrix(env.data = df.clim,by.interval = F)
wzxhzdk:0
> * Different statistics can be used, as in *summaryWTH*. The statistic argument is used to select between *mean* or *quantile*. If the selection is made in quantile, the complementary argument *prob* is used to choose the quantiles to be used.
wzxhzdk:1
> * We can perform a Quality Control (QC) based on the maximum sd tolered.
wzxhzdk:2
> * We can perform a Quality Control (QC) based on the maximum sd tolered.
wzxhzdk:3
> * Create for specific variables. To do this, insert the name of the variables in the *id.var* argument.
wzxhzdk:4
> * Or even combine with summaryWTH by using *is.processed=TRUE*:
```{r}
data<-summaryWTH(df.clim,env.id = 'env',statistic = 'quantile')
W_matrix(env.data = data,is.processed = T)
wzxhzdk:5
> * Typologies can be defined across different time intervals by setting the argument *by.interval = TRUE* (generic time intervals)
wzxhzdk:6
> * Typologies can be defined across **specific** time intervals by setting the argument *by.interval = TRUE* and and defining the time windows (in days after begining of the data) using the function *time.window*
```{r}
env_typing(env.data = df.clim,env.id = 'env',var.id='T2M',by.interval = T,time.window = c(0,15,35,65,90,120))
wzxhzdk:7
> * The two-way table of typologies can be plotted based on this [code](https://raw.githubusercontent.com/allogamous/EnvRtype/master/plot.R)
### **Example 1 : Air Temperature (**T2M**)**
**Option 1: facet by developmental stages**
