[R-br] Comparações múltiplas em modelos lineares generalizados

[Gilenio Borges Fernandes]

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
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130702/2f787582/attachment.html>