[R-br] Desdobramento de contrastes GLM via glht()

[ASANTOS]

Boa noite pessoal,

     Gostaria de saber se existe alguma forma de se atribuírem índices 
ao teste realizado com a função  glht() do pacote multcomp, exemplificado:

require(multcomp)
y1<-c(rpois(100, lambda = 5),rpois(100, lambda = 10),rpois(100, lambda = 
5),rpois(100, lambda = 7))##Geração da variável resposta
trat<-c(rep(1,100),rep(2,100),rep(3,100),rep(4,100))##Geração dos 
tratamentos
repeticao<-rep(1:100,4)##Imaginado um DIC colocando as repetições
trat<-as.factor(trat)
m.1<-glm(y1~trat+repeticao,family="poisson")##Ajustando GLM Poisson
anova.m.1<-anova(m.1,test="Chi")
anova.m.1##Verificando a significancia dos tratamentos e se não há 
efeito das repetições
summary(glht(m.1,linfct=mcp(trat="Tukey")))##Desdobrando os tratamentos

Como a função utiliza o Tukey acho bem complicado como no caso simulado, 
decidir entre tratamentos que possuem índice duplo, por exemplo, alguém 
conhece alguma forma de calcular a decisão de distribuição dos índices,
Obrigado,
Alexandre