从R中使用GageRR / GageRRDesign生成的数据帧中提取

时间:2019-06-10 14:20:59

标签: r dataframe

如何在使用R中的gageRR函数生成的数据框中提取“ VarCompContrib”列?

这是针对测量系统的GageRR分析。我正在尝试制作一个非常用户友好的程序,其他人可以在其中输入所需的信息,例如操作员,零件和测量的数量以及测量本身,然后输出正确的分析结果。稍后,我将使用if语句执行“分析”部分,但实际上在管理由gageRR生成的数据帧时遇到了麻烦。

library(MASS)
library(Rsolnp)
library(qualityTools)

design = gageRRDesign(Operators=3, Parts=10, Measurements=2, randomize=FALSE)

response(design) = c(23,22,22,22,22,25,23,22,23,22,20,22,22,22,24,25,27,28,      
23,24,23,24,24,22,22,22,24,23,22,24,20,20,25,24,22,24,21,20,21,22,21,22,21,
                     21,24,27,25,27,23,22,25,23,23,22,22,23,25,21,24,23)
gdo=gageRR(design)

plot(gdo)

我希望在VarCompContrib下获得一个7位数的列向量

1 个答案:

答案 0 :(得分:0)

对于初学者来说,您可以使用gdo来查看str(gdo)的结构。从那里,我们看到Varcomp是一个插槽,因此我们可以使用gdo@Varcomp访问它,并将其转换为data.frame

library(qualityTools)

design <- gageRRDesign(Operators = 3, Parts = 10, Measurements = 2, randomize = FALSE)

response(design) <- c(
  23,22,22,22,22,25,23,22,23,22,20,22,22,22,24,25,27,28,23,24,23,24,24,22,22,22,24,23,22,24,
  20,20,25,24,22,24,21,20,21,22,21,22,21,21,24,27,25,27,23,22,25,23,23,22,22,23,25,21,24,23
)

gdo <- gageRR(design)

data.frame(gdo@Varcomp)
#   totalRR repeatability reproducibility         a a_b     bTob totalVar
# 1 1.66441      1.209028       0.4553819 0.4553819   0 1.781211 3.445621