[R-br] Contrastes - MLG

classic Classic list List threaded Threaded
2 messages Options
Reply | Threaded
Open this post in threaded view
|

[R-br] Contrastes - MLG

R-br mailing list

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é 

_______________________________________________
R-br mailing list
[hidden email]
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.
Reply | Threaded
Open this post in threaded view
|

Re: [R-br] Contrastes - MLG

R-br mailing list
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) <[hidden email]> 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", 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
[hidden email]
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.

_______________________________________________
R-br mailing list
[hidden email]
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.