---

title: “segundo3 Anova” author: “zoilomaximo” date: “5/11/2020” output: html_document

library(readxl)
tela <- read_excel("tela.xlsx")
summary(tela)

##   sustancia           muestra           resistencia   
##  Length:20          Length:20          Min.   :0.400  
##  Class :character   Class :character   1st Qu.:1.275  
##  Mode  :character   Mode  :character   Median :1.750  
##                                        Mean   :1.960  
##                                        3rd Qu.:2.250  
##                                        Max.   :4.400

attach(tela)
sustancia <- as.factor(sustancia)
summary(sustancia)

## Quimico 1 Quimico 2 Quimico 3 Quimico 4 
##         5         5         5         5

plot tratamientos

Diagrama de cajas para los tratamientos

bloques

ggplot(tela, aes(x = muestra, y = resistencia)) +
  geom_boxplot(fill = "grey80", colour = "blue") +
  scale_x_discrete() + xlab("Muestra de tela") +
  ylab("Resistencia")

Supuestos a cumplir en anova

anova.tela = aov(resistencia ~ sustancia+muestra,
                 data=tela)
summary(anova.tela)

##             Df Sum Sq Mean Sq F value   Pr(>F)    
## sustancia    3 18.044   6.015   75.89 4.52e-08 ***
## muestra      4  6.693   1.673   21.11 2.32e-05 ***
## Residuals   12  0.951   0.079                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Normalidad en base a residuales

si p.valor > 0,05 se cumple supuesto Con base en la prueba Shapiro wilk se cumple el supuesto

e<-anova.tela$residuals ##### residuales

shapiro.test(e) #### Normalidad

## 
##  Shapiro-Wilk normality test
## 
## data:  e
## W = 0.98077, p-value = 0.9437

hist(e, freq=FALSE)
curve(dnorm(x,mean(e), sd(e)), xlim=c(-0.7,0.7), add=TRUE, col=2)

homogeneidad de varianzas

library(carData)
library(car)

leveneTest(e ~ as.factor(sustancia), data = tela, center = "median") #####homogeneidad

## Levene's Test for Homogeneity of Variance (center = "median")
##       Df F value Pr(>F)
## group  3   0.516 0.6771
##       16

Anova DBCA

library(carData)
library(car)

analisis post anova

library(agricolae) #### PAQUETE DE PRUEBAS POST-ANOVA
HSD.test(anova.tela,"sustancia", alpha=0.05,
         console=TRUE, group=FALSE)

## 
## Study: anova.tela ~ "sustancia"
## 
## HSD Test for resistencia 
## 
## Mean Square Error:  0.07925 
## 
## sustancia,  means
## 
##           resistencia       std r Min Max
## Quimico 1        1.14 0.4037326 5 0.5 1.6
## Quimico 2        1.76 0.7924645 5 0.4 2.4
## Quimico 3        1.38 0.4764452 5 0.6 1.8
## Quimico 4        3.56 0.9449868 5 2.0 4.4
## 
## Alpha: 0.05 ; DF Error: 12 
## Critical Value of Studentized Range: 4.19866 
## 
## Comparison between treatments means
## 
##                       difference pvalue signif.        LCL         UCL
## Quimico 1 - Quimico 2      -0.62 0.0204       * -1.1485978 -0.09140218
## Quimico 1 - Quimico 3      -0.24 0.5523         -0.7685978  0.28859782
## Quimico 1 - Quimico 4      -2.42 0.0000     *** -2.9485978 -1.89140218
## Quimico 2 - Quimico 3       0.38 0.1973         -0.1485978  0.90859782
## Quimico 2 - Quimico 4      -1.80 0.0000     *** -2.3285978 -1.27140218
## Quimico 3 - Quimico 4      -2.18 0.0000     *** -2.7085978 -1.65140218