Re: [R-br] Matrix de contrastes em fatorial DBC com a função glht!

Prezado Walmes! Como sempre, muito obrigado! A função apcMatrix é muitíssimo útil. Uma sugestão! Embora eu saiba que tu não gostes de letrinhas, talvez seria interessante criar alguma função semelhante a função "cld", já que a mesma só trabalha caso seja indicado no argumento "linfct" a função "mcp" com o nome do tratamento de interesse. Em casos mais gerais como foi neste exemplo, a função "cld" não trabalha. Abraços!

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.

Poderia fornecer o arquivo mycld.R para rodar? 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. --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

source("https://dl.dropboxusercontent.com/u/48140237/mycld.R") À disposição. Walmes.
participantes (3)
-
Ivan Bezerra Allaman
-
Mauro Sznelwar
-
Walmes Zeviani