Parcial 2: ## Factorial Simple al Azar. ## Factorial Completo en Arreglo Completamente al Azar. ## Factorial Completo en Bloques. ______________________________________________________ ## Factorial Completo en Parcelas Divididas. ## Factorial Incompleto al Azar.

Trabajo Final (INVESTIGAR). ## Factorial Incompleto en Bloques con Covariables ## Script Block Design.

Taller - Final (Quiz): ## Diseño en Medidas Repetidas. ## Diseño en Lattice. ## Analisis de Perfiles. ## Diseño de Efectos Aleatorios y Mixtos. ____________________________________________________________________________

#Factorial Completo en Parcelas Divididas. (Split Plot Design) - De los mas comunes en Agricultura. Deben existir 2 Factores como minimo.

Similar a hacer 2 veces el experimento, ya que requiere 2 aleatorizaciones. - Riego - Fertilizacion - Densidad de Siembra - Maquinaria – Como cada parcela es diferente, se recomienda usarla como un Factor de Bloqueo. (Variedad Espacial, no controlables).

# 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
xy = expand.grid(x = seq(6), y = seq(4))
xy = sample_frac(xy)
xy$f1 = gl(4, 6, 24, paste0('V',1:4))
xy$f2 = gl(3, 2, 24, paste0('R',1:3))
xy$rep = gl(2, 1, 24, paste0('r',1:2))
xy$name = paste0(xy$f1, xy$f2, xy$rep)
xy
##    x y f1 f2 rep   name
## 1  2 2 V1 R1  r1 V1R1r1
## 2  2 4 V1 R1  r2 V1R1r2
## 3  3 4 V1 R2  r1 V1R2r1
## 4  4 1 V1 R2  r2 V1R2r2
## 5  1 4 V1 R3  r1 V1R3r1
## 6  1 3 V1 R3  r2 V1R3r2
## 7  5 2 V2 R1  r1 V2R1r1
## 8  4 3 V2 R1  r2 V2R1r2
## 9  4 2 V2 R2  r1 V2R2r1
## 10 2 3 V2 R2  r2 V2R2r2
## 11 6 2 V2 R3  r1 V2R3r1
## 12 5 3 V2 R3  r2 V2R3r2
## 13 3 2 V3 R1  r1 V3R1r1
## 14 3 3 V3 R1  r2 V3R1r2
## 15 1 1 V3 R2  r1 V3R2r1
## 16 6 4 V3 R2  r2 V3R2r2
## 17 3 1 V3 R3  r1 V3R3r1
## 18 5 1 V3 R3  r2 V3R3r2
## 19 1 2 V4 R1  r1 V4R1r1
## 20 2 1 V4 R1  r2 V4R1r2
## 21 6 1 V4 R2  r1 V4R2r1
## 22 4 4 V4 R2  r2 V4R2r2
## 23 6 3 V4 R3  r1 V4R3r1
## 24 5 4 V4 R3  r2 V4R3r2
library(ggplot2)
ggplot(xy)+
  aes(x,y,label=name, fill=f2)+
  geom_tile(color='white')+
  geom_label(fill='white')+
  labs(title = 'Completamente al Azar')

No es posible Aleatorizar el riego, por mera disposicion de las tuberias en riego. - Requiere una doble aleatoriacion. Se genera un segundo Factor con la Variedad (Subparcelas).

Se organiza el primer factor de forma consistente, y se aleatoriza el segundo.

library(dplyr)
xy = expand.grid(y = seq(4), x = seq(6))
f2 = gl(3, 8, 24, paste0('R',1:3))
lf1 = paste0('V',1:4)
f1 = c(sample(lf1),sample(lf1),
       sample(lf1),sample(lf1),
       sample(lf1),sample(lf1))
rep = rep(rep(paste0('r',1:2), each=4), 3)

data = data.frame(xy, f1, f2, rep)
data$name = with(data, paste0(f1, rep))
library(ggplot2)
ggplot(data)+
  aes(x,y,label=name, fill=f1)+
  geom_tile(color='white')+
  geom_text()+
  facet_wrap(~f2, scales = 'free')+
  theme(axis.text = element_blank())

En Parcelas divididas, se suman las aleatorizaciones mas los errores de cada repeticion (2), y la interaccion entre ambos.

set.seed(123)
data$biom= rnorm(24,8,2)

Analisis Descriptivo.

ggplot(data)+
  aes(f2,biom)+
  geom_boxplot()

ggplot(data)+
  aes(f1,biom)+
  geom_boxplot()

ggplot(data)+
  aes(f2,biom,fill=f1)+
  geom_boxplot()

mod1= aov(biom ~ f2*f1 + Error(f2:rep),data)
## Warning in aov(biom ~ f2 * f1 + Error(f2:rep), data): Error() model is singular
summary(mod1)
## 
## Error: f2:rep
##           Df Sum Sq Mean Sq F value Pr(>F)  
## f2         2 13.327   6.664   8.078  0.062 .
## Residuals  3  2.475   0.825                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Error: Within
##           Df Sum Sq Mean Sq F value Pr(>F)  
## f1         3  27.59   9.196   3.676 0.0562 .
## f2:f1      6  18.68   3.113   1.244 0.3685  
## Residuals  9  22.52   2.502                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

No se Rechaza la Hipotesis Nula (Que dice que NO hay interaccion), por ende NO existe relacion entre la biomasa y el riego.