
7 Ago
2019
7 Ago
'19
19:47
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é