Untitled

set.seed(123)

datos = data.frame(
gen = gl(3,4,12, c('g1','g2','g3')),
rep = gl(4,1,12, c('r1','r2','r3','r4')),
rto = c(3.5, 3.6, 4.2,
3.8, 3.9, 3.9,
3.6, 4.1, 4.3,
3.5, 3.8, 4.3)
)

datos

TABLA DE ANALISIS DE VARIANZA

boxplot(rto ~ gen, data=datos)

boxplot(rto ~ rep, data=datos)

###Tabla de analisis de varianza mod = aov(rto ~ gen, datos) summary(mod)

Ejemplo 2

set.seed(123)
datos = data.frame(
  fert = gl(5,12,60, c('ctrl','D1','D2','D3','D4')),
  AF = sort(rnorm(60, 10, 0.8))
)
boxplot(AF ~ fert, datos)

mod2 = aov(AF ~ fert, datos)
summary(mod2)

##             Df Sum Sq Mean Sq F value Pr(>F)    
## fert         4 28.541   7.135   142.5 <2e-16 ***
## Residuals   55  2.754   0.050                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1