[R-br] Duplo fatorial desbalanceado - Soluções simples de análise

[Jair]

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