Resultados de los tratamientos para el estrés hídrico en rábano - Diseño factorial completo en arreglo completamente al azar

library(dplyr)

## 
## Attaching package: 'dplyr'

## The following objects are masked from 'package:stats':
## 
##     filter, lag

## The following objects are masked from 'package:base':
## 
##     intersect, setdiff, setequal, union

Tratamiento <- c("1", "2", "3", "4")

Riego <- c("90", "60", "60", "60")

Vit.C <- c("0", "0", "500", "1000")

df <- data.frame(Tratamiento, Riego, Vit.C)

df %>% rename ("Riego (%)" = Riego, "Vitamina C (ppm)" = Vit.C)

# Cultivo de rábano.
set.seed(120)

# Factor 1.
R1 <- gl(1, 30, labels = c("90%"))
R2 <- gl(1, 90, labels = c("60%"))


Riego2 <- c(R1, R2)

# Factor 2.
VC1 <- gl(1, 60, labels = c("0ppm"))
VC2 <- gl(2, 30, labels = c("500ppm", "1000ppm"))


Vit.C2 <- c(VC1, VC2)

# Respuesta ante el estrés hídrico.

# Contenido relativo de agua.
CRA <- rnorm(n=120, mean = 3, sd = 0.3)

# No. de estomas abiertos.
No.EA <- rnorm(n=120, mean = 3, sd = 0.3)

# No. de estomas cerrados.
No.EC <- rnorm(n=120, mean = 3, sd = 0.3)

# Contenido de clorofila.
Cont.Clor <- rnorm(n=120, mean = 3, sd = 0.3)

# Altura del follaje.
Alt.F <- rnorm(n=120, mean = 3, sd = 0.3)

# Peso fresco.
PF <- rnorm(n=120, mean = 3, sd = 0.3)

# Peso seco.
PS <- rnorm(n=120, mean = 3, sd = 0.3)

# Area Foliar.
AF <- rnorm(n=120, mean = 3, sd = 0.3)

# Pérdida de electrolitos.
PE <- rnorm(n=120, mean = 3, sd = 0.3)

Tratamientos <- data.frame(Riego2, Vit.C2, CRA, No.EA, No.EC, Cont.Clor, Alt.F, PF, PS, AF, PE); Tratamientos

library(collapsibleTree)
collapsibleTreeSummary(df = Tratamientos, hierarchy = c('Riego2', 'Vit.C2', 'CRA'), collapsed = F, fontSize = 16)