PaulaEg2.8.1: Compara<c3><a7><c3><a3>o de Tipos de Snack

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

Description

Dados de um experimento desenvolvido pelo Departamento de Nutri<c3><a7><c3><a3>o da Faculdade de Sa<c3><ba>de P<c3><ba>blica da USP em que 5 tipos diferentes de um novo snack, com baixo teor de gordura e de <c3><a1>cidos graxos, foram comparados ao longo de 20 semanas. Nesse novo produto a gordura vegetal hidrogenada, respons<c3><a1>vel pela fixa<c3><a7><c3><a3>o do aroma do produto, foi substitu<c3><ad>da, totalmente ou parcialmente, por <c3><b3>leo de canola. Ao todo foram produzidas 750 observa<c3><a7><c3><b5>es, referentes a 15 avalia<c3><a7><c3><b5>es para cada tipo de snack a cada 2 semanas.

Format

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

semana

Semana da avalia<c3><a7><c3><a3>o.

tipo

Tipo de snack avaliado. Os n<c3><ad>veis representam as seguintes configura<c3><a7><c3><b5>es: A: 22% de gordura, 0% de <c3><b3>leo de canola, B: 0% de gordura, 22% de <c3><b3>leo de canola, C: 17% de gordura, 5% de <c3><b3>leo de canola, D: 11% de gordura, 11% de <c3><b3>leo de canola e E: 5% de gordura, 17% de <c3><b3>leo de canola.

fnpc

For<c3><a7>a necess<c3><a1>ria para o cisalhamento.

Source

PAULA (2004), Exemplo 2.8.1, p<c3><a1>g. 150; Exemplo 2.9.3, p<c3><a1>g. 169.

Examples

data(PaulaEg2.8.1)
str(PaulaEg2.8.1)

# Experimento balanceado, 15 observa<c3><a7><c3><b5>es para cada tipo em cada
# semana
ftable(PaulaEg2.8.1[, c("tipo", "semana")])
xtabs(fnpc ~ tipo + semana, data = PaulaEg2.8.1)

library(lattice)

bwplot(fnpc ~ tipo | factor(semana), pch = "|",
       data =  PaulaEg2.8.1,
       as.table = TRUE,
       strip = strip.custom(strip.names = TRUE,
                            var.name = "semana"))

# Estat<c3><ad>sticas descritivas
resumo <- aggregate(fnpc ~ tipo + semana,
                    data = PaulaEg2.8.1,
                    FUN = function(x) {
                        c("M<c3><a9>dia" = mean(x),
                          "D.Padr<c3><a3>o" = sd(x),
                          "C.Varia<c3><a7><c3><a3>o" = sd(x)/mean(x),
                          "n" = length(x))
                    })
ftable(xtabs(fnpc ~ tipo + semana, data = resumo))

xyplot(fnpc[, "M<c3><a9>dia"] ~ semana,
       groups = tipo,
       data = resumo,
       type = c("l", "g"),
       auto.key = list(
           points = FALSE,
           lines = TRUE,
           title = "snack",
           cex.title = 1.1,
           corner = c(0.1, 0.9)))

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