Prezados(as) 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? Veja o exemplo seguinte do modelo logit com 4 tratamentos e uma covariável onde fiz os intervalos simultâneos na resposta média: Trat = rep(gl(4, 4, label = c('T1','T2', 'T3', 'T4')),3) 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, 1,0,0,1,0) 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, 1,1,1,0,0,1,0) d1=data.frame(Trat,y_ant,y) attach(d1) gmod1 <- glm(y ~ Trat+y_ant, data = d1,family = binomial()) summary(gmod1) # # Comparacao de pares de tratamentos (no preditor linear) library(multcomp) comp1=glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE)) summary(comp1, test = adjusted("Westfall")) # Wastfall (1997) # # Teste e IC's para diferencas de pares de tratamentos no preditor linear gmod1_ci <- confint(glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE))) summary(gmod1_ci) gmod1_ci$confint # # Comparacao de pares de tratamentos (nas prob condicionais ajustadas) gmod1_ci <- confint(glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE))) gmod1_ci$confint <- apply(gmod1_ci$confint, 1, binomial()$linkinv) gmod1_ci$confint # gmod1_ci <- confint(glht(gmod1, linfct = mcp(Trat = "Tukey",covariate_average = TRUE))) gmod1_ci$confint <- apply(gmod1_ci$confint, 2, binomial()$linkinv) gmod1_ci$confint Desde já Grato, -- Gilenio Borges Fernandes Professor Associado III Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 URL: http://lattes.cnpq.br/6764860618464860