如何在使用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位数的列向量
答案 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