2

我有一个简短的 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 图表创建后立即进行一些过程能力分析（在第 5 页上的此处(PDF) 进行了描述）。为了创建分析图表，我需要将之前创建的 XBar 图表结果中的 LCL 和 UCL 结果存储为变量。有什么办法可以做到这一点吗？

4

1 回答 1

5

?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]

于 2011-05-17T22:06:28.620 回答