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:[email protected]
[email protected]
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
