Análisis de datos

Tabla de contenido

1. Grupo 1
2. Grupo 2
3. Grupo 3
4. Grupo 4
5. Grupo 5
6. Grupo 6
7. Grupo 7
8. Grupo 8

Grupo 1

MuestreoGrupo1Apio <- read_excel("/home/kevin/Trabajos/Manejos/Muestreo1/MuestreoGrupo1.xlsx", sheet = "Apio")
MuestreoGrupo1Avena <- read_excel("/home/kevin/Trabajos/Manejos/Muestreo1/MuestreoGrupo1.xlsx", sheet = "Avena")

Estadísticas Apio

Datos_Ap1 <- MuestreoGrupo1Apio %>%
          group_by(Unidad_espacial) %>%
          summarise(media_A = round(mean(Altura),3),
                    media_C = round(mean(Cobertura),3),
                    media_S = round(mean(Severidad),3),
                    desv_A = round(sd(Altura),3),
                    desv_C = round(sd(Cobertura),3),
                    desv_S = round(sd(Severidad),3),
                    varianza_A = round(var(Altura),3),
                    varianza_C = round(var(Cobertura),3),
                    varianza_S = round(var(Severidad),3),
                    moda_Ins = mode_f(Fenologia),
                    moda_Ins = mode_f(Incidencia),
                    curtosis_A = round(kurtosis(Altura),3),
                    curtosis_C = round(kurtosis(as.numeric(Cobertura)),3),
                    curtosis_S = round(kurtosis(Severidad),3))
datatable(Datos_Ap1, autoHideNavigation = T)

Gráficas

Altura

MuestreoGrupo1Apio %>%
  group_by(Unidad_espacial)%>%
  ggplot(aes(y=Altura, fill = Unidad_espacial))+
  geom_boxplot()+
  facet_wrap(~Unidad_espacial, scale="free")

Violin

MuestreoGrupo1Apio %>%
  group_by(Unidad_espacial)%>%
  ggplot(aes(x=0,y=Altura, fill = Unidad_espacial))+
  geom_violin()+
  facet_wrap(~Unidad_espacial, scale="free")

Cobertura

MuestreoGrupo1Apio %>%
  group_by(Unidad_espacial)%>%
  ggplot(aes(y=Cobertura, fill = Unidad_espacial))+
  geom_boxplot()+
  facet_wrap(~Unidad_espacial, scale="free")

Severidad

MuestreoGrupo1Apio %>%
  group_by(Unidad_espacial)%>%
  ggplot(aes(y=Severidad, fill = Unidad_espacial))+
  geom_boxplot()+
  facet_wrap(~Unidad_espacial, scale="free")

Violin

MuestreoGrupo1Apio %>%
  group_by(Unidad_espacial)%>%
  ggplot(aes(x=0,y=Severidad, fill = Unidad_espacial))+
  geom_violin()+
  facet_wrap(~Unidad_espacial, scale="free")

Estadísticas Avena

Datos_Av1 = MuestreoGrupo1Avena %>%
  summarise(media_a = mean(Avena),
            mediana_L = median(Lolium),
            moda_H = mode_f(Holcus),
            moda_C = mode_f(Cobertura),
            media_D = mean(Dicotiledonea, na.rm = T))
datatable(Datos_Av1, style = 'bootstrap4')

Grupo 2

MuestreoGrupo2Apio <- read_excel("/home/kevin/Trabajos/Manejos/Muestreo1/MuestreoGrupo2.xlsx", sheet = "Apio")
MuestreoGrupo2Avena <- read_excel("/home/kevin/Trabajos/Manejos/Muestreo1/MuestreoGrupo2.xlsx", sheet = "Avena")

Estadísticas Apio

Datos_Ap2 <- MuestreoGrupo2Apio %>%
              summarise(media_A = round(mean(`Altura (cm)`),3),
                        media_C = round(mean(`Cobertura(%)`),3),
                        
                        desv_A = round(sd(`Altura (cm)`),3),
                        desv_C = round(sd(`Cobertura(%)`),3),
                        
                        varianza_A = round(var(`Altura (cm)`),3),
                        varianza_C = round(var(`Cobertura(%)`),3),
                        moda_F = mode_f(Fenología),
                        moda_I = mode_f(Incidencia),
                        curtosis_A = round(kurtosis(`Altura (cm)`),3),
                        curtosis_C = round(kurtosis(`Cobertura(%)`),3))
datatable(Datos_Ap2, style = 'bootstrap4')

Gráficas

Altura

MuestreoGrupo2Apio %>%
  ggplot(aes(x = 0, y=`Altura (cm)`))+
  geom_boxplot(width = 1.5, fill = "brown", varwidth = F)+
  stat_summary(fun.y = mean, color = 'yellow')

MuestreoGrupo2Apio %>%
  ggplot(aes(x = 0, y=`Altura (cm)`))+
  geom_violin(fill="brown")+
  stat_summary(fun.y = mean, color = 'yellow')+
  stat_summary(fun.y = median, color = 'blue')

Cobertura

MuestreoGrupo2Apio %>%
  ggplot(aes(x = 0, y=`Cobertura(%)`))+
  geom_boxplot(width = 1.5, fill = "brown", varwidth = F)+
  stat_summary(fun = mean, color = 'yellow')

MuestreoGrupo2Apio %>%
  ggplot(aes(x = 0, y=`Cobertura(%)`))+
  geom_violin(fill="brown")+
  stat_summary(fun = mean, color = 'yellow')+
  stat_summary(fun = median, color = 'blue')

Estadísticas Avena

Datos_Av2 <- MuestreoGrupo2Avena %>%
              summarise(media_L = round(mean(`Lollium (%)`, na.rm = T),3),
                        media_C = round(mean(`Cobertura(%)`),3),
                        
                        desv_L = round(sd(`Lollium (%)`, na.rm=T),3),
                        desv_C = round(sd(`Cobertura(%)`),3),
                        
                        varianza_L = round(var(`Lollium (%)`, na.rm = T),3),
                        varianza_C = round(var(`Cobertura(%)`),3),
                        
                        curtosis_L = round(kurtosis(`Lollium (%)`, na.rm=T),3),
                        curtosis_C = round(kurtosis(`Cobertura(%)`),3))
datatable(Datos_Av2, style = 'bootstrap4')

Gráficas

MuestreoGrupo2Avena %>%
  ggplot(aes(x = 0, y = `Cobertura(%)`))+
  geom_boxplot(fill = "cyan4")+
  stat_summary(fun = mean, color = "red")

MuestreoGrupo2Avena %>%
  ggplot(aes(x = 0, y = `Lollium (%)`))+
  geom_boxplot(fill = "cyan4")+
  stat_summary(fun = mean, color = "red")

MuestreoGrupo2Avena %>%
  ggplot(aes(x = 0, y = `Cobertura(%)`))+
  geom_violin(fill = "cyan4")+
  stat_summary(fun = mean, color = "red")+
  stat_summary(fun = median, color = "yellow")

MuestreoGrupo2Avena %>%
  ggplot(aes(x = 0, y = `Lollium (%)`))+
  geom_violin(fill = "cyan4")+
  stat_summary(fun = mean, color = "red")+
  stat_summary(fun = median, color = "yellow")

Grupo 3

MuestreoGrupo3Apio <- read_excel("/home/kevin/Trabajos/Manejos/Muestreo1/MuestreoGrupo3.xlsx", sheet = "Apio")
MuestreoGrupo3Avena <- read_excel("/home/kevin/Trabajos/Manejos/Muestreo1/MuestreoGrupo3.xlsx", sheet = "Avena")

Estadística Apio

Datos_Ap3 <- MuestreoGrupo3Apio %>%
  group_by(Estrato, Muestra) %>%
  summarise(media_A = mean(`Altura (cm)`),
            media_S = mean(`Severidad (%)`),
            media_C = mean(`Cobertura de apio`),
            media_M = mean(`Cobertura Malezas (%)`),
            desv_A = sd(`Altura (cm)`),
            desv_S = sd(`Severidad (%)`),
            desv_C = sd(`Cobertura de apio`),
            desv_M = sd(`Cobertura Malezas (%)`),
            var_A = var(`Altura (cm)`),
            var_S = var(`Severidad (%)`),
            var_C = var(`Cobertura de apio`),
            var_M = var(`Cobertura Malezas (%)`),
            mode_Fe = mode_f(Fenología))

## `summarise()` has grouped output by 'Estrato'. You can override using the `.groups` argument.

datatable(Datos_Ap3, style = 'bootstrap4')

Gráficas

MuestreoGrupo3Apio %>%
  group_by(Estrato)%>%
  ggplot(aes(x= 0, y=`Altura (cm)`, fill = Estrato))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~Estrato, scale="free")

MuestreoGrupo3Apio %>%
  group_by(Estrato)%>%
  ggplot(aes(x=0,y=`Severidad (%)`, fill = Estrato))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~Estrato, scale="free")

MuestreoGrupo3Apio %>%
  group_by(Estrato)%>%
  ggplot(aes(x = 0, y=`Cobertura de apio`, fill = Estrato))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~Estrato, scale="free")

MuestreoGrupo3Apio %>%
  group_by(Estrato)%>%
  ggplot(aes(x = 0, y=`Cobertura Malezas (%)`, fill = Estrato))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~Estrato, scale="free")

MuestreoGrupo3Apio %>%
  group_by(Estrato)%>%
  ggplot(aes(x= 0, y=`Altura (cm)`, fill = Estrato))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~Estrato, scale="free")

MuestreoGrupo3Apio %>%
  group_by(Estrato)%>%
  ggplot(aes(x= 0, y=`Severidad (%)`, fill = Estrato))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~Estrato, scale="free")

MuestreoGrupo3Apio %>%
  group_by(Estrato)%>%
  ggplot(aes(x= 0, y=`Cobertura de apio`, fill = Estrato))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~Estrato, scale="free")

MuestreoGrupo3Apio %>%
  group_by(Estrato)%>%
  ggplot(aes(x= 0, y=`Cobertura Malezas (%)`, fill = Estrato))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~Estrato, scale="free")

Estadística Avena

Datos_Av3 <- MuestreoGrupo3Avena %>%
  group_by(Estrato) %>%
  summarise(media_A = mean(`Avena (%)`),
            media_S = mean(`Dicoti (%)`),
            media_C = mean(`Lolium (%)`),
            media_M = mean(`Pennisetum (%)`),
            desv_A = round(sd(`Avena (%)`),2),
            desv_S = round(sd(`Dicoti (%)`),2),
            desv_C = round(sd(`Lolium (%)`),2),
            desv_M = round(sd(`Pennisetum (%)`),2),
            var_A = round(var(`Avena (%)`),2),
            var_S = round(var(`Dicoti (%)`),2),
            var_C = round(var(`Lolium (%)`),2),
            var_M = round(var(`Pennisetum (%)`),2))
datatable(Datos_Av3, style = 'bootstrap4')

Gráficas

MuestreoGrupo3Avena %>%
  group_by(Estrato)%>%
  ggplot(aes(x=0,y=`Avena (%)`, fill = Estrato))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~Estrato, scale="free")

MuestreoGrupo3Avena %>%
  group_by(Estrato)%>%
  ggplot(aes(x=0,y=`Dicoti (%)`, fill = Estrato))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~Estrato, scale="free")

MuestreoGrupo3Avena %>%
  group_by(Estrato)%>%
  ggplot(aes(x=0,y=`Lolium (%)`, fill = Estrato))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~Estrato, scale="free")

MuestreoGrupo3Avena %>%
  group_by(Estrato)%>%
  ggplot(aes(x=0,y=`Pennisetum (%)`, fill = Estrato))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~Estrato, scale="free")

MuestreoGrupo3Avena %>%
  group_by(Estrato)%>%
  ggplot(aes(x= 0, y=`Avena (%)`, fill = Estrato))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~Estrato, scale="free")

MuestreoGrupo3Avena %>%
  group_by(Estrato)%>%
  ggplot(aes(x= 0, y=`Dicoti (%)`, fill = Estrato))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~Estrato, scale="free")

MuestreoGrupo3Avena %>%
  group_by(Estrato)%>%
  ggplot(aes(x= 0, y=`Lolium (%)`, fill = Estrato))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~Estrato, scale="free")

MuestreoGrupo3Avena %>%
  group_by(Estrato)%>%
  ggplot(aes(x= 0, y=`Pennisetum (%)`, fill = Estrato))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~Estrato, scale="free")

Grupo 4

MuestreoGrupo4Apio <- read_excel("/home/kevin/Trabajos/Manejos/Muestreo1/MuestreoGrupo4.xlsx", sheet = "Apio")
MuestreoGrupo4Avena <- read_excel("/home/kevin/Trabajos/Manejos/Muestreo1/MuestreoGrupo4.xlsx", sheet = "Avena")

Estadística Apio

Datos_Ap4 = MuestreoGrupo4Apio %>%
  group_by(Punto) %>%
  summarise(media_I = mean(`Incidencia (%)`),
            media_S = mean(`Severidad (%)`),
            media_A = mean(`Altura (cm)`),
            desv_I = round(sd(`Incidencia (%)`),2),
            desv_S = round(sd(`Severidad (%)`),2),
            des_A = round(sd(`Altura (cm)`),2),
            var_I = round(var(`Incidencia (%)`),2),
            var_S = round(var(`Severidad (%)`),2),
            var_A = round(var(`Altura (cm)`),2),
            cur_A = round(kurtosis(`Incidencia (%)`),2),
            cur_S = round(kurtosis(`Severidad (%)`),2),
            cur_A = round(kurtosis(`Altura (cm)`),2),
            moda_H = mode_f(`Número de Hojas`),
            moda_M = mode_f(Malezas))
datatable(Datos_Ap4, style = 'bootstrap4')

Gráficas

MuestreoGrupo4Apio %>%
  group_by(Punto)%>%
  ggplot(aes(x=0,y=`Incidencia (%)`, fill = Punto))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~Punto, scale="free")

MuestreoGrupo4Apio %>%
  group_by(Punto)%>%
  ggplot(aes(x=0,y=`Severidad (%)`, fill = Punto))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~Punto, scale="free")

MuestreoGrupo4Apio %>%
  group_by(Punto)%>%
  ggplot(aes(x=0,y=`Altura (cm)`, fill = Punto))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~Punto, scale="free")

MuestreoGrupo4Apio%>%
  group_by(Punto)%>%
  ggplot(aes(x= 0, y=`Severidad (%)`, fill = Punto))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~Punto, scale="free")

MuestreoGrupo4Apio%>%
  group_by(Punto)%>%
  ggplot(aes(x= 0, y=`Altura (cm)`, fill = Punto))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~Punto, scale="free")

Estadística Avena

Datos_Av4 = MuestreoGrupo4Avena %>%
  summarise(moda_M = mode_f(Malezas),
            moda_A = mode_f(Avena))
datatable(Datos_Av4, style = 'bootstrap4')

Grupo 5

MuestreoGrupo5Apio <- read_excel("/home/kevin/Trabajos/Manejos/Muestreo1/MuestreoGrupo5.xlsx", sheet = "Apio")
MuestreoGrupo5Avena <- read_excel("/home/kevin/Trabajos/Manejos/Muestreo1/MuestreoGrupo5.xlsx", sheet = "Avena")

## New names:
## * `` -> ...2
## * `` -> ...4
## * `` -> ...5
## * `` -> ...6
## * `` -> ...7
## * ...

Estadística Apio

Datos_Ap5 = MuestreoGrupo5Apio %>%
  group_by(Muestra)%>%
  summarise(media_A = mean(altura),
            desv_A = sd(altura),
            var_A = var(altura),
            moda_S = mode_f(severidad),
            moda_ES = mode_f(`Escala de severidad`),
            moda_I = mode_f(incidencia),
            moda_AC = mode_f(`Abierto/ cerrado`),
            moda_E = mode_f(`Escala de cobertura calles`),
            moda_P = mode_f(plagas),
            moda_EP = mode_f(`Escala Plagas`),
            moda_F = mode_f(`fenología (tallos)`))
datatable(Datos_Ap5, style = 'bootstrap4')

Gráficas

MuestreoGrupo5Apio %>%
  group_by(Muestra)%>%
  ggplot(aes(x=0,y=altura, fill = Muestra))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~Muestra, scale="free")

MuestreoGrupo5Apio%>%
  group_by(Muestra)%>%
  ggplot(aes(x= 0, y=altura, fill = Muestra))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~Muestra, scale="free")

Grupo 6

MuestreoGrupo6Apio <- read_excel("/home/kevin/Trabajos/Manejos/Muestreo1/MuestreoGrupo6.xlsx", sheet = "apio 20")
MuestreoGrupo6Avena <- read_excel("/home/kevin/Trabajos/Manejos/Muestreo1/MuestreoGrupo6.xlsx", sheet = "avena")

Estadística Apio

Datos_Ap6 = MuestreoGrupo6Apio %>%
  group_by(muestreo) %>%
  summarise(moda_H = mode_f(hojas),
            moda_E = mode_f(enfermas),
            media_A = round(mean(altura),2),
            media_S = round(mean(severidad),2),
            media_M = round(mean(malezas),2),
            desv_A = round(sd(altura),2),
            desv_S = round(sd(severidad),2),
            desv_M = round(sd(malezas),2),
            var_A = round(var(altura),2),
            var_S = round(var(severidad),2),
            var_M = round(var(malezas),2),
            cur_A = round(kurtosis(altura),2),
            cur_S = round(kurtosis(severidad),2),
            cur_M = round(kurtosis(malezas),2))
datatable(Datos_Ap6, style = 'bootstrap4')

Gráficas

MuestreoGrupo6Apio %>%
  group_by(muestreo)%>%
  ggplot(aes(x=0,y=altura, fill = muestreo))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~muestreo, scale="free")

MuestreoGrupo6Apio %>%
  group_by(muestreo)%>%
  ggplot(aes(x=0,y=severidad, fill = muestreo))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~muestreo, scale="free")

MuestreoGrupo6Apio %>%
  group_by(muestreo)%>%
  ggplot(aes(x=0,y=malezas, fill = muestreo))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~muestreo, scale="free")

MuestreoGrupo6Apio%>%
  group_by(muestreo)%>%
  ggplot(aes(x= 0, y=altura, fill = muestreo))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~muestreo, scale="free")

MuestreoGrupo6Apio%>%
  group_by(muestreo)%>%
  ggplot(aes(x= 0, y=severidad, fill = muestreo))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~muestreo, scale="free")

MuestreoGrupo6Apio%>%
  group_by(muestreo)%>%
  ggplot(aes(x= 0, y=malezas, fill = muestreo))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~muestreo, scale="free")

Estadística Avena

Datos_Av6 = MuestreoGrupo6Avena %>%
  group_by(muestreo) %>%
  summarise(moda_H = mode_f(hojas),
            moda_I =mode_f(incidencia),
            media_A = round(mean(altura),2),
            media_S = round(mean(malezas),2),
            desv_A = round(sd(altura),2),
            desv_M = round(sd(malezas),2),
            var_A = round(var(altura),2),
            var_M = round(var(malezas),2),
            cur_A = round(kurtosis(altura),2),
            cur_M = round(kurtosis(malezas),2))
datatable(Datos_Ap6, style = 'bootstrap4')

Gráficas

MuestreoGrupo6Avena %>%
  group_by(muestreo)%>%
  ggplot(aes(x=0,y=altura, fill = muestreo))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~muestreo, scale="free")

MuestreoGrupo6Avena %>%
  group_by(muestreo)%>%
  ggplot(aes(x=0,y=malezas, fill = muestreo))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~muestreo, scale="free")

MuestreoGrupo6Avena%>%
  group_by(muestreo)%>%
  ggplot(aes(x= 0, y=altura, fill = muestreo))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~muestreo, scale="free")

MuestreoGrupo6Avena%>%
  group_by(muestreo)%>%
  ggplot(aes(x= 0, y=malezas, fill = muestreo))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~muestreo, scale="free")

Grupo 7

MuestreoGrupo7Apio <- read_excel("/home/kevin/Trabajos/Manejos/Muestreo1/MuestreoGrupo7.xlsx", sheet = "DataApio")
MuestreoGrupo7Avena <- read_excel("/home/kevin/Trabajos/Manejos/Muestreo1/MuestreoGrupo7.xlsx", sheet = "DataAvena")

Estadística Apio

Datos_Ap7 = MuestreoGrupo7Apio %>%
  group_by(ESTRATO) %>%
  summarise(media_A = round(mean(ALTURACM),2),
            media_H = round(mean(NO_HOJAS),2),
            media_HAF = round(mean(NO_HOJAS_AF),2),
            media_I = round(mean(INCIDENCIA),2),
            media_C = round(mean(PO_COBERTURA),2),
            desv_A = round(sd(ALTURACM),2),
            desv_H = round(sd(NO_HOJAS),2),
            desv_HAF = round(sd(NO_HOJAS_AF),2),
            desv_I = round(sd(INCIDENCIA),2),
            desv_C = round(sd(PO_COBERTURA),2),
            var_A = round(var(ALTURACM),2),
            var_H = round(var(NO_HOJAS),2),
            var_HAF = round(var(NO_HOJAS_AF),2),
            var_I = round(var(INCIDENCIA),2),
            var_C = round(var(PO_COBERTURA),2),
            cur_A = round(kurtosis(ALTURACM),2),
            cur_H = round(kurtosis(NO_HOJAS),2),
            cur_HAF = round(kurtosis(NO_HOJAS_AF),2),
            cur_I = round(kurtosis(INCIDENCIA),2),
            cur_C = round(kurtosis(PO_COBERTURA),2),
            moda_S = mode_f(SEVERIDAD))
datatable(Datos_Ap7, style = 'bootstrap4')

Gráficas

MuestreoGrupo7Apio %>%
  group_by(ESTRATO)%>%
  ggplot(aes(x=0,y=ALTURACM, fill = ESTRATO))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Apio %>%
  group_by(ESTRATO)%>%
  ggplot(aes(x=0,y=NO_HOJAS, fill = ESTRATO))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Apio %>%
  group_by(ESTRATO)%>%
  ggplot(aes(x=0,y=NO_HOJAS_AF, fill = ESTRATO))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Apio %>%
  group_by(ESTRATO)%>%
  ggplot(aes(x=0,y=INCIDENCIA, fill = ESTRATO))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Apio %>%
  group_by(ESTRATO)%>%
  ggplot(aes(x=0,y=PO_COBERTURA, fill = ESTRATO))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Apio%>%
  group_by(ESTRATO)%>%
  ggplot(aes(x= 0, y=ALTURACM, fill = ESTRATO))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Apio%>%
  group_by(ESTRATO)%>%
  ggplot(aes(x= 0, y=NO_HOJAS, fill = ESTRATO))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Apio%>%
  group_by(ESTRATO)%>%
  ggplot(aes(x= 0, y=NO_HOJAS_AF, fill = ESTRATO))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Apio%>%
  group_by(ESTRATO)%>%
  ggplot(aes(x= 0, y=INCIDENCIA, fill = ESTRATO))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Apio%>%
  group_by(ESTRATO)%>%
  ggplot(aes(x= 0, y=PO_COBERTURA, fill = ESTRATO))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~ESTRATO, scale="free")

Estadística Avena

Datos_Av7 = MuestreoGrupo7Avena %>%
  group_by(ESTRATO) %>%
  summarise(media_A = round(mean(PO_COBERTURA_A),2),
            media_LI = round(mean(PO_COBERTURA_LI),2),
            media_AV = round(mean(Altura_Avena),2),
            media_AM = round(mean(ALTURA_MALEZA),2),
            desv_A = round(sd(PO_COBERTURA_A),2),
            desv_LI = round(sd(PO_COBERTURA_LI),2),
            desv_AV = round(sd(Altura_Avena),2),
            desv_AM = round(sd(ALTURA_MALEZA),2),
            var_A = round(var(PO_COBERTURA_A),2),
            var_LI = round(var(PO_COBERTURA_LI),2),
            var_AV = round(var(Altura_Avena),2),
            var_AM = round(var(ALTURA_MALEZA),2),
            cur_A = round(kurtosis(PO_COBERTURA_A),2),
            cur_LI = round(kurtosis(PO_COBERTURA_LI),2),
            cur_AV = round(kurtosis(Altura_Avena),2),
            cur_AM = round(kurtosis(ALTURA_MALEZA),2))
datatable(Datos_Av7, style = 'bootstrap4')

Gráficas

MuestreoGrupo7Avena %>%
  group_by(ESTRATO)%>%
  ggplot(aes(x=0,y=PO_COBERTURA_A, fill = ESTRATO))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Avena %>%
  group_by(ESTRATO)%>%
  ggplot(aes(x=0,y=PO_COBERTURA_LI, fill = ESTRATO))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Avena %>%
  group_by(ESTRATO)%>%
  ggplot(aes(x=0,y=Altura_Avena, fill = ESTRATO))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Avena %>%
  group_by(ESTRATO)%>%
  ggplot(aes(x=0,y=ALTURA_MALEZA, fill = ESTRATO))+
  geom_boxplot()+
  stat_summary(fun = mean, color = "red") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Avena%>%
  group_by(ESTRATO)%>%
  ggplot(aes(x= 0, y=PO_COBERTURA_A, fill = ESTRATO))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Avena%>%
  group_by(ESTRATO)%>%
  ggplot(aes(x= 0, y=PO_COBERTURA_LI, fill = ESTRATO))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Avena%>%
  group_by(ESTRATO)%>%
  ggplot(aes(x= 0, y=Altura_Avena, fill = ESTRATO))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~ESTRATO, scale="free")

MuestreoGrupo7Avena%>%
  group_by(ESTRATO)%>%
  ggplot(aes(x= 0, y=ALTURA_MALEZA, fill = ESTRATO))+
  geom_violin()+
  stat_summary(fun = mean, color = "red") +
  stat_summary(fun = median, color = "yellow") +
  facet_wrap(~ESTRATO, scale="free")

Grupo 8