Pessoal, gostaria de sugestões de soluções simples para análise do experimento a seguir. #------------------------------------------------------------------ #Experimento fatorial com dois fatores em DIC #Fator 1: Solo (4 tipos de solo) #Fator 2: Profundidade (3 profundidades) #Variável resposta: Teor de cromo #Desbalanceado: não foi possível detectar a variável resposta #para alguns solos em alguma profundidade. #------------------------------------------------------------------ rm(list = ls()) solo <- rep(c("S1","S2","S3","S4"), each=15) #solos prof <- rep(rep(c("P1","P2","P3"), each=5),4) #profundidades rep <- rep(1:5, 12) #repetição cromo <- c(30.40, 4.28, 9.27, 14.46, 10.54, 34.52, 6.36, 16.89, 18.92, 12.92, 44.63, 18.53, 23.42, 23.37, 18.43, 57.54, 132.40, 129.17, 135.37, 157.58, 67.89, 154.09, 162.63, 155.19, 171.08, 70.82, 122.19, 164.69, 179.06, 160.65, 121.25, 37.54, 30.26, 46.60, NA, 81.91, 26.35, 34.91, 46.03, 13.42, NA, NA, NA, NA, NA, 1.80, 4.33, 8.57, 3.75, 2.31, 4.12, 10.72, 24.70, 6.50, 3.55, 57.91, 1.58, NA, 97.41, 56.73) dados <- data.frame(solo, prof, rep, cromo) str(dados) #dados<-na.omit(dados) #str(dados) require(lattice) xyplot(cromo~prof|solo, data = dados) xyplot(cromo~solo|prof, data = dados) library(car) #Type III SS ANOVA names(dados) m1 <- Anova(lm(cromo ~ solo * prof, dados), type = "3") m1 require(easyanova) dado <- dados[,c(1:2,4)] m2 <- ea2(data=dado, design = 1, alpha = 0.05, list = FALSE, p.adjust=1, plot=2) m2 #--------------------------------------------------------------------- Obrigado desde já, Jair