## 3.4 Describe the data

Once the data have been collected, a first step is to describe them with `plot()`

.
Seven types of plot, through the `plot_type`

argument are possible:

- presence abscence matrix that represent the combinaison of germplasm \(\times\) location
- histogramm
- barplot, where sd error are displayed
- boxplot
- interaction
- biplot
- radar
- raster
- map

Then you must choose which factor to represent on the x axis (`x_axis`

argument),
the factor to display in color (`in_col`

argument), and of course the variables to describe (`vec_variables`

argument).

It is possible to tune the number of factor displayed (`nb_parameters_per_plot_x_axis`

and `nb_parameters_per_plot_in_col`

arguments) and the size of the labels regarding biplot and radar (`labels_on`

and `labels_size`

arguments).

Note that descriptive plots can be done based on version within the data set. See section 3.8 formore details.

### 3.4.1 Format the data

Get two data set to look at some examples

```
data("data_model_GxE")
data_model_GxE = format_data_PPBstats(data_model_GxE, type = "data_agro")
```

`## data has been formated for PPBstats functions.`

```
data("data_model_bh_GxE")
data_model_bh_GxE = format_data_PPBstats(data_model_bh_GxE, type = "data_agro")
```

```
## Warning in format_data_PPBstats.data_agro(data): Column "long" is needed to
## get map and not present in data.
```

```
## Warning in format_data_PPBstats.data_agro(data): Column "lat" is needed to
## get map and not present in data.
```

`## data has been formated for PPBstats functions.`

### 3.4.2 presence abscence matrix

The presence absence matrix may be different from experimental design planned because of NA. The plot represents the presence/absence matrix of G \(\times\) E combinations.

```
p = plot(
data_model_GxE, plot_type = "pam",
vec_variables = c("y1", "y2")
)
names(p)
```

`## [1] "y1" "y2"`

`p$y1`

A score of 3 is for a given germplasm replicated three times in a given environement.

```
p = plot(
data_model_bh_GxE, plot_type = "pam",
vec_variables = c("y1", "y2")
)
p$y1
```

Here there are lots of 0 meaning that a lot of germplasm are no in at least two locations. A score of 1 is for a given germplasm in a given location. A score of 2 is for a given germplasm replicated twice in a given location.

### 3.4.3 histogramm

```
p = plot(
data_model_GxE, plot_type = "histogramm",
vec_variables = c("y1", "y2")
)
p$y1
```

`## $`-NA|-NA``

`## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.`

### 3.4.4 barplot

```
p = plot(
data_model_GxE, plot_type = "barplot",
vec_variables = c("y1", "y2"),
x_axis = "germplasm"
)
```

Note that for each element of the following list, there are as many graph as needed with `nb_parameters_per_x_axis`

parameters per graph.

`names(p$y1)`

`## [1] "germplasm-1|-NA" "germplasm-2|-NA" "germplasm-3|-NA" "germplasm-4|-NA"`

`p$y1$`germplasm-1|-NA``

```
p = plot(
data_model_GxE, plot_type = "barplot",
vec_variables = c("y1", "y2"),
x_axis = "germplasm",
in_col = "location"
)
```

Note that for each element of the following list, there are as many graph as needed with `nb_parameters_per_x_axis`

and `nb_parameters_per_in_col`

parameters per graph.

`names(p$y1)`

```
## [1] "germplasm-1|location-1" "germplasm-2|location-1"
## [3] "germplasm-3|location-1" "germplasm-4|location-1"
```

`p$y1$`germplasm-1|location-1``

### 3.4.5 boxplot

```
p = plot(
data_model_GxE, plot_type = "boxplot",
vec_variables = c("y1", "y2"),
x_axis = "germplasm"
)
```

Note that for each element of the following list, there are as many graph as needed with `nb_parameters_per_x_axis`

parameters per graph.

`names(p$y1)`

`## [1] "germplasm-1|-NA" "germplasm-2|-NA" "germplasm-3|-NA" "germplasm-4|-NA"`

`p$y1$`germplasm-1|-NA``

```
p = plot(
data_model_GxE, plot_type = "boxplot",
vec_variables = c("y1", "y2"),
x_axis = "germplasm",
in_col = "location"
)
```

Note that for each element of the following list, there are as many graph as needed with `nb_parameters_per_x_axis`

and `nb_parameters_per_in_col`

parameters per graph.

`names(p$y1)`

```
## [1] "germplasm-1|location-1" "germplasm-2|location-1"
## [3] "germplasm-3|location-1" "germplasm-4|location-1"
```

`p$y1$`germplasm-1|location-1``

### 3.4.6 interaction

```
p = plot(
data_model_GxE, plot_type = "interaction",
vec_variables = c("y1", "y2"),
x_axis = "germplasm",
in_col = "location"
)
```

Note that for each element of the following list, there are as many graph as needed with `nb_parameters_per_x_axis`

and `nb_parameters_per_in_col`

parameters per graph.

`names(p$y1)`

```
## [1] "germplasm-1|location-1" "germplasm-2|location-1"
## [3] "germplasm-3|location-1" "germplasm-4|location-1"
```

`p$y1$`germplasm-1|location-1``

It is also possible to have on the `x_axis`

the date in julian day that have been automatically calculated from `format_data_PPBstats()`

.
Note that this is possible only for `plot_type = "histogramm"`

, `"barplot"`

, `"boxplot"`

and `"interaction"`

.

```
p = plot(
data_model_GxE, plot_type = "interaction",
vec_variables = c("y1", "y2"),
x_axis = "date_julian",
in_col = "location"
)
```

```
## Warning in plot_descriptive_data(x, plot_type, x_axis, in_col,
## vec_variables, : x_axis = "date_julian" is a special feature that will
## display julian day for a given variable automatically calculated from
## format_data_PPBstats().
```

`p$y1$`y1$date_julian-1|location-1``

```
## geom_path: Each group consists of only one observation. Do you need to
## adjust the group aesthetic?
```

### 3.4.7 biplot

```
p = plot(
data_model_GxE, plot_type = "biplot",
vec_variables = c("y1", "y2", "y3"),
in_col = "germplasm", labels_on = "germplasm"
)
```

The name of the list correspond to the pairs of variables displayed.
Note that for each element of the following list, there are as many graph as needed with `nb_parameters_per_in_col`

parameters per graph.

`names(p)`

`## [1] "y1 - y2" "y1 - y3" "y2 - y3"`

`p$`y1 - y2`$`-NA|germplasm-1``

### 3.4.8 radar

```
p = plot(
data_model_GxE, plot_type = "radar",
vec_variables = c("y1", "y2", "y3"),
in_col = "location"
)
p
```

### 3.4.9 raster

Raster plot can be done for factor variables.
Note than when there are no single value for a given `x_axis`

, colums `block`

, `X`

and `Y`

are added in order to have single value.

```
p = plot(
data_model_GxE,
plot_type = "raster",
vec_variables = c("desease", "vigor"),
x_axis = "germplasm"
)
```

```
## Warning in fun_raster(data, vec_variables, x_axis,
## nb_parameters_per_plot_x_axis): There are no single value for each x_axis,
## therefore block, X and Y colums have been added in order to have single
## value.
```

`p$`germplasm-block-X-Y-9|-NA``

### 3.4.10 map

You can display map with location if you have data with latitude and longitude for each location. When using map, do not forget to use credit : Map tiles by Stamen Design, under CC BY 3.0. Data by OpenStreetMap, under ODbL.

```
p = plot(
data_model_GxE, plot_type = "map", labels_on = "location"
)
p$map
```

and add pies for a given variables

```
p = plot(
data_model_GxE, vec_variables = c("y1", "desease"),
plot_type = "map"
)
p$pies_on_map_y1
```

`p$pies_on_map_desease`