Opcional:

Nesta seção você poderá realizar a ANOVA usando operações aritméticas no R. Utilize o arquivo de dados sobre linhagens como exemplo.

n <- length(resp)
n
nt <- length(levels(trat))
nt

correcao <- ((sum(resp))^2)/n
correcao

gltot <- n - 1
gltra <- nt - 1
glres <- n - nt

sqtot <- sum(resp^2) - correcao
sqtot

trtot <- tapply(resp, trat, sum)
sqtra <- (sum((trtot)^2))/6 - correcao
sqtra

sqres <- sqtot - sqtra
sqres

qmtra <- sqtra/gltra
qmres <- sqres/glres

fval <- qmtra/qmres
fval
pval <- pf(fval, gltra, glres, lower.tail = F)
pval

## Definindo o quadro da ANOVA

qav <- matrix(NA, nr=3, nc=5)
dimnames(qav) <- list(c("Tratamentos", "Residuo", "Total"),
c("GL","SQ","QM","valorF","valorP"))
qav

## Montando o Quadro da ANOVA

qav[1,] <- c(gltra, glres, qmtra, fval, pval)
qav[2,1:3] <- c(sqtra, sqres, qmres)
qav[3,1:2] <- c(gltot, sqtot)
qav