GRUPO SULFATOS: LISBETH BELIZARIO, JULITZA CUTIPA, RUSSELL QUISPE

file.choose()

[1] "C:\\Users\\Lisbeth\\Downloads\\DISEÑO EXPERIMENTAL (TRABAJOS)\\SULFATOS.xlsx"

ruta_propulsora <- "C:\\Users\\Lisbeth\\Downloads\\DISEÑO EXPERIMENTAL (TRABAJOS)\\SULFATOS.xlsx"

excel_sheets(ruta_propulsora)

[1] "Hoja1"

Propul<-read_excel(ruta_propulsora)

print(head(Propul))

view(Propul)

attach(Propul)

names(Propul)

 [1] "Año"             "Tipo de agua"    "COD. Muestra"    "Distrito"       
 [5] "Provincia"       "Aspecto"         "Color"           "Olor"           
 [9] "Sabor"           "sulfatos (mg/L)"

summary(Propul)

      Año       Tipo de agua       COD. Muestra         Distrito        
 Min.   :2012   Length:1718        Length:1718        Length:1718       
 1st Qu.:2014   Class :character   Class :character   Class :character  
 Median :2017   Mode  :character   Mode  :character   Mode  :character  
 Mean   :2017                                                           
 3rd Qu.:2021                                                           
 Max.   :2027                                                           
  Provincia           Aspecto             Color               Olor          
 Length:1718        Length:1718        Length:1718        Length:1718       
 Class :character   Class :character   Class :character   Class :character  
 Mode  :character   Mode  :character   Mode  :character   Mode  :character  
                                                                            
                                                                            
                                                                            
    Sabor           sulfatos (mg/L)   
 Length:1718        Length:1718       
 Class :character   Class :character  
 Mode  :character   Mode  :character  
                                      
                                      
                                      

str(Propul)

tibble [1,718 × 10] (S3: tbl_df/tbl/data.frame)
 $ Año            : num [1:1718] 2012 2012 2012 2012 2012 ...
 $ Tipo de agua   : chr [1:1718] "Agua potable" "Agua potable" "Agua potable" "Agua potable" ...
 $ COD. Muestra   : chr [1:1718] "2903" "2936" "3104" "2977" ...
 $ Distrito       : chr [1:1718] "Chupa" "0" "Ilave" "Orurillo" ...
 $ Provincia      : chr [1:1718] "Azangaro" "Carabaya" "Ilave" "Melgar" ...
 $ Aspecto        : chr [1:1718] "Liquido" "Liquido" "Solido" "Liquido" ...
 $ Color          : chr [1:1718] "Incoloro" "Incoloro" "Café" "Incoloro" ...
 $ Olor           : chr [1:1718] "Inodoro" "Inodoro" "Inodoro" "Inodoro" ...
 $ Sabor          : chr [1:1718] "Insipido" "Insipido" "Insipido" "Insipido" ...
 $ sulfatos (mg/L): chr [1:1718] "4" "4" "4" "760" ...

YEAR <- factor(Propul$`Año`) 

TYPEOFWATER <- factor(Propul$`Tipo de agua`)

COD <- factor(Propul$`COD. Muestra`)

DISTRIC <- factor(Propul$`Distrito`)

PROVINCE <- factor(Propul$`Provincia`)

ASPECT <- factor(Propul$`Aspecto`)

COLOR <- factor(Propul$`Color`)

SMELL <- factor(Propul$`Olor`)

FLAVOR <- factor(Propul$`Sabor`)

SULFATES <- factor(Propul$`sulfatos (mg/L)`)

SULFATES1<-as.numeric(SULFATES)

par(mfrow=c(1,1))

boxplot(split(SULFATES1,YEAR),xlab="Año", ylab="Sulfatos") 

boxplot(split(SULFATES1,TYPEOFWATER),xlab="Tipodeagua", ylab="Sulfatos") 

boxplot(split(SULFATES1,COD),xlab="COD.Muestra", ylab="Sulfatos") 

boxplot(split(SULFATES1,DISTRIC),xlab="Distrito", ylab="Sulfatos") 

boxplot(split(SULFATES1,PROVINCE),xlab="Provincia", ylab="Sulfatos")

boxplot(split(SULFATES1,ASPECT),xlab="Aspecto", ylab="Sulfatos")

boxplot(split(SULFATES1,COLOR),xlab="Color", ylab="Sulfatos")

boxplot(split(SULFATES1,SMELL),xlab="Olor", ylab="Sulfatos")

boxplot(split(SULFATES1,FLAVOR),xlab="Sabor", ylab="Sulfatos")

#Análisis de varianza usando la función (aov) Analysis of Variance

prop.aov<-aov(SULFATES1 ~ YEAR+TYPEOFWATER+COD+DISTRIC+PROVINCE+ASPECT+COLOR+SMELL+FLAVOR+SULFATES)

anova(prop.aov)

Warning: ANOVA F-tests on an essentially perfect fit are unreliable

Analysis of Variance Table

Response: SULFATES1
             Df   Sum Sq Mean Sq    F value    Pr(>F)    
YEAR          9  5584148  620461 5.8219e+27 < 2.2e-16 ***
TYPEOFWATER  86 17795307  206922 1.9416e+27 < 2.2e-16 ***
COD         935 45584760   48754 4.5746e+26 < 2.2e-16 ***
DISTRIC      74  3988277   53896 5.0571e+26 < 2.2e-16 ***
PROVINCE     12   598659   49888 4.6811e+26 < 2.2e-16 ***
ASPECT        3    65074   21691 2.0353e+26 < 2.2e-16 ***
COLOR         9   432157   48017 4.5055e+26 < 2.2e-16 ***
SMELL         7   419673   59953 5.6255e+26 < 2.2e-16 ***
FLAVOR        3    29946    9982 9.3661e+25 < 2.2e-16 ***
SULFATES    232 12446708   53650 5.0340e+26 < 2.2e-16 ***
Residuals   326        0       0                         
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1