1### R code from vignette source 'qualityTools.rnw'
2### Encoding: UTF-8
3
4###################################################
5### code chunk number 1: qualityTools.rnw:176-177
6###################################################
7library(qualityTools)
8
9
10###################################################
11### code chunk number 2: qualityTools.rnw:181-185
12###################################################
13#create artificial defect data set
14defects = c(rep("E", 62), rep("B", 15), rep("F", 3), rep("A", 10),
15rep("C",20), rep("D", 10))
16paretoChart(defects)
17
18
19###################################################
20### code chunk number 3: qualityTools.rnw:262-266
21###################################################
22x = c(9.991, 10.013, 10.001, 10.007, 10.010, 10.013, 10.008, 10.017, 10.005,
2310.005, 10.002, 10.017, 10.005, 10.002,  9.996, 10.011, 10.009, 10.006,
2410.008, 10.003, 10.002, 10.006, 10.010, 9.992, 10.013)
25cg(x, target = 10.003, tolerance = c(9.903, 10.103))
26
27
28###################################################
29### code chunk number 4: qualityTools.rnw:288-296
30###################################################
31#create a gage RnR design
32design = gageRRDesign(Operators=3, Parts=10, Measurements=2, randomize=FALSE)
33#set the response
34response(design) = c(23,22,22,22,22,25,23,22,23,22,20,22,22,22,24,25,27,28,
3523,24,23,24,24,22,22,22,24,23,22,24,20,20,25,24,22,24,21,20,21,22,21,22,21,
3621,24,27,25,27,23,22,25,23,23,22,22,23,25,21,24,23)
37#perform Gage RnR
38gdo = gageRR(design)
39
40
41###################################################
42### code chunk number 5: GageRR
43###################################################
44#visualization of Gage RnR
45plot(gdo)
46
47
48###################################################
49### code chunk number 6: qualityTools.rnw:306-307
50###################################################
51#visualization of Gage RnR
52plot(gdo)
53
54
55###################################################
56### code chunk number 7: NormPcr
57###################################################
58set.seed(1234)
59#generate some data
60norm = rnorm(20, mean = 20)
61#generate some data
62weib = rweibull(20, shape = 2, scale = 8)
63#process capability
64pcr(norm, "normal", lsl = 17, usl = 23)
65
66
67###################################################
68### code chunk number 8: WeibPcr
69###################################################
70#process cabapility
71pcr(weib, "weibull", usl = 20)
72
73
74###################################################
75### code chunk number 9: qualityTools.rnw:406-407
76###################################################
77set.seed(1234)
78#generate some data
79norm = rnorm(20, mean = 20)
80#generate some data
81weib = rweibull(20, shape = 2, scale = 8)
82#process capability
83pcr(norm, "normal", lsl = 17, usl = 23)
84
85
86###################################################
87### code chunk number 10: qualityTools.rnw:412-413
88###################################################
89#process cabapility
90pcr(weib, "weibull", usl = 20)
91
92
93###################################################
94### code chunk number 11: qualityTools.rnw:425-426
95###################################################
96pcr(weib, "weibull", usl = 20)
97
98
99###################################################
100### code chunk number 12: qqPlot
101###################################################
102par(mfrow = c(1,2))
103qqPlot(weib, "weibull"); qqPlot(weib, "normal")
104
105
106###################################################
107### code chunk number 13: qualityTools.rnw:437-438
108###################################################
109par(mfrow = c(1,2))
110qqPlot(weib, "weibull"); qqPlot(weib, "normal")
111
112
113###################################################
114### code chunk number 14: ppPlot
115###################################################
116par(mfrow = c(1,2))
117ppPlot(norm, "weibull"); ppPlot(norm, "normal")
118
119
120###################################################
121### code chunk number 15: qualityTools.rnw:452-453
122###################################################
123par(mfrow = c(1,2))
124ppPlot(norm, "weibull"); ppPlot(norm, "normal")
125
126
127###################################################
128### code chunk number 16: qualityTools.rnw:491-497
129###################################################
130set.seed(1234)
131fdo = facDesign(k = 3, centerCube = 4) #fdo - factorial design object
132names(fdo) = c("Factor 1", "Factor 2", "Factor 3") #optional
133lows(fdo) = c(80, 120, 1) #optional
134highs(fdo) = c(120, 140, 2) #optional
135summary(fdo) #information about the factorial design
136
137
138###################################################
139### code chunk number 17: qualityTools.rnw:503-505
140###################################################
141#set first value
142yield = simProc(x1 = 120, x2 = 140, x3 = 2)
143
144
145###################################################
146### code chunk number 18: qualityTools.rnw:517-521
147###################################################
148yield = c(simProc(120,140, 1),simProc(80,140, 1),simProc(120,140, 2),
149simProc(120,120, 1),simProc(90,130, 1.5),simProc(90,130, 1.5),
150simProc(80,120, 2),simProc(90,130, 1.5),simProc(90,130, 1.5),
151simProc(120,120, 2),simProc(80,140, 2),simProc(80,120, 1))
152
153
154###################################################
155### code chunk number 19: qualityTools.rnw:527-528
156###################################################
157response(fdo) = yield	#assign yield to the factorial design object
158
159
160###################################################
161### code chunk number 20: effPlot
162###################################################
163effectPlot(fdo, classic = TRUE)
164
165
166###################################################
167### code chunk number 21: iaPlot
168###################################################
169interactionPlot(fdo)
170
171
172###################################################
173### code chunk number 22: qualityTools.rnw:545-546
174###################################################
175effectPlot(fdo, classic = TRUE)
176
177
178###################################################
179### code chunk number 23: qualityTools.rnw:551-552
180###################################################
181interactionPlot(fdo)
182
183
184###################################################
185### code chunk number 24: qualityTools.rnw:564-566
186###################################################
187lm.1 = lm(yield ~ A*B*C, data = fdo)
188summary(lm.1)
189
190
191###################################################
192### code chunk number 25: ParetPlot
193###################################################
194par(mfrow = c(1,2))
195paretoPlot(fdo)
196normalPlot(fdo)
197
198
199###################################################
200### code chunk number 26: qualityTools.rnw:579-580
201###################################################
202par(mfrow = c(1,2))
203paretoPlot(fdo)
204normalPlot(fdo)
205
206
207###################################################
208### code chunk number 27: wirePlot
209###################################################
210par(mfrow = c(1,2))
211wirePlot(A, B, yield, data = fdo)
212contourPlot(A, B, yield, data = fdo)
213
214
215###################################################
216### code chunk number 28: qualityTools.rnw:595-596
217###################################################
218par(mfrow = c(1,2))
219wirePlot(A, B, yield, data = fdo)
220contourPlot(A, B, yield, data = fdo)
221
222
223###################################################
224### code chunk number 29: qualityTools.rnw:612-613
225###################################################
226fdo.frac = fracDesign(k = 3, gen = "C = AB", centerCube = 4)
227
228
229###################################################
230### code chunk number 30: qualityTools.rnw:618-619
231###################################################
232summary(fdo.frac)
233
234
235###################################################
236### code chunk number 31: qualityTools.rnw:632-633
237###################################################
238aliasTable(fdo.frac)
239
240
241###################################################
242### code chunk number 32: qualityTools.rnw:638-639
243###################################################
244confounds(fdo.frac)
245
246
247###################################################
248### code chunk number 33: fracChoose
249###################################################
250fracChoose()
251
252
253###################################################
254### code chunk number 34: qualityTools.rnw:649-650
255###################################################
256fracChoose()
257
258
259###################################################
260### code chunk number 35: qualityTools.rnw:660-661
261###################################################
262fdo1 = facDesign(k = 3, centerCube = 2, replicates = 2)
263
264
265###################################################
266### code chunk number 36: qualityTools.rnw:668-671
267###################################################
268set.seed(1234)
269y2 = rnorm(12, mean = 20)
270response(fdo) = data.frame(yield, y2)
271
272
273###################################################
274### code chunk number 37: multiresp
275###################################################
276par(mfrow = c(1,2))
277wirePlot(A, B, yield, data = fdo, form = "yield~A+B+C+A*B")
278contourPlot(A, B, y2, data = fdo, form = "y2~A+B+C+A*B")
279
280
281###################################################
282### code chunk number 38: qualityTools.rnw:683-684
283###################################################
284par(mfrow = c(1,2))
285wirePlot(A, B, yield, data = fdo, form = "yield~A+B+C+A*B")
286contourPlot(A, B, y2, data = fdo, form = "y2~A+B+C+A*B")
287
288
289###################################################
290### code chunk number 39: multiresp2
291###################################################
292par(mfrow = c(1,2))
293wirePlot(A,B,y2, data = fdo, factors = list(C=-1), form = "y2~A*B*C")
294wirePlot(A,B,y2, data = fdo, factors = list(C=1), form = "y2~A*B*C")
295
296
297###################################################
298### code chunk number 40: qualityTools.rnw:699-700
299###################################################
300par(mfrow = c(1,2))
301wirePlot(A,B,y2, data = fdo, factors = list(C=-1), form = "y2~A*B*C")
302wirePlot(A,B,y2, data = fdo, factors = list(C=1), form = "y2~A*B*C")
303
304
305###################################################
306### code chunk number 41: qualityTools.rnw:708-711
307###################################################
308fits(fdo) = lm(yield ~ A+B, data = fdo)
309fits(fdo) = lm(y2 ~ A*B*C, data = fdo)
310fits(fdo)
311
312
313###################################################
314### code chunk number 42: qualityTools.rnw:718-720
315###################################################
316sao =steepAscent(factors=c("A","B"),response="yield",data=fdo,steps=20)
317sao
318
319
320###################################################
321### code chunk number 43: steepAsc
322###################################################
323predicted = simProc(sao[,5], sao[,6])
324response(sao) = predicted
325plot(sao, type = "b", col = 2)
326
327
328###################################################
329### code chunk number 44: qualityTools.rnw:732-733
330###################################################
331predicted = simProc(sao[,5], sao[,6])
332response(sao) = predicted
333plot(sao, type = "b", col = 2)
334
335
336###################################################
337### code chunk number 45: qualityTools.rnw:747-753
338###################################################
339#set the seed for randomization of the runs
340set.seed(1234)
341fdo2 = facDesign(k = 2, centerCube = 3)
342names(fdo2) = c("Factor 1", "Factor 2")
343lows(fdo2) = c(134, 155)
344highs(fdo2) = c(155, 175)
345
346
347###################################################
348### code chunk number 46: qualityTools.rnw:758-762
349###################################################
350yield = c(simProc(134,175), simProc(144.5,165.5), simProc(155,155),
351simProc(144.5,165.5), simProc(155,175), simProc(144.5,165.5),
352simProc(134,155))
353response(fdo2) = yield
354
355
356###################################################
357### code chunk number 47: qualityTools.rnw:769-771
358###################################################
359rsdo = starDesign(data = fdo2)
360rsdo
361
362
363###################################################
364### code chunk number 48: qualityTools.rnw:776-779
365###################################################
366yield2 = c(yield, simProc(130,165), simProc(155,165), simProc(144,155),
367simProc(144,179),simProc(144,165),simProc(144,165),simProc(144,165))
368response(rsdo) = yield2
369
370
371###################################################
372### code chunk number 49: qualityTools.rnw:784-785
373###################################################
374lm.3 = lm(yield2 ~ A*B + I(A^2) + I(B^2), data = rsdo)
375
376
377###################################################
378### code chunk number 50: respSur
379###################################################
380par(mfrow=c(1,2))
381wirePlot(A,B,yield2,form="yield2~A*B+I(A^2)+I(B^2)",data=rsdo,theta=-70)
382contourPlot(A,B,yield2,form="yield2~A*B+I(A^2)+I(B^2)",data=rsdo)
383
384
385###################################################
386### code chunk number 51: qualityTools.rnw:797-798
387###################################################
388par(mfrow=c(1,2))
389wirePlot(A,B,yield2,form="yield2~A*B+I(A^2)+I(B^2)",data=rsdo,theta=-70)
390contourPlot(A,B,yield2,form="yield2~A*B+I(A^2)+I(B^2)",data=rsdo)
391
392
393###################################################
394### code chunk number 52: qualityTools.rnw:807-811
395###################################################
396A = seq(40, 210, length = 100)
397B = seq(90, 190, length = 100)
398C = seq(90, 190, length = 100)
399filled.contour(A, B,outer(A,B, simProc, noise = FALSE), xlab = "Factor 1", ylab = "Factor 2", color.palette = colorRampPalette(c("#00007F", "blue", "#007FFF", "cyan","#7FFF7F", "yellow", "#FF7F00", "red", "#7F0000")))
400
401
402###################################################
403### code chunk number 53: qualityTools.rnw:821-822
404###################################################
405fdo = rsmDesign(k = 3, alpha = 1.633, cc = 0, cs = 6)
406
407
408###################################################
409### code chunk number 54: qualityTools.rnw:827-828
410###################################################
411fdo = randomize(fdo, so = TRUE)
412
413
414###################################################
415### code chunk number 55: desTab
416###################################################
417rsdo = rsmChoose()
418
419
420###################################################
421### code chunk number 56: qualityTools.rnw:838-839
422###################################################
423rsdo = rsmChoose()
424
425
426###################################################
427### code chunk number 57: qualityTools.rnw:849-851
428###################################################
429fdo3 = facDesign(k = 6)
430rsdo = starDesign(alpha = "orthogonal", data = fdo3)
431
432
433###################################################
434### code chunk number 58: qualityTools.rnw:860-861
435###################################################
436randomize(fdo, random.seed = 123)
437
438
439###################################################
440### code chunk number 59: qualityTools.rnw:865-866
441###################################################
442randomize(fdo, so = TRUE)
443
444
445###################################################
446### code chunk number 60: qualityTools.rnw:903-906
447###################################################
448d1 = desirability(y1, 120, 170, scale = c(1,1), target = "max")
449d3 = desirability(y3, 400, 600, target = 500)
450d1
451
452
453###################################################
454### code chunk number 61: desirb
455###################################################
456par(mfrow = c(1,2))
457plot(d1, col = 2); plot(d3, col = 2)
458
459
460###################################################
461### code chunk number 62: qualityTools.rnw:917-918
462###################################################
463par(mfrow = c(1,2))
464plot(d1, col = 2); plot(d3, col = 2)
465
466
467###################################################
468### code chunk number 63: qualityTools.rnw:931-939
469###################################################
470ddo = rsmDesign(k = 3, alpha = 1.633, cc = 0, cs = 6)
471ddo = randomize(ddo, so = TRUE)
472#optional
473names(ddo) = c("silica", "silan", "sulfur")
474#optional
475highs(ddo) = c(1.7, 60, 2.8)
476#optional
477lows(ddo) = c(0.7, 40, 1.8)
478
479
480###################################################
481### code chunk number 64: qualityTools.rnw:944-952
482###################################################
483y1 = c(102, 120, 117, 198, 103, 132, 132, 139, 102, 154, 96, 163, 116,
484 153, 133, 133, 140, 142, 145, 142)
485y2 = c(900, 860, 800, 2294, 490, 1289, 1270, 1090, 770, 1690, 700, 1540,
486 2184, 1784, 1300, 1300, 1145, 1090, 1260, 1344)
487y3 = c(470, 410, 570, 240, 640, 270, 410, 380, 590, 260, 520, 380, 520,
488 290, 380, 380, 430, 430, 390, 390)
489y4 = c(67.5, 65, 77.5, 74.5, 62.5, 67, 78, 70, 76, 70, 63, 75, 65, 71,
490 70, 68.5, 68, 68, 69, 70)
491
492
493###################################################
494### code chunk number 65: qualityTools.rnw:957-958
495###################################################
496response(ddo) = data.frame(y1, y2, y3, y4)[c(5,2,3,8,1,6,7,4,9:20),]
497
498
499###################################################
500### code chunk number 66: qualityTools.rnw:963-965
501###################################################
502d2 = desirability(y2, 1000, 1300, target = "max")
503d4 = desirability(y4, 60, 75, target = 67.5)
504
505
506###################################################
507### code chunk number 67: qualityTools.rnw:970-971
508###################################################
509desires(ddo)=d1; desires(ddo)=d2; desires(ddo)=d3; desires(ddo)=d4
510
511
512###################################################
513### code chunk number 68: qualityTools.rnw:976-980
514###################################################
515fits(ddo) = lm(y1 ~ A+B+C+A:B+A:C+B:C+I(A^2)+I(B^2)+I(C^2), data = ddo)
516fits(ddo) = lm(y2 ~ A+B+C+A:B+A:C+B:C+I(A^2)+I(B^2)+I(C^2), data = ddo)
517fits(ddo) = lm(y3 ~ A+B+C+A:B+A:C+B:C+I(A^2)+I(B^2)+I(C^2), data = ddo)
518fits(ddo) = lm(y4 ~ A+B+C+A:B+A:C+B:C+I(A^2)+I(B^2)+I(C^2), data = ddo)
519
520
521###################################################
522### code chunk number 69: qualityTools.rnw:985-986
523###################################################
524optimum(ddo, type = "optim")
525
526
527###################################################
528### code chunk number 70: qualityTools.rnw:995-1003
529###################################################
530mdo = mixDesign(3,2, center = FALSE, axial = FALSE, randomize = FALSE,
531replicates  = c(1,1,2,3))
532names(mdo) = c("polyethylene", "polystyrene", "polypropylene")
533
534#set response (i.e. yarn elongation)
535elongation = c(11.0, 12.4, 15.0, 14.8, 16.1, 17.7, 16.4, 16.6, 8.8, 10.0,
53610.0, 9.7, 11.8, 16.8, 16.0)
537response(mdo) = elongation
538
539
540###################################################
541### code chunk number 71: qualityTools.rnw:1008-1009
542###################################################
543mdo
544
545
546###################################################
547### code chunk number 72: respAcont
548###################################################
549par(mfrow=c(1,2))
550contourPlot3(A, B, C, elongation, data = mdo, form = "quadratic")
551wirePlot3(A, B, C, elongation, data=mdo, form="quadratic", theta=-170)
552
553
554###################################################
555### code chunk number 73: qualityTools.rnw:1021-1022
556###################################################
557par(mfrow=c(1,2))
558contourPlot3(A, B, C, elongation, data = mdo, form = "quadratic")
559wirePlot3(A, B, C, elongation, data=mdo, form="quadratic", theta=-170)
560
561
562###################################################
563### code chunk number 74: qualityTools.rnw:1040-1046
564###################################################
565set.seed(1234)
566tdo = taguchiDesign("L9_3")
567values(tdo) = list(A  = c(20, 40, 60), B = c("mateial 1", "material 2",
568"material 3"), C = c(1,2,3))
569names(tdo) = c("Factor 1", "Factor 2", "Factor 3", "Factor 4")
570summary(tdo)
571
572
573###################################################
574### code chunk number 75: taguchi
575###################################################
576response(tdo) = rnorm(9)
577effectPlot(tdo, col = 2)
578
579
580###################################################
581### code chunk number 76: qualityTools.rnw:1057-1058
582###################################################
583response(tdo) = rnorm(9)
584effectPlot(tdo, col = 2)
585
586
587###################################################
588### code chunk number 77: qualityTools.rnw:1095-1096
589###################################################
590sessionInfo()
591
592
593