在R中訪問控制圖結果？

Question

我有一個簡短的R腳本來加載一堆數據並將它繪製在XBar圖表中。使用下面的代碼，我可以繪製數據並查看各種統計信息。

library(qcc) 
tir<-read.table("data.dat", header=T,,sep="\t") 
names(tir) 
attach(tir) 
rand <- sample(tir) 
xbarchart <- qcc(rand[1:100,],type="R") 
summary(xbarchart) 

我希望能夠做一些過程能力分析（描述here（5頁PDF））創建X條圖之後。爲了創建分析圖表，我需要將前面創建的XBar圖表結果中的LCL和UCL結果存儲爲變量。有什麼辦法可以做到這一點？

Answer 1

我將使用?qcc幫助文件中的示例回答您的問題。

x <- c(33.75, 33.05, 34, 33.81, 33.46, 34.02, 33.68, 33.27, 33.49, 33.20, 
    33.62, 33.00, 33.54, 33.12, 33.84) 

xbarchart <- qcc(x, type="xbar.one", std.dev = "SD") 

一個有用的功能檢查的變量和函數結果的結構是str()，簡稱結構。

str(xbarchart) 

List of 11 
$ call  : language qcc(data = x, type = "xbar.one", std.dev = "SD") 
$ type  : chr "xbar.one" 
$ data.name : chr "x" 
$ data  : num [1:15, 1] 33.8 33 34 33.8 33.5 ... 
    ..- attr(*, "dimnames")=List of 2 
    .. ..$ Group : chr [1:15] "1" "2" "3" "4" ... 
    .. ..$ Samples: NULL 
$ statistics: Named num [1:15] 33.8 33 34 33.8 33.5 ... 
    ..- attr(*, "names")= chr [1:15] "1" "2" "3" "4" ... 
$ sizes  : int [1:15] 1 1 1 1 1 1 1 1 1 1 ... 
$ center : num 33.5 
$ std.dev : num 0.342 
$ nsigmas : num 3 
$ limits : num [1, 1:2] 32.5 34.5 
    ..- attr(*, "dimnames")=List of 2 
    .. ..$ : chr "" 
    .. ..$ : chr [1:2] "LCL" "UCL" 
$ violations:List of 2 
    ..$ beyond.limits : int(0) 
    ..$ violating.runs: num(0) 
- attr(*, "class")= chr "qcc" 

您會注意到此列表中倒數第二個元素稱爲$ limits，幷包含LCL和UCL的兩個值。

它是簡單的，以提取這些元素：

limits <- xbarchart$limits 
limits 

     LCL  UCL 
32.49855 34.54811 

因此LCL <- limits[1]和UCL <- limits[2]