Na realidade tem como você enganar a glht() para usar a cld() mesmo sem passar para o argumento linfct= algo produzido pela mcp(). Foi preciso abrir a função para encontrar como fazer. Essa é a maravilha do software livre código aberto! Você pode visitar essa análise e ver como fiz, basicamente criei elementos não nulos na lista retornada pela glht().

http://www.leg.ufpr.br/~walmes/analises/GBTognon/ivg_tridentata.html

A parte que destaco da análise é essa

Xc <- apc(Xm, lev=lev)

g0 <- summary(glht(m2, linfct=Xc), test=adjusted(type="fdr"))

source("/home/walmes/Dropbox/Public/mycld.R")

## g0

g0$focus <- "Tratamentos"

g0$Type <- "Tukey"

grid$cld <- my.cld(g0)$mcletters$Letters

Se não me falhe a memória, essa função my.cld eu refiz e renomeei para cld2. Está disponível em https://github.com/walmes/wzRfun/blob/master/R/cld2.R. Alguma adaptação será necessária.

Esse outro script http://www.leg.ufpr.br/~walmes/cursoR/geneticaEsalq/script20.html usa a cld2().

À disposição.

Walmes.