Parcial_2_3_Valentina_Ortiz

# Rendimiento por bloque y tratamiento
rendimientos <- c(
  5.5, 5.5, 6.9, 7.2,   # Bloque 1: C A B D
  5.8, 5.3, 5.6, 4.1,   # Bloque 2: D B C A
  6.9, 5.3, 6.6, 4.5,   # Bloque 3: D A B C
  5.0, 7.2, 6.1, 7.0,   # Bloque 4: A B C D
  6.2, 5.7, 4.8, 5.8    # Bloque 5: B C A D
)

# Tratamientos según el orden de los datos (respetando los bloques)
tratamientos <- c(
  "C", "A", "B", "D",
  "D", "B", "C", "A",
  "D", "A", "B", "C",
  "A", "B", "C", "D",
  "B", "C", "A", "D"
)

# Bloques: 5 bloques con 4 tratamientos cada uno
bloques <- factor(rep(1:5, each = 4))

# Convertimos tratamiento a factor
tratamientos <- factor(tratamientos)

# Crear dataframe
datos <- data.frame(
  Rendimiento = rendimientos,
  Tratamiento = tratamientos,
  Bloque = bloques
)

=============================

HIPÓTESIS ESTADÍSTICAS

=============================

H0: Todos los sistemas de cultivo (tratamientos) tienen el mismo efecto sobre el rendimiento

H1: Al menos un sistema de cultivo tiene un efecto diferente

# =============================
# ANOVA CON BLOQUES
# =============================

modelo <- aov(Rendimiento ~ Tratamiento + Bloque, data = datos)
summary(modelo)

##             Df Sum Sq Mean Sq F value   Pr(>F)    
## Tratamiento  3  8.946  2.9820  11.528 0.000756 ***
## Bloque       4  3.520  0.8800   3.402 0.044346 *  
## Residuals   12  3.104  0.2587                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

# =============================
# PRUEBA POST ANOVA: LSD
# =============================

# Si no lo has hecho antes:
# install.packages("agricolae")

library(agricolae)

# LSD con bloque incluido en el modelo
lsd_result <- LSD.test(modelo, "Tratamiento", group = TRUE)

# Mostrar los grupos formados
print(lsd_result$groups)

##   Rendimiento groups
## D        6.54      a
## B        6.44      a
## C        5.48      b
## A        4.94      b

# =============================
# GRÁFICO DE BARRAS CON GRUPOS LSD
# =============================

bar.group(lsd_result$groups,
          ylim = c(0, 8),
          col = c("lightblue", "lightgreen", "orange", "salmon"),
          main = "Rendimiento de col crespa por sistema de cultivo",
          xlab = "Sistema de cultivo",
          ylab = "Rendimiento (t/ha)")

=============================

CONCLUSIÓN DEL ENSAYO

=============================

Según la prueba LSD, el sistema de cultivo D presentó el mayor rendimiento

y fue estadísticamente superior al sistema A.

Se sugiere considerar D como opción recomendada en condiciones similares.