labestData: Biblioteca de Dados para Aprendizado de Estatística

Description Format Source Examples

Experimento conduzido na EMBRAPA Milho e Sorgo envolvendo a avalia<c3><a7><c3><a3>o de 16 cultivares de sorgo conduzido no delineamento de l<c3><a1>tice quadrado balanceado.

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

bloc: Fator de 4 n<c3><ad>veis, usado para controle local.
cult: Fator de 16 n<c3><ad>veis. Tratamento aplicado para verificar a produ<c3><a7><c3><a3>o de sorgo.
prod: Produ<c3><a7><c3><a3>o de sorgo, medidos em kg/parcela nas parcelas.

RAMALHO et al. (2005), p<c3><a1>g. 165.

library(lattice)

data(RamalhoTb11.1)

str(RamalhoTb11.1)

addmargins(xtabs(~cult + bloc, data = RamalhoTb11.1))

xyplot(prod ~ cult, data = RamalhoTb11.1,
       xlab = "Cultivares de sorgo",
       ylab = expression("Produ<c3><a7><c3><a3>o de gr<c3><a3>os"~(kg~parcela^{-1})))

k <- nlevels(RamalhoTb11.1$cult)
a <- seq(0, 2 * pi, length.out = k + 1)[-(k + 1)]
par(mfrow = c(2, 2))
col <- 1
for (b in levels(RamalhoTb11.1$bloc)) {
    plot(sin(a), cos(a), asp = 1,
         xlim = c(-1.1, 1.1),
         ylim = c(-1.1, 1.1),
         axes = FALSE, xlab = NA, ylab = NA)
    mtext(paste("Bloco", b))
    i <- unique(as.integer(subset(RamalhoTb11.1, bloc == b)$cult))
    cb <- combn(x = i, m = 2)
    segments(x0 = sin(a[cb[1, ]]), y0 = cos(a[cb[1, ]]),
             x1 = sin(a[cb[2, ]]), y1 = cos(a[cb[2, ]]),
             col = col)
    text(x = 1.08 * sin(a[i]), y = 1.08 * cos(a[i]),
         labels = levels(RamalhoTb11.1$cult)[i])
    col <- col + 1
}