Walmes,

           Muito mais eficiente que fazer na mão, problema resolvido e abordagem adotada.

Obrigado,

Alexandre


Em 09/05/2014 11:14, walmes . escreveu:

##-----------------------------------------------------------------------------
## No caso de desdobrar interação dupla.

M <- LSmatrix(g0, effect=c("a","b"))
grid <- attr(M, "grid")

## Matriz de contrastes entre níveis de `focus` dentro dos níveis de
## `split`.
split <- c("b")
focus <- "a"
spl <- interaction(grid[,split])
i <- 1:nrow(grid)
l <- split(i, f=spl)
contr <- lapply(l,
                function(row){
                    ## Matriz de contrastes par a par.
                    a <- apc(M[row,], lev=levels(d[,focus]))
                    ## Prefixo no nome das linhas.
                    rownames(a) <- paste(spl[row[1]],
                                         rownames(a), sep="/")
                    return(a)
                })
contr <- do.call(rbind, contr)
contr

## Constrastes.
summary(glht(g0, linfct=contr),
        test=adjusted(type="fdr"))

##-----------------------------------------------------------------------------
## Representando em um gráfico média com IC.

grid <- attr(M, "grid")
grid <- sapply(grid, function(x) if(is.character(x)) factor(x) else x)
means <- confint(glht(g0, linfct=M), calpha=univariate_calpha())
grid <- cbind(data.frame(grid), means$confint)
str(grid)

segplot(a~lwr+upr|b, data=grid, draw=FALSE, centers=Estimate)

##-----------------------------------------------------------------------------

-- 
======================================================================
Alexandre dos Santos
Proteção Florestal 
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial 
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr@yahoo.com.br 
        alexandre.santos@cas.ifmt.edu.br 
Lattes: http://lattes.cnpq.br/1360403201088680                   
======================================================================