Trabajo Práctico de Control de Calidad

mcarboni-cc-tp.R

# ACTIVIDAD 1
library(qcc)
library(SixSigma)
library(moments)

granulos_str <- "SUBGRUPO,MEDIA 
1,214.18 
2,213.48 
3,213.98 
4,214.12 
5,214.46 
6,213.38 
7,231.56 
8,214.08 
9,231.72 
10,214.64 
11,213.92 
12,213.96 
13,214.20 
14,213.74 
15,214.26 
16,214.18 
17,214.00 
18,213.60 
19,214.20 
20,214.38 
21,213.78 
22,213.74 
23,213.32 
24,214.02 
25,214.24"

granulos=read.table(text=granulos_str,header=T,dec=".",sep=",")

summary(granulos$MEDIA)

par(mfrow=c(1,1))
hist(granulos$MEDIA,nclass=4,col=6, xlab = 'Media en mg', ylab = 'Frecuencia', main = 'Media de Muestra de Cuatro Gránulos de PVC')
par(mfrow=c(1,2))
boxplot(granulos$MEDIA,col=3, ylab = 'Media en mg', main = 'Boxplot de las muestras')
qqnorm(granulos$MEDIA)
qqline(granulos$MEDIA)

par(mfrow=c(1,1))
granulos.chart <- qcc(data=granulos$MEDIA, type="xbar.one")
summary(granulos.chart)

granulos.limpio <- granulos[-c(7,9),]
par(mfrow=c(1,3))
hist(granulos.limpio$MEDIA,nclass=4,col=6, xlab = 'Media en mg', ylab = 'Frecuencia', main = 'Media de Muestra de Cuatro Gránulos de PVC')
boxplot(granulos.limpio$MEDIA,col=3, ylab = 'Media en mg', main = 'Boxplot de las muestras')
qqnorm(granulos.limpio$MEDIA)
qqline(granulos.limpio$MEDIA)
granulos.limpio.chart <- qcc(data=granulos.limpio$MEDIA, type="xbar.one")
summary(granulos.limpio.chart)

process.capability(object = granulos.limpio.chart,spec.limits = c(212,218) )

# ACTIVIDAD 2

presion_grupos_str = "SUBGRUPO;a1;a2;a3
1;315,6;319,2;303,8
2;318,8;309,2;321,4
3;311,2;312,1;342,9
4;322;321,1;329,1
5;315,2;327,4;300,6
6;310,3;319,8;338,5
7;320,6;315,9;318,3
8;322,2;303,6;323,4
9;329,1;306,7;312,4
10;322,4;318,8;299,7
11;326,2;310,1;338,5
12;328,8;325;322
13;328,8;306,3;305,6
14;318,7;320,8;310,3
15;326,7;316,7;327,3
16;313,4;307,4;329,5
17;337,3;312,9;324,4
18;316,3;314,1;323
19;327,2;338,2;340,9
20;337,8;343;337,4
21;309,2;321,7;310,5
22;314,3;321,6;318
23;318,9;322,2;333,5
24;303,7;326,3;337,1
25;319,3;338,8;320,9
26;317;327,4;312,5
27;310,6;318,5;336,7
28;319,5;326;333,2
29;308,6;321,7;306
30;316,2;321,6;328,5"

presion_grupos=read.table(text=presion_grupos_str,header=T,dec=",",sep=";")

presion_medias <- rowMeans(presion_grupos[,2:4])

par(mfrow=c(1,3))
hist(presion_medias,nclass=4,col=6, xlab = 'Media en psi', ylab = 'Frecuencia', main = 'Resistencia de los artículos')
boxplot(presion_medias,col=3, ylab = 'Media en psi', main = 'Boxplot de las medias')
qqnorm(presion_medias)
qqline(presion_medias)

skewness(presion_medias)

ps <- presion_grupos[,2:4]
presion_grupos.chart <- qcc(data = ps, type = "xbar")
summary(presion_grupos.chart)

process.capability(object = presion_grupos.chart,spec.limits = c(300,NA) )

presion_grupos.limpio <- presion_grupos[-20,]
ps <- presion_grupos.limpio[,2:4]
presion_grupos.chart <- qcc(data = ps, type = "xbar")
summary(presion_grupos.chart)

process.capability(object = presion_grupos.chart,spec.limits = c(300,NA) )

presion_grupos.limpio <- presion_grupos[-c(19,20),]
ps <- presion_grupos.limpio[,2:4]
presion_grupos.chart <- qcc(data = ps, type = "xbar")
summary(presion_grupos.chart)

process.capability(object = presion_grupos.chart,spec.limits = c(300,NA) )

presion_medias <- rowMeans(presion_grupos.limpio[,2:4])
presion_medias.chart <- qcc(data=presion_medias, type="xbar.one")
process.capability(object = presion_medias.chart,spec.limits = c(300,NA) )