labestData: Biblioteca de Dados para Aprendizado de Estatística

Description Format Source Examples

Com o interesse de testar a hip<c3><b3>tese de igualdade entre duas variedades de milho, foram mensuradas as vari<c3><a1>veis aleat<c3><b3>rias produtividade e altura das plantas em cada uma das variedades.

Um data.frame com 11 observa<c3><a7><c3><b5>es e 3 vari<c3><a1>veis.

varie: Fator que indica a qual variedade a observa<c3><a7><c3><a3>o pertence. Variedade A, com 6 observa<c3><a7><c3><b5>es e B com 5.
prod: Produtividade mensurada em t ha^{-1} (tonelada por hectare).
altura: Altura m<c3><a9>dia das plantas em metros.

FERREIRA (2011), Exemplo 7.4 p<c3><a1>g. 302.

data(FerreiraEg7.4)
str(FerreiraEg7.4)

aggregate(prod ~ varie, data = FerreiraEg7.4, summary)
aggregate(altura ~ varie, data = FerreiraEg7.4, summary)
by(FerreiraEg7.4[2:3], FerreiraEg7.4[1], cov)

layout(1)
with(FerreiraEg7.4, {
     par(mar = c(4, 5, 4, 5))
     plot.default(y = prod, xlab = "", ylab = "",
                  x = jitter(as.numeric(varie), factor = 0.5) - 0.1,
                  xlim = c(0.5, 2.5),
                  col = 4,
                  axes = FALSE)
     box()
     axis(side = 2, col.axis = 4)
     mtext(side = 2, text = "Produtividade em t/ha",
           line = 3, col = 4)
     par(new = TRUE, mar = c(4, 5, 4, 5))
     plot.default(y = altura, xlab = "", ylab = "",
                  x = jitter(as.numeric(varie), factor = 0.5) + 0.1,
                  xlim = c(0.5, 2.5),
                  col = 2,
                  axes = FALSE)
     axis(side = 4, col.axis = 2)
     axis(side = 1, at = 1:2,
          labels = c("Variedade A", "Variedade B"))
     mtext(side = 4, text = "Altura em m",
           line = 3, col = 2)
})