#Gage R&R Study #J. Hemmingway
require(qualityTools)## Loading required package: qualityTools## Loading required package: Rsolnp## Loading required package: MASS##
## Attaching package: 'qualityTools'## The following object is masked from 'package:stats':
##
## sigmagageDesign<-read.csv('gagerr.csv', header = TRUE)
design = gageRRDesign ( Operators = 2, Parts = 6, Measurements= 15, randomize = F) ; design## StandardOrder RunOrder Operator Part Measurement
## 1 1 1 A A NA
## 2 2 2 B A NA
## 3 3 3 A B NA
## 4 4 4 B B NA
## 5 5 5 A C NA
## 6 6 6 B C NA
## 7 7 7 A D NA
## 8 8 8 B D NA
## 9 9 9 A E NA
## 10 10 10 B E NA
## 11 11 11 A F NA
## 12 12 12 B F NA
## 13 13 13 A A NA
## 14 14 14 B A NA
## 15 15 15 A B NA
## 16 16 16 B B NA
## 17 17 17 A C NA
## 18 18 18 B C NA
## 19 19 19 A D NA
## 20 20 20 B D NA
## 21 21 21 A E NA
## 22 22 22 B E NA
## 23 23 23 A F NA
## 24 24 24 B F NA
## 25 25 25 A A NA
## 26 26 26 B A NA
## 27 27 27 A B NA
## 28 28 28 B B NA
## 29 29 29 A C NA
## 30 30 30 B C NA
## 31 31 31 A D NA
## 32 32 32 B D NA
## 33 33 33 A E NA
## 34 34 34 B E NA
## 35 35 35 A F NA
## 36 36 36 B F NA
## 37 37 37 A A NA
## 38 38 38 B A NA
## 39 39 39 A B NA
## 40 40 40 B B NA
## 41 41 41 A C NA
## 42 42 42 B C NA
## 43 43 43 A D NA
## 44 44 44 B D NA
## 45 45 45 A E NA
## 46 46 46 B E NA
## 47 47 47 A F NA
## 48 48 48 B F NA
## 49 49 49 A A NA
## 50 50 50 B A NA
## 51 51 51 A B NA
## 52 52 52 B B NA
## 53 53 53 A C NA
## 54 54 54 B C NA
## 55 55 55 A D NA
## 56 56 56 B D NA
## 57 57 57 A E NA
## 58 58 58 B E NA
## 59 59 59 A F NA
## 60 60 60 B F NA
## 61 61 61 A A NA
## 62 62 62 B A NA
## 63 63 63 A B NA
## 64 64 64 B B NA
## 65 65 65 A C NA
## 66 66 66 B C NA
## 67 67 67 A D NA
## 68 68 68 B D NA
## 69 69 69 A E NA
## 70 70 70 B E NA
## 71 71 71 A F NA
## 72 72 72 B F NA
## 73 73 73 A A NA
## 74 74 74 B A NA
## 75 75 75 A B NA
## 76 76 76 B B NA
## 77 77 77 A C NA
## 78 78 78 B C NA
## 79 79 79 A D NA
## 80 80 80 B D NA
## 81 81 81 A E NA
## 82 82 82 B E NA
## 83 83 83 A F NA
## 84 84 84 B F NA
## 85 85 85 A A NA
## 86 86 86 B A NA
## 87 87 87 A B NA
## 88 88 88 B B NA
## 89 89 89 A C NA
## 90 90 90 B C NA
## 91 91 91 A D NA
## 92 92 92 B D NA
## 93 93 93 A E NA
## 94 94 94 B E NA
## 95 95 95 A F NA
## 96 96 96 B F NA
## 97 97 97 A A NA
## 98 98 98 B A NA
## 99 99 99 A B NA
## 100 100 100 B B NA
## 101 101 101 A C NA
## 102 102 102 B C NA
## 103 103 103 A D NA
## 104 104 104 B D NA
## 105 105 105 A E NA
## 106 106 106 B E NA
## 107 107 107 A F NA
## 108 108 108 B F NA
## 109 109 109 A A NA
## 110 110 110 B A NA
## 111 111 111 A B NA
## 112 112 112 B B NA
## 113 113 113 A C NA
## 114 114 114 B C NA
## 115 115 115 A D NA
## 116 116 116 B D NA
## 117 117 117 A E NA
## 118 118 118 B E NA
## 119 119 119 A F NA
## 120 120 120 B F NA
## 121 121 121 A A NA
## 122 122 122 B A NA
## 123 123 123 A B NA
## 124 124 124 B B NA
## 125 125 125 A C NA
## 126 126 126 B C NA
## 127 127 127 A D NA
## 128 128 128 B D NA
## 129 129 129 A E NA
## 130 130 130 B E NA
## 131 131 131 A F NA
## 132 132 132 B F NA
## 133 133 133 A A NA
## 134 134 134 B A NA
## 135 135 135 A B NA
## 136 136 136 B B NA
## 137 137 137 A C NA
## 138 138 138 B C NA
## 139 139 139 A D NA
## 140 140 140 B D NA
## 141 141 141 A E NA
## 142 142 142 B E NA
## 143 143 143 A F NA
## 144 144 144 B F NA
## 145 145 145 A A NA
## 146 146 146 B A NA
## 147 147 147 A B NA
## 148 148 148 B B NA
## 149 149 149 A C NA
## 150 150 150 B C NA
## 151 151 151 A D NA
## 152 152 152 B D NA
## 153 153 153 A E NA
## 154 154 154 B E NA
## 155 155 155 A F NA
## 156 156 156 B F NA
## 157 157 157 A A NA
## 158 158 158 B A NA
## 159 159 159 A B NA
## 160 160 160 B B NA
## 161 161 161 A C NA
## 162 162 162 B C NA
## 163 163 163 A D NA
## 164 164 164 B D NA
## 165 165 165 A E NA
## 166 166 166 B E NA
## 167 167 167 A F NA
## 168 168 168 B F NA
## 169 169 169 A A NA
## 170 170 170 B A NA
## 171 171 171 A B NA
## 172 172 172 B B NA
## 173 173 173 A C NA
## 174 174 174 B C NA
## 175 175 175 A D NA
## 176 176 176 B D NA
## 177 177 177 A E NA
## 178 178 178 B E NA
## 179 179 179 A F NA
## 180 180 180 B F NAresponse(design) = gageDesign$Measure
gdo = gageRR(design, tolerance = .70)##
## AnOVa Table - crossed Design
## Df Sum Sq Mean Sq F value Pr(>F)
## Operator 1 0.177 0.177 261.30 <2e-16 ***
## Part 5 23.297 4.659 6860.11 <2e-16 ***
## Operator:Part 5 0.195 0.039 57.43 <2e-16 ***
## Residuals 168 0.114 0.001
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## ----------
##
## Gage R&R
## VarComp VarCompContrib Stdev StudyVar StudyVarContrib
## totalRR 0.004773 0.03006 0.0691 0.415 0.1734
## repeatability 0.000679 0.00428 0.0261 0.156 0.0654
## reproducibility 0.004094 0.02578 0.0640 0.384 0.1606
## Operator 0.001539 0.00969 0.0392 0.235 0.0984
## Operator:Part 0.002555 0.01609 0.0505 0.303 0.1269
## Part to Part 0.154011 0.96994 0.3924 2.355 0.9849
## totalVar 0.158784 1.00000 0.3985 2.391 1.0000
## P/T Ratio
## totalRR 0.592
## repeatability 0.223
## reproducibility 0.548
## Operator 0.336
## Operator:Part 0.433
## Part to Part 3.364
## totalVar 3.416
##
## ---
## * Contrib equals Contribution in %
## **Number of Distinct Categories (truncated signal-to-noise-ratio) = 8summary(gdo)##
## Operators: 2 Parts: 6
## Measurements: 15 Total: 180
## ----------
##
## AnOVa Table - crossed Design
## Df Sum Sq Mean Sq F value Pr(>F)
## Operator 1 0.177 0.177 261.30 <2e-16 ***
## Part 5 23.297 4.659 6860.11 <2e-16 ***
## Operator:Part 5 0.195 0.039 57.43 <2e-16 ***
## Residuals 168 0.114 0.001
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## ----------
##
## Gage R&R
## VarComp VarCompContrib Stdev StudyVar StudyVarContrib
## totalRR 0.004773 0.03006 0.0691 0.415 0.1734
## repeatability 0.000679 0.00428 0.0261 0.156 0.0654
## reproducibility 0.004094 0.02578 0.0640 0.384 0.1606
## Operator 0.001539 0.00969 0.0392 0.235 0.0984
## Operator:Part 0.002555 0.01609 0.0505 0.303 0.1269
## Part to Part 0.154011 0.96994 0.3924 2.355 0.9849
## totalVar 0.158784 1.00000 0.3985 2.391 1.0000
## P/T Ratio
## totalRR 0.592
## repeatability 0.223
## reproducibility 0.548
## Operator 0.336
## Operator:Part 0.433
## Part to Part 3.364
## totalVar 3.416
##
## ---
## * Contrib equals Contribution in %
## **Number of Distinct Categories (truncated signal-to-noise-ratio) = 8plot(gdo)