Homework 13

2. Syringe Caps

library(qcc)
data.dir = "http://media.pearsoncmg.com/cmg/pmmg_mml_shared/mathstatsresources/Akritas"
cap = read.table(file.path(data.dir,"SqcSyringeL.txt"),header=T)

cap.groups = qcc.groups(cap$x,cap$sample)

#?stats.xbar
# "UWAVE-R" is the default 


### This is wrong
qcc(cap.groups,type="xbar",std.dev="UWAVE-R")$call

## qcc(data = cap.groups, type = "xbar", std.dev = "UWAVE-R")

# a
qcc(cap.groups[c(1:15),],type="xbar",newdata=cap.groups[16:47,],std.dev="UWAVE-R")$call

## qcc(data = cap.groups[c(1:15), ], type = "xbar", std.dev = "UWAVE-R", 
##     newdata = cap.groups[16:47, ])

# b
qcc(cap.groups[c(1,3:15),],type="xbar",newdata=cap.groups[16:47,],std.dev="UWAVE-R")$call

## qcc(data = cap.groups[c(1, 3:15), ], type = "xbar", std.dev = "UWAVE-R", 
##     newdata = cap.groups[16:47, ])

#d
qcc(cap.groups[c(1:15),],type="xbar",newdata=cap.groups[16:47,],std.dev="UWAVE-S")$call

## qcc(data = cap.groups[c(1:15), ], type = "xbar", std.dev = "UWAVE-S", 
##     newdata = cap.groups[16:47, ])

qcc(cap.groups[c(1,3:15),],type="xbar",newdata=cap.groups[16:47,],std.dev="UWAVE-S")$call

## qcc(data = cap.groups[c(1, 3:15), ], type = "xbar", std.dev = "UWAVE-S", 
##     newdata = cap.groups[16:47, ])

4. construction labor X-chart

data.dir = "http://media.pearsoncmg.com/cmg/pmmg_mml_shared/mathstatsresources/Akritas"
rate = read.table(file.path(data.dir,"SqcLaborProd.txt"),header=T)$x

#a
qcc(rate[1:30],type="xbar.one",newdata=rate[31:length(rate)])$call

## qcc(data = rate[1:30], type = "xbar.one", newdata = rate[31:length(rate)])

#b
qcc(rate[c(1:11,14:30)],type="xbar.one",newdata=rate[31:length(rate)])$call

## qcc(data = rate[c(1:11, 14:30)], type = "xbar.one", newdata = rate[31:length(rate)])