访问控制图导致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)

我希望能够在创建XBar图表后立即进行一些过程能力分析(在此处描述(PDF),第5页).为了创建分析图表,我需要将之前创建的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]