par(mfrow = c(1,1))

Dia0

dia0<-aov(proyecto$ dia0 ~ proyecto $Planta)

shapiro.test(dia0$residuals)

bartlett.test(proyecto$ dia0, proyecto $Planta)

TukeyHSD(dia0)

summary(dia0)

boxplot(proyecto$ dia0 ~ proyecto $Planta, xlab = “Tratamiento”, ylab = “Longitud de Tallo (cm)”)

Dia1

dia1<-aov(proyecto$ dia1 ~ proyecto $Planta)

shapiro.test(dia1$residuals)

fligner.test(proyecto$ dia1, proyecto $Planta)

kruskal.test(proyecto$ dia1 ~ proyecto $Planta)

pairwise.t.test(proyecto$ dia1, proyecto$Planta)

boxplot(proyecto$ dia1 ~ proyecto $Planta,xlab = “Tratamiento”, ylab = “Longitud de Tallo (cm)”)

Dia2 dia2<-aov(proyecto$ dia2 ~ proyecto $Planta)

shapiro.test(dia2$residuals)

bartlett.test(proyecto$ dia2, proyecto $Planta)

summary(dia2)

boxplot(proyecto$ dia2 ~ proyecto $Planta,xlab = “Tratamiento”, ylab = “Longitud de Tallo (cm)”)

Dia3 dia3<-aov(proyecto$ dia3 ~ proyecto $Planta)

shapiro.test(dia3$residuals)

fligner.test(proyecto$ dia3, proyecto $Planta)

kruskal.test(proyecto$ dia3 ~ proyecto $Planta)

boxplot(proyecto$ dia3 ~ proyecto $Planta, xlab = “Tratamiento”, ylab = “Longitud de Tallo (cm)”)

Dia4 dia4<-aov(proyecto$ dia4 ~ proyecto $Planta)

shapiro.test(dia4$residuals)

fligner.test(proyecto$ dia4, proyecto $Planta)

kruskal.test(proyecto$ dia4 ~ proyecto $Planta)

boxplot(proyecto$ dia4 ~ proyecto $Planta,xlab = “Tratamiento”, ylab = “Longitud de Tallo (cm)”)

Peso Raiz pesor<-aov(proyecto$ pesoraiz ~ proyecto $Planta)

shapiro.test(pesor$residuals)

bartlett.test(proyecto$ pesoraiz, proyecto $Planta)

summary(pesor)

TukeyHSD(pesor)

boxplot(proyecto$ pesoraiz ~ proyecto $Planta,xlab = “Tratamiento”, ylab = “Peso seco de Raíz (g)”)

Peso Tallo pesot<-aov(proyecto$ pesotallo ~ proyecto $Planta)

shapiro.test(pesot$residuals)

bartlett.test(proyecto$ pesotallo, proyecto $Planta)

summary(pesot)

boxplot(proyecto$ pesotallo ~ proyecto $Planta,,xlab = “Tratamiento”, ylab = “Peso seco de Tallo (g)”)

Hojas hojas<-aov(proyecto$ hojas ~ proyecto $Planta)

shapiro.test(hojas$residuals)

fligner.test(proyecto$ hojas, proyecto $Planta)

kruskal.test(proyecto$ hojas ~ proyecto $Planta)

boxplot(proyecto$ hojas~ proyecto $Planta,xlab = “Tratamiento”, ylab = “Número de Hojas”)

Clorosis por Chi cloro<-matrix(c(10,10,17,3,13,7,10,10),nrow=2)

cloro

colnames(cloro)= c(“control”, “tratamientoA”, “tratamientoB”, “tratamientoC”)

rownames(cloro)=c(“Ausencia”, “Presencia”)

cloro

chisq.test(cloro)

prop.table(cloro)

barplot(cloro, legend = rownames(cloro), beside=F, axis.lty = 1,xlab = “Clorosis”, ylab = “Cantidad de plantas”)

Necrosis por Chi necro<-matrix(c(19,1,16,4,17,3,13,7),nrow=2)

necro

colnames(necro)= c(“control”, “tratamientoA”, “tratamientoB”, “tratamientoC”)

rownames(necro) =c(“Ausencia”, “Presencia”)

necro

chisq.test(necro)

prop.table(necro)

fisher.test(necro)

barplot(necro, legend = rownames(necro), beside=F, axis.lty = 1,xlab = “Necrosis”, ylab = “Cantidad de plantas”)