labestData: Biblioteca de Dados para Aprendizado de Estatística

Description Format Source Examples

Experimento da avalia<c3><a7><c3><a3>o da produ<c3><a7><c3><a3>o de 40 cultivares de sorgo conduzido no delineamento alfa-l<c3><a1>tice.

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

bloc: Fator de 8 n<c3><ad>veis, usado para controle local. Os blocos s<c3><a3>o incompletos, ou seja, o tamanho do bloco <c3><a9> menor que o n<c3><ba>mero de n<c3><ad>veis de cultivares.
cult: Fator de 40 n<c3><ad>veis que s<c3><a3>o as cultivares de sorgo.
prod: Produ<c3><a7><c3><a3>o de gr<c3><a3>os de sorgo, medidos em kg/parcela.

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

library(lattice)

data(RamalhoEg11.10)

str(RamalhoEg11.10)

addmargins(xtabs(~cult + bloc, data = RamalhoEg11.10))

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

k <- nlevels(RamalhoEg11.10$cult)
a <- seq(0, 2 * pi, length.out = k + 1)[-(k + 1)]
par(mfrow = c(2, 4))
col <- 1
for (b in levels(RamalhoEg11.10$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(RamalhoEg11.10, 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(RamalhoEg11.10$cult)[i])
    col <- col + 1
}