<div dir="ltr"><div>Prezados(as)</div><div><br></div><div>Comparações múltiplas (inferência simultânea) feitas através de funções paramétricas com o pacote multcomp do R são válidas para modelos lineares ou para o preditor linear dos modelos lineares generalizados. No caso de modelo de Poisson (ligacao log), modelo logit, etc., os resultados (conclusões) nos preditores lineares podem ser estendidos para a resposta média?</div>
<div><br></div><div><br></div><div>Veja o exemplo seguinte do modelo logit com 4 tratamentos e uma covariável onde fiz os intervalos simultâneos na resposta média:</div><div><br></div><div>Trat  = rep(gl(4, 4, label = c('T1','T2', 'T3', 'T4')),3)</div>
<div>y=c(1,0,1,1,0,1,0,0,1,1,0,1,1,0,0,0,0,0,0,1,1,1,0,0,0,1,1,1,1,0,0,1,0,0,0,1,1,1,0,0,1,1,1,</div><div>    1,0,0,1,0)</div><div>y_ant=c(0,0,0,1,0,1,0,0,0,0,1,0,0,0,0,1,0,1,0,1,1,0,1,0,0,1,1,1,1,0,0,0,0,0,0,1,1,1,0,0,1,</div>
<div>        1,1,1,0,0,1,0)</div><div>d1=data.frame(Trat,y_ant,y)</div><div>attach(d1)</div><div>gmod1 <- glm(y ~ Trat+y_ant, data = d1,family = binomial())</div><div>summary(gmod1)</div><div>#</div><div># Comparacao de pares de tratamentos (no preditor linear)</div>
<div>library(multcomp)</div><div>comp1=glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE))</div><div>summary(comp1, test = adjusted("Westfall"))  # Wastfall (1997)</div><div>#</div><div>
# Teste e IC's para diferencas de pares de tratamentos no preditor linear</div><div>gmod1_ci <- confint(glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE)))</div><div>summary(gmod1_ci)</div>
<div>gmod1_ci$confint</div><div>#</div><div># Comparacao de pares de tratamentos (nas prob condicionais ajustadas)</div><div>gmod1_ci <- confint(glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE)))</div>
<div>gmod1_ci$confint <- apply(gmod1_ci$confint, 1, binomial()$linkinv)</div><div>gmod1_ci$confint</div><div>#</div><div>gmod1_ci <- confint(glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE)))</div>
<div>gmod1_ci$confint <- apply(gmod1_ci$confint, 2, binomial()$linkinv)</div><div>gmod1_ci$confint</div><div><br></div><div><br></div><div>Desde já Grato,</div><div><br></div>-- <br>Gilenio Borges Fernandes<div>Professor Associado III<br>
Universidade Federal da Bahia<br>Instituto de Matemática<br>Departamento de Estatística<br>Av. Adhemar de Barros, s/n – Ondina.<br>40.170-110 - Salvador - BA, Brasil<br>Tel.: (071)3283-6280/6336  Fax:  (071)3283-6276<br>URL: <a href="http://lattes.cnpq.br/6764860618464860" target="_blank">http://lattes.cnpq.br/6764860618464860</a><br>
</div>
</div>