Boa noite pessoal,
Gostaria de saber se existe alguma forma de se atribuírem índices
ao teste realizado com a função glht() do pacote multcomp, exemplificado:
require(multcomp)
y1<-c(rpois(100, lambda = 5),rpois(100, lambda = 10),rpois(100, lambda =
5),rpois(100, lambda = 7))##Geração da variável resposta
trat<-c(rep(1,100),rep(2,100),rep(3,100),rep(4,100))##Geração dos
tratamentos
repeticao<-rep(1:100,4)##Imaginado um DIC colocando as repetições
trat<-as.factor(trat)
m.1<-glm(y1~trat+repeticao,family="poisson")##Ajustando GLM Poisson
anova.m.1<-anova(m.1,test="Chi")
anova.m.1##Verificando a significancia dos tratamentos e se não há
efeito das repetições
summary(glht(m.1,linfct=mcp(trat="Tukey")))##Desdobrando os tratamentos
Como a função utiliza o Tukey acho bem complicado como no caso simulado,
decidir entre tratamentos que possuem índice duplo, por exemplo, alguém
conhece alguma forma de calcular a decisão de distribuição dos índices,
Obrigado,
Alexandre
_______________________________________________
R-br mailing list
[email protected]
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.
