# ========================================
# Control Estadístico de Calidad
# Gráfica c
# Dr. Carlos Téllez Martínez
# Marzo 2015
# Tecnológico de Monterrey Campus Guadalajara
# ========================================

# Lectura de datos
Grafica_c1 <- read.csv("/Carlos Tellez Martinez/SkyDrive/Clases/Control Estadístico del Proceso/Material 2013/Tema 2/Videos R/Grafica c/Grafica_c1.csv")
View(Grafica_c1)
Grafica_c1
##    Muestra Disconformidades
## 1        1               21
## 2        2               24
## 3        3               16
## 4        4               12
## 5        5               15
## 6        6                5
## 7        7               28
## 8        8               20
## 9        9               31
## 10      10               25
## 11      11               20
## 12      12               24
## 13      13               16
## 14      14               19
## 15      15               10
## 16      16               17
## 17      17               13
## 18      18               22
## 19      19               18
## 20      20               39
## 21      21               30
## 22      22               24
## 23      23               16
## 24      24               19
## 25      25               17
## 26      26               15
attach(Grafica_c1)
names(Grafica_c1)
## [1] "Muestra"          "Disconformidades"
# Cálculo de los límites de control
# Promedio de c
Promedio_c <- mean(Disconformidades)
Promedio_c
## [1] 19.84615
# Límites de control
LSC <- Promedio_c+3*sqrt(Promedio_c)
LSC
## [1] 33.21086
LC <- Promedio_c
LC
## [1] 19.84615
LIC <-Promedio_c-3*sqrt(Promedio_c)
LIC
## [1] 6.481447
if(LIC<0){LIC <- 0}
LIC
## [1] 6.481447
# Gráfica de control
opts= "b"
plot(Disconformidades, main="Gráfica c", xlab="Observación", xaxt='n', ylim=c(LIC, LSC+10))
axis(side=1, at=seq(1,length(Muestra),1) ,labels=seq(1,length(Muestra),1))
lines(Disconformidades, type=opts)
text(Disconformidades, labels=Muestra, pos=2) # Se etiquetan las muestras
abline(h=LSC, col="blue")
abline(h=LC, col="red")
abline(h=LIC, col="blue")

# Uso del paquete qcc
library("qcc")
## Package 'qcc', version 2.6
## Type 'citation("qcc")' for citing this R package in publications.
Grafica_c <- qcc(Disconformidades,type="c")

# Puntos fuera de control
Grafica_c$violations
## $beyond.limits
## [1] 20  6
## 
## $violating.runs
## numeric(0)