R语言实现质量控制图

上次费了半天劲写出了Xbar-R代码之后第二天就发现了R语言中的SixSigma和qcc程序包。各种质量控制图应有尽有，画的还比我好看#手动微笑。

> library(qcc)  #加载qcc程序包

> qcc(a[1:20,], type="R") 

> a <- a[-7,]  #删去第7组异常值

> qcc(a[1:19,], type="R")  #注意删去之后样本数变为19

此时R图判稳。接下来对Xbar图操作只要换为type="xbar"即可。最终得到的控制用控制图如下：

关于组均值方差之类的统计量也一应俱全：

pcc函数可以实现的控制图如下。有兴趣的可以看Six Sigma with R这本书。