Diferenças

Aqui você vê as diferenças entre duas revisões dessa página.

--- pessoais:walmes:cursorufgd [2011/09/27 10:18] – [section 17] walmes
+++ pessoais:walmes:cursorufgd [2012/04/03 22:14] (atual) – [section 16] walmes
@@ Linha 373: / Linha 373: @@
 quadro <- cbind(FV=rownames(anova(m1)), GL=anova(m1)[,"Df"], QMR)
 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 ====