PaulaEx4.6.15: N<c3><ba>mero de Infec<c3><a7><c3><b5>es de Ouvido em...

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

Description

Dados referentes a um estudo realizado em 1990 com recrutas americanos em que a vari<c3><a1>vel de interesse era o n<c3><ba>mero de infec<c3><a7><c3><b5>es de ouvido. Al<c3><a9>m disso, foram coletadas as seguintes informa<c3><a7><c3><b5>es sobre os recrutas: h<c3><a1>bito de nadar, local em que costuma nadar, idade e sexo.

Format

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

habito

Fator com dois n<c3><ad>veis que indica o h<c3><a1>bito em nadar do recruta (ocasional ou frequente).

local

Fator com dois n<c3><ad>veis que indica o local onde o recruta costuma nadar, (praia ou piscina).

idade

Idade do recruta, em anos categorizados em tr<c3><aa>s n<c3><ad>veis (15-19, 20-24 e 25-29).

sexo

Sexo (F para feminino e M para masculino).

ninfec

N<c3><ba>mero de infec<c3><a7><c3><b5>es de ouvido diagnosticadas pelo pr<c3><b3>prio recruta.

Source

PAULA (2004), Exerc<c3><ad>cio 4.6.15, p<c3><a1>g. 346.

References

Hand, D. J, Daly, F., Lunn, A. D., McConway, K. J., Ostrowski, E. (1994). A Handbook of Small Data Sets. Chapman and Hall, London.

Examples

data(PaulaEx4.6.15)
str(PaulaEx4.6.15)

xt <- xtabs(ninfec ~ ., data = PaulaEx4.6.15)
ftable(prop.table(xt))
plot(xt, xlab = "H<c3><a1>bito", main = "")

mv <- aggregate(ninfec ~ .,
                data = PaulaEx4.6.15,
                FUN = function(x) {
                    c(mu = mean(x), var = var(x))
                })

library(lattice)

# Rela<c3><a7><c3><a3>o M<c3><a9>dia Vari<c3><a2>ncia
xyplot(ninfec[, "var"] ~ ninfec[, "mu"],
        data = mv,
        ylab = "Vari<c3><a2>ncia amostral",
        xlab = "M<c3><a9>dia amostral",
        panel = function(x, y) {
            panel.xyplot(x, y, type = c("p", "r"), grid = TRUE)
            panel.abline(a = 0, b = 1, lty = 2)
        })

library(latticeExtra)

useOuterStrips(
    xyplot(ninfec ~ idade | habito + local,
           groups = sexo,
           data = PaulaEx4.6.15,
           jitter.x = TRUE, jitter.y = TRUE,
           type = c("p", "g", "a"),
           ylab = "N<c3><ba>mero de infec<c3><a7><c3><b5>es",
           auto.key = list(space = "right", cex.title = 1,
                           title = "Sexo"))
    )

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