labestData: Biblioteca de Dados para Aprendizado de Estatística

Description Format Source Examples

Ensaio realizado pelo Departamento de Produ<c3><a7><c3><a3>o Animal da Secretaria de Agricultura paulista, que teve por objetivo estudar a influ<c3><aa>ncia da case<c3><ad>na iodada na nutri<c3><a7><c3><a3>o de vacas leiteiras. Consideram-se 4 tratamentos (doses de caseina iodada) e fizeram-se 3 repeti<c3><a7><c3><b5>es, em blocos casualizados com cada bloco contendo 4 parcelas, cada uma constitu<c3><ad>da por uma vaca. As vacas eram todas mesti<c3><a7>as Flamengo \times Caracu e todas paridas em maio, junho ou julho de 1959. Os blocos foram organizados tendo em vista a produ<c3><a7><c3><a3>o de leite das vacas em um per<c3><ad>odo pr<c3><a9>-experimental de duas semanas, em que todas recebiam uma mesma ra<c3><a7><c3><a3>o. No primeiro bloco ficaram as vacas de maior produ<c3><a7><c3><a3>o, no segundo as 4 seguintes e no terceiro as de menor produ<c3><a7><c3><a3>o.

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

bloc: Fator de 3 n<c3><ad>veis qualitativos, usado para controle local.
caseina: Dosagem de case<c3><ad>na iodada, em gramas.
prod: Produ<c3><a7><c3><a3>o de leite das vacas, em kg, no per<c3><ad>odo experimental de 8 semanas (de 20 de Outubro a 15 Dezembro de 1959).

A observa<c3><a7><c3><a3>o da testemunha (0 de caseina) no bloco 3 foi perdida. O valor correspondente foi estimado por m<c3><a9>todos de imputa<c3><a7><c3><a3>o de parcela perdida discutido na se<c3><a7><c3><a3>o 5.4 de PIMENTEL-GOMES (2009).

PIMENTEL-GOMES (2009), p<c3><a1>gina 269.

library(lattice)

data(PimentelPg269)
str(PimentelPg269)

xyplot(prod ~ caseina, data = PimentelPg269,
       groups = bloc, type = "o",
       auto.key = list(title = "Bloco", cex.title = 1.1,
                       columns = 3),
       ylab = "Produ<c3><a7><c3><a3>o de Leite (kg)",
       xlab = "Dose de Case<c3><ad>na Iodada (gramas)")

# Parcela perdida que foi estimada.
subset(PimentelPg269, caseina == 0 & bloc == "3")

i <- with(PimentelPg269, caseina == 0 & bloc == "3")
PimentelPg269$prod[i] <- NA