Velocidadagitacion.docx

velocidadagitacion<-c(5,5,5,5,10,10,10,10,15,15,15,15,20,20,20,20) bloquehorno<-c(1,2,3,4,1,2,3,4,1,2,3,4,1,2,3,4) textura<-c(8,4,5,6,14,5,6,9,14,6,9,2,17,9,3,6)

#Ho: textura i = textura j #Hi: textura i =/ textura j

velocidadagitacion<-factor(velocidadagitacion) bloquehorno<-factor(bloquehorno) #tabla anova modelo<-lm(textura~bloquehorno+velocidadagitacion) ANOVA<-aov(modelo) summary(ANOVA) qf(p=0.05,df1 = 3,df2 = 9,lower.tail = F) pf(q=0.853,df1 = 3,df2 = 9,lower.tail = F) library(agricolae) grupos.velocidadagitacion<-LSD.test(y=ANOVA,trt = "velocidadagitacion",group = T,console = T) bar.group(x=grupos.velocidadagitacion$groups) # estadisticas de LSD grupos.velocidadagitacion$statistics # prueba de rangos multiples de duncan grupos.velocidadagitacion.duncan<-duncan.test(y=ANOVA,trt = "velocidadagitacion",group = T,console = T) bar.group(x=grupos.velocidadagitacion.duncan$groups)

# validacion

qqnorm(rstandard(modelo)) qqline(rstandard(modelo)) shapiro.test(rstandard(modelo)) bartlett.test(textura~velocidadagitacion)