[R-br] Contrastes - MLG

[Cesar Rabak]

André,

A função glht está usando a ordem "alfabética" dos seus fatores, como a
função não deve ter meios de mudar isso (você já *olhou* no FM, não
*olhou*?) 😉,
a única forma seria fazer uma cópia do objeto retornado 'p' neste teu
exemplo, e inverter os sinais de *todos* os contrastes e colocar os
*labels* manualmente
no gráfico.

HTH
--
Cesar Rabak


On Wed, Aug 7, 2019 at 4:50 PM Andre Oliveira por (R-br) <
r-br em listas.c3sl.ufpr.br> wrote:

> Boa tarde!
>
> Tem como inverter a ordem do contraste que vem por padrão na função glht
> do pacote multcomp? Digo inverter para bater com o padrão do resultado do
> summary(fit) onde coloca t1-t2 e não t2 - t1 como fica no objeto p.
>
> dados <- read.table("http://www.leg.ufpr.br/~paulojus/dados/exemplo01.txt
> <https://l.facebook.com/l.php?u=http%3A%2F%2Fwww.leg.ufpr.br%2F%7Epaulojus%2Fdados%2Fexemplo01.txt%3Ffbclid%3DIwAR2krDEvOI3rtXw4xASIrG0AxtCOmoRmbY0gt0SGuBk-7g3PhQcb3fS8o9o&h=AT3yPMcJkTSq7ku8iquqWH8fXOHRd76dpt80tMGemb1IX4c-yDBvcK2kBDS9NehkvTxIwEZkPdQWWtusjgFnfpytgzMKEyTEuIfdENjzZl7RO8VxxlRP3fhEUIwedOiNIZlwR5o_3H0>",
> head = T)
>
> attach(dados)
>
> boxplot(resp~trat, col=8)
>
> hist(resp, col=8)
>
> fit=glm(resp~trat, family = Gamma)
>
> summary(fit)
>
> require(multcomp)
>
> p <- glht(fit, linfct = mcp(trat = "Tukey"))
> plot(p, main= "", xlab = "")
>
> obg
>
> att.
> André
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
> Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça
> código mínimo reproduzível.
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20190807/54b9507a/attachment.html>