setwd("F:/2WS10") # change working directory

diamdata <- read.table("sutherlanddata.r", header=TRUE)

# initial analysis


qcc(diamdata[1:40,c("X1","X2","X3","X4","X5")],sizes=5,type="R")

qcc(diamdata[1:40,c("X1","X2","X3","X4","X5")],sizes=5,type="xbar")

# improved analysis after removing operator special cause

qcc(diamdata[41:80,c("X1","X2","X3","X4","X5")],sizes=5,type="R")

qcc(diamdata[41:80,c("X1","X2","X3","X4","X5")],sizes=5,type="xbar")


# final analysis after warm up special cause removed


qcc(diamdata[41:80,c("X1","X2","X3","X4","X5")],sizes=5,type="R")

qcc(diamdata[81:120,c("X1","X2","X3","X4","X5")],sizes=5,type="xbar")