<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Você pode obter as estimativas de cada casela de diversas formas: 1) removendo intercepto da fórmula, 2) mudando o tipo de constraste, 3) calculando as funções lineares por fora com matrizes apropriadas e 4) fazendo o item 3 com a popMeans() e/ou popMatrix(), que considero a opção mais simples. Veja o CMR.<br>
<br><span style="font-family:courier new,monospace">##-----------------------------------------------------------------------------<br><br>## dados artificiais<br>da <- expand.grid(A=gl(4,4), B=gl(3,1))<br>da$y <- rpois(nrow(da), lambda=10)<br>
<br>## constrates usados por padrão<br>options()$contrasts<br><br>## removendo o intercepto, primeiro fator com estimativas por nível<br>m0 <- glm(y~-1+A+B, data=da, family=poisson)<br>summary(m0)<br>m0 <- glm(y~-1+B+A, data=da, family=poisson)<br>
summary(m0)<br><br>## mudando a opção de contraste localmente <br>m0 <- glm(y~A+B, data=da, family=poisson,<br>          contrasts=list(A=contr.helmert, B=contr.SAS))<br>summary(m0)<br><br>M <- unique(model.matrix(m0))<br>
M<br><br>## comportamento padrão<br>m0 <- glm(y~A+B, data=da, family=poisson)<br>summary(m0)<br><br>M <- unique(model.matrix(m0))<br>M<br><br>##-----------------------------------------------------------------------------<br>
<br>## estimativas para cada combinação A*B<br>M%*%coef(m0)<br><br>## para soma direta de matrizes<br>source("<a href="http://fisher.osu.edu/~schroeder_9/AMIS900/blockdiag.R">http://fisher.osu.edu/~schroeder_9/AMIS900/blockdiag.R</a>")<br>
<br>mA <- contrasts(da$A)<br>mB <- contrasts(da$B)<br>M <- list(mA, mB)<br>M <- cbind(1, blockdiag(M))<br>M <- unique(M)<br>M<br><br>M%*%coef(m0)<br>coef(m0)<br><br>require(doBy)<br><br>## efeito marginal de A (ou seja, na média dos níveis de B)<br>
popMeans(m0, effect="A")<br><br>## efeito marginal de B (ou seja, na média dos níveis de A)<br>popMeans(m0, effect="B")<br><br>## efeito de A condicional à B=1<br>popMeans(m0, effect="A", at=list(B="1"))<br>
<br>## aqui é para ver as matrizes envolvidas<br>popMatrix(m0, effect="A")<br>popMatrix(m0, effect="A", at=list(B="1"))<br><br>##-----------------------------------------------------------------------------<br>
</span><br>À disposição.<br>Walmes.<br></div></div>