[R-br] Contraste com glm com distribuição de poisson em um fatorial 2x2 [RESOLVIDO]

ASANTOS alexandresantosbr em yahoo.com.br
Sexta Maio 9 14:19:35 BRT 2014


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 em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140509/0114f9d5/attachment.html>


Mais detalhes sobre a lista de discussão R-br