Cómo determinar si hay diferencias significativas en más de dos tratamientos en RStudio, lenguaje R?

 Para pruebas parametricas: 

renmaizaov <- aov(totparcelas$renmaiz ~ factorparcelas)

summary(renmaizaov) #valor p < 0.05

shapiro.test(renmaizaov$residuals) #valor p > 0.05

TukeyHSD(renmaizaov)

Para pruebas no parametricas:

shapiro.test(Pextarvensesaov$residuals) #valor p < 0.05

Pextarvensesaov <- kruskal.test(totparcelas$Pextarvenses ~ factorparcelas) #valor p < 0.05

dunnTest(totparcelas$Pextarvenses, factorparcelas, method = "bonferroni")