RamalhoEx8.1: Produ<c3><a7><c3><a3>o de Gr<c3><a3>os de Arroz em...

In pet-estatistica/labestData: Biblioteca de Dados para Aprendizado de Estatística

Description

Experimentos para avaliar a produ<c3><a7><c3><a3>o de gr<c3><a3>os de cultivares de arroz, conduzidos em tr<c3><aa>s locais do Estado de Minas Gerais (Lambari, Lavras e Felixl<c3><a2>ndia) em delineamento inteiramente casualizado.

Format

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

cult

Fator de n<c3><ad>veis nominais que representa as cultivares de arroz.

local

Fator de 3 n<c3><ad>veis nominais, munic<c3><ad>pios de Minas Gerais onde foram realizados os experimentos.

rept

N<c3><ba>meros inteiros que identificam as repeti<c3><a7><c3><b5>es em cada local.

prod

Produ<c3><a7><c3><a3>o de gr<c3><a3>os de arroz, em kg ha^{-1}.

Source

RAMALHO et al. (2005), Exerc<c3><ad>cio 8.1, p<c3><a1>g. 132.

Examples

data(RamalhoEx8.1)

str(RamalhoEx8.1)

with(RamalhoEx8.1,
     tapply(prod, list(cult = cult, local = local),
            FUN = mean))

aggregate(prod ~ cult + local, data = RamalhoEx8.1,
          FUN = function(x) { c(mean = mean(x), var = var(x)) })

library(lattice)

xyplot(prod ~ cult | local, data = RamalhoEx8.1,
       as.table = TRUE, layout = c(NA, 1), type = c("p", "a"),
       xlab = "Cultivar de arroz",
       ylab = expression("Produ<c3><a7><c3><a3>o"~(kg~ha^{-1})))

# Ordenar cultivares pela m<c3><a9>dia dentro de cada local.
a <- by(data = RamalhoEx8.1, INDICES = RamalhoEx8.1$local,
        FUN = function(d) {
            with(d, reorder(interaction(cult, local, drop = TRUE),
                            prod))
        })
a <- levels(unlist(a)); a

RamalhoEx8.1 <- transform(RamalhoEx8.1,
                          cult.loc = factor(interaction(cult, local),
                                            levels = a))

xyplot(prod ~ cult.loc | local, data = RamalhoEx8.1,
       as.table = TRUE, layout = c(NA, 1), type = c("p", "a"),
       scales = list(x = "free"),
       xlab = "Cultivar de arroz (ordenadas por local)",
       ylab = expression("Produ<c3><a7><c3><a3>o"~(kg~ha^{-1})),
       xscale.components = function(...) {
           ans <- xscale.components.default(...)
           ans$bottom$labels$labels <-
               gsub(x = ans$bottom$labels$labels,
                    pattern = "\\..*$",
                    replacement = "")
           return(ans)
       })

pet-estatistica/labestData documentation built on May 25, 2019, 12:47 a.m.