Pessoal Preciso fazer um teste de Dunnet para verificar se existe diferença entre os tratamentos (níveis 0,002, 0,004 e 0,006) e a testemunha (dose 0,000). Esse experimento é um aditivo colocado em um alimento, mas é preciso que os consumidores (julgadores) não percebam a diferença. Ou seja, não seja diferente da testemunha. Realize a análise, mas o tete de Dunnet sai comparando cada julgador, gostaria que a saída, me desse, cada tratamento em relação a testemunha que é a dose zero do banco em anexo. Usei o pacote multicomp para fazer a comparação. Desde já agradeço dados3 <- read.table("prova3.txt", h=T, dec=",") attach(dados3) names(dados3) levels(trat) trat <- factor(trat) prova3.av <- aov(resp~julga+trat) anova(prova3.av) require(multcomp) summary(glht(prova3.av)) João Peterson Pereira Gardin Doutor em Fisiologia Vegetal, Eng. Agr., Pesquisador Fones: (49) 3566-0054 e 9933-4290 MSN: petersongardin@yahoo.com.br - SKYPE: petersongardin Endereço: Rua João Zardo, 1660, Bairro Campo Experimental, Videira - SC - CEP 89560000
Se trocar a ordem dos termos na formula do modelo? Não tem como especificar níveis de qual fator você quer comparar na glht? Quando são contrastes de Tukey eu uso assim glht(m0, linfct=mcp("fatorescolhido"=Tukey)) # algo nesse sentido deve haver algo semelhante. Não tenho como ver a documentação agora. Além do mais, você não precisa obrigatóriamente do Dunnet, você pode checar essas hipóteses por contraste ortogonal, e isso sai no summary() do modelo (todos os níveis vs o de referência) se o nível de referência for a dose zero e usar o contraste tipo contr.treatment (que é o padrão). Peterson, se você passar o dput() dos dados ao inves dele em anexo seu CMR fica mais convidativo, são poucos dados, e para esse caso, pode usar dados simulados. À disposição. Walmes. -- ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (2)
-
Peterson Gardin -
Walmes Zeviani