RamalhoEg11.4: Produ<c3><a7><c3><a3>o de Gr<c3><a3>os de Milho

Description Format Source Examples

Description

Experimento conduzido na EMBRAPA Milho e Sorgo em 1997, envolvendo a avalia<c3><a7><c3><a3>o de produ<c3><a7><c3><a3>o de gr<c3><a3>os de 36 cultivares de milho, conduzido no delineamento de l<c3><a1>tice quadrado parcialmente balanceado.

Format

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

bloc

Fator de 6 n<c3><ad>veis, usado para controle local.

cult

Fator de 36 n<c3><ad>veis. Tratamento aplicado para verificar a produ<c3><a7><c3><a3>o de gr<c3><a3>os de milho.

prod

Produ<c3><a7><c3><a3>o de gr<c3><a3>os de milho, medidos em kg/parcela nas parcelas.

Source

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

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
library(lattice)

data(RamalhoEg11.4)

str(RamalhoEg11.4)

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

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

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

pet-estatistica/labestData documentation built on May 24, 2017, 3:35 a.m.