Introdução

Este relatório analisa o tempo de dissolução de comprimidos de Omeprazol em diferentes meios (Água, Limão, Vinagre e Bicarbonato).

# 1. Definindo os dados e tratamentos
y <- c(60,69,70,61,70, 40,41,41,42,42, 44,41,42,44,43, 45,45,46,49,48)
trat <- factor(rep(c("Agua", "Limao", "Vinagre", "Bicarbonato"), each=5))
dados <- data.frame(trat, y)

# 2. Mostrando as médias
tapply(y, trat, mean)
##        Agua Bicarbonato       Limao     Vinagre 
##        66.0        46.6        41.2        42.8
boxplot(y ~ trat, data = dados, 
        col = c("skyblue","lightgreen","salmon","khaki"),
        main = "Dissolução do Omeprazol por Meio",
        xlab = "Tratamentos", 
        ylab = "Tempo (segundos)")

modelo <- aov(y ~ trat, data = dados)
summary(modelo)
##             Df Sum Sq Mean Sq F value   Pr(>F)    
## trat         3 1969.7   656.6   84.18 5.16e-10 ***
## Residuals   16  124.8     7.8                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
# Teste de Normalidade (Shapiro-Wilk)
shapiro.test(residuals(modelo))
## 
##  Shapiro-Wilk normality test
## 
## data:  residuals(modelo)
## W = 0.94624, p-value = 0.3136
if(!require(agricolae)) install.packages("agricolae")
## Carregando pacotes exigidos: agricolae
library(agricolae)

teste_tukey <- HSD.test(modelo, "trat", group = TRUE)
print(teste_tukey$groups)
##                y groups
## Agua        66.0      a
## Bicarbonato 46.6      b
## Vinagre     42.8     bc
## Limao       41.2      c