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 ======================================================================