Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
Ambos lados da revisão anterior Revisão anterior | |||
pessoais:walmes:cursorufgd [2011/09/27 10:18] walmes [section 17] |
pessoais:walmes:cursorufgd [2012/04/03 22:14] (atual) walmes [section 16] |
||
---|---|---|---|
Linha 373: | Linha 373: | ||
quadro <- cbind(FV=rownames(anova(m1)), GL=anova(m1)[,"Df"], QMR) | quadro <- cbind(FV=rownames(anova(m1)), GL=anova(m1)[,"Df"], QMR) | ||
write.table(quadro, file="resumoanova.xls", quote=FALSE, row.names=FALSE, sep="\t") | write.table(quadro, file="resumoanova.xls", quote=FALSE, row.names=FALSE, sep="\t") | ||
- | </code> | ||
+ | require(agricolae) | ||
+ | |||
+ | # aplica o teste de Tukey (aqui só para efeito principal de potássio) | ||
+ | # em caso de interação o procedimento é outro | ||
+ | tukey <- lapply(ajustes, | ||
+ | function(a){ | ||
+ | qmr <- anova(a)[,"Mean Sq"] | ||
+ | dfr <- anova(a)[,"Df"] | ||
+ | da <- a$model | ||
+ | t0 <- with(da, HSD.test(r, potassio, dfr, qmr)) | ||
+ | return(t0) | ||
+ | }) | ||
+ | </code> | ||
==== Lista de afazeres ==== | ==== Lista de afazeres ==== | ||