
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", 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é

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@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@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.
participantes (2)
-
Andre Oliveira
-
Cesar Rabak