labestData: Biblioteca de Dados para Aprendizado de Estatística

Description Format Source Examples

Experimento que estudou a produ<c3><a7><c3><a3>o de milho como fun<c3><a7><c3><a3>o de fontes e doses de aduba<c3><a7><c3><a3>o nitrogenada. O experimento foi instalado em delineamento de blocos casualizados com arranjo de tratamento em parcela subdividida, sendo as doses casualizadas dentro dos n<c3><ad>veis de adubo.

Um data.frame com 48 observa<c3><a7><c3><b5>es e 4 vari<c3><a1>veis, em que

adubo: Fator de n<c3><ad>veis nominais que s<c3><a3>o os adubos fontes de nitrog<c3><aa>nio.
dose: Fator de n<c3><ad>veis m<c3><a9>tricos que s<c3><a3>o as doses aplicadas, em kg de nitrog<c3><aa>nio por hectare.
bloco: Fator de n<c3><ad>veis nominais que representa os blocos do experimento.
prod: Produ<c3><a7><c3><a3>o de milho, em kg ha ^{-1}

BANZATTO; KRONKA (2013), Quadro 6.2.5, p<c3><a1>g. 144.

library(lattice)

data(BanzattoQd6.2.5)

str(BanzattoQd6.2.5)

ftable(with(BanzattoQd6.2.5,
            tapply(prod,
                   list(adubo = adubo,
                        dose = dose,
                        bloco = bloco),
                   FUN = mean)))

xyplot(prod ~ dose, groups = adubo, data = BanzattoQd6.2.5,
       as.table = TRUE, jitter.x = TRUE, type = c("p", "a"),
       auto.key = list(title = "Adubos fontes de nitrog<c3><aa>nio",
                       cex.title = 1, columns = 2),
       xlab = expression(
           "Quantidade aplicada de nitrog<c3><aa>nio"~(kg~ha^{-1})),
       ylab = expression("Produ<c3><a7><c3><a3>o de milho"~(kg~ha^{-1})))

xyplot(prod ~ dose | adubo, data = BanzattoQd6.2.5,
       as.table = TRUE, jitter.x = TRUE,
       xlab = expression(
           "Quantidade aplicada de nitrog<c3><aa>nio"~(kg~ha^{-1})),
       ylab = expression("Produ<c3><a7><c3><a3>o de milho"~(kg~ha^{-1})),
       panel = function(x, y, ...){
           panel.xyplot(x, y, ...)
           m0 <- lm(y ~ poly(x, degree = 2))
           panel.curve(predict(m0, newdata = list(x = x)))
       })