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 ==== | ||