# UNIVERSIDAD NACIONAL MAYOR DE SAN MARCOS
# UNIVERSIDAD DEL PERU, DECANA DE AMERICA
# FACULTAD DE INGENIERIA INDUSTRIAL
# CURSO: ESTADISTICA INDUSTRIAL
# EJERCICIO 03
library(car)
## Warning: package 'car' was built under R version 4.1.3
## Loading required package: carData
## Warning: package 'carData' was built under R version 4.1.3
library(ggplot2)
## Warning: package 'ggplot2' was built under R version 4.1.3
library(readxl)
## Warning: package 'readxl' was built under R version 4.1.3
# Ingreso de Datos
Datos <- read_excel("E:/UNMSM LIMA/ESTADISTICA INDUSTRIAL/UNIDAD I/TAREA 02/Datos.xlsx")
Datos
## # A tibble: 32 x 2
## Resistencia Plastico
## <dbl> <chr>
## 1 135 p1
## 2 175 p1
## 3 97 p1
## 4 169 p1
## 5 213 p1
## 6 171 p1
## 7 115 p1
## 8 143 p1
## 9 275 p2
## 10 170 p2
## # ... with 22 more rows
# Convertir de character a factor
Datos$Plastico <- factor(Datos$Plastico)
Datos$Plastico
## [1] p1 p1 p1 p1 p1 p1 p1 p1 p2 p2 p2 p2 p2 p2 p2 p2 p3 p3 p3 p3 p3 p3 p3 p3 p4
## [26] p4 p4 p4 p4 p4 p4 p4
## Levels: p1 p2 p3 p4
# Modelo
mod <- aov(Resistencia ~ Plastico, data = Datos)
mod
## Call:
## aov(formula = Resistencia ~ Plastico, data = Datos)
##
## Terms:
## Plastico Residuals
## Sum of Squares 69072.12 37410.75
## Deg. of Freedom 3 28
##
## Residual standard error: 36.55268
## Estimated effects may be unbalanced
# Resultado
summary(mod)
## Df Sum Sq Mean Sq F value Pr(>F)
## Plastico 3 69072 23024 17.23 1.55e-06 ***
## Residuals 28 37411 1336
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
# Grafico-Cajas
boxplot(Resistencia~Plastico,data=Datos, col=c("red","green","blue","orange"), main="Diagrama de la Calidad de los plasticos utilizados")
mtext("By:BSP", side = 3, adj = 1, family = "mono")
