Module 4 Solved Practice Problems
Exam Practice - Work through each problem by hand first using scratch paper and a scientific calculator. Click "View Step-by-Step Solution" to verify your arithmetic, ANOVA partition tables, and final statistical decisions.
Problem 1: Completely Randomized Design (CRD)
Problem Statement
Three fertilizer compounds (A , B , C A, B, C A , B , C ) are applied to 12 homogeneous experimental plots (t = 3 , r = 4 , N = 12 t = 3, r = 4, N = 12 t = 3 , r = 4 , N = 12 ). The measured crop yields (in kg) are:
Fertilizer A: 18 , 20 , 19 , 22 18, 20, 19, 22 18 , 20 , 19 , 22
Fertilizer B: 24 , 22 , 21 , 23 24, 22, 21, 23 24 , 22 , 21 , 23
Fertilizer C: 26 , 25 , 27 , 28 26, 25, 27, 28 26 , 25 , 27 , 28
Test at the α = 0.05 \alpha = 0.05 α = 0.05 significance level whether there is a statistically significant difference in mean yield among the three fertilizers. (Given: Critical F 0.05 ( 2 , 9 ) ≈ 4.26 F_{0.05}(2, 9) \approx 4.26 F 0.05 ( 2 , 9 ) ≈ 4.26 ).
View Step-by-Step Solution
Null Hypothesis (H 0 H_0 H 0 ): μ A = μ B = μ C \mu_A = \mu_B = \mu_C μ A = μ B = μ C (The mean yields of all three fertilizers are equal).
Alternative Hypothesis (H 1 H_1 H 1 ): At least one fertilizer mean yield is different.
2. Totals & Correction Factor (C F CF C F )
T A = 18 + 20 + 19 + 22 = 79 T_A = 18 + 20 + 19 + 22 = 79 T A = 18 + 20 + 19 + 22 = 79
T B = 24 + 22 + 21 + 23 = 90 T_B = 24 + 22 + 21 + 23 = 90 T B = 24 + 22 + 21 + 23 = 90
T C = 26 + 25 + 27 + 28 = 106 T_C = 26 + 25 + 27 + 28 = 106 T C = 26 + 25 + 27 + 28 = 106
Grand Total G = 79 + 90 + 106 = 275 \text{Grand Total } G = 79 + 90 + 106 = 275 Grand Total G = 79 + 90 + 106 = 275
C F = G 2 N = 275 2 12 = 75625 12 ≈ 6302.0833 CF = \frac{G^2}{N} = \frac{275^2}{12} = \frac{75625}{12} \approx 6302.0833 C F = N G 2 = 12 27 5 2 = 12 75625 ≈ 6302.0833
3. Sum of Squares Computations
∑ y i j 2 = 18 2 + 20 2 + 19 2 + 22 2 + 24 2 + 22 2 + 21 2 + 23 2 + 26 2 + 25 2 + 27 2 + 28 2 = 6413 \sum y_{ij}^2 = 18^2 + 20^2 + 19^2 + 22^2 + 24^2 + 22^2 + 21^2 + 23^2 + 26^2 + 25^2 + 27^2 + 28^2 = 6413 ∑ y ij 2 = 1 8 2 + 2 0 2 + 1 9 2 + 2 2 2 + 2 4 2 + 2 2 2 + 2 1 2 + 2 3 2 + 2 6 2 + 2 5 2 + 2 7 2 + 2 8 2 = 6413
SST = ∑ y i j 2 − C F = 6413 − 6302.0833 = 110.9167 \text{SST} = \sum y_{ij}^2 - CF = 6413 - 6302.0833 = 110.9167 SST = ∑ y ij 2 − C F = 6413 − 6302.0833 = 110.9167
SSTreat = 79 2 + 90 2 + 106 2 4 − C F = 6241 + 8100 + 11236 4 − 6302.0833 = 6394.25 − 6302.0833 = 92.1667 \text{SSTreat} = \frac{79^2 + 90^2 + 106^2}{4} - CF = \frac{6241 + 8100 + 11236}{4} - 6302.0833 = 6394.25 - 6302.0833 = 92.1667 SSTreat = 4 7 9 2 + 9 0 2 + 10 6 2 − C F = 4 6241 + 8100 + 11236 − 6302.0833 = 6394.25 − 6302.0833 = 92.1667
SSE = SST − SSTreat = 110.9167 − 92.1667 = 18.7500 \text{SSE} = \text{SST} - \text{SSTreat} = 110.9167 - 92.1667 = 18.7500 SSE = SST − SSTreat = 110.9167 − 92.1667 = 18.7500
4. CRD ANOVA Table Source df \text{df} df SS \text{SS} SS MS \text{MS} MS Calculated F F F Tabulated F 0.05 F_{0.05} F 0.05 Fertilizers (Treatments) 3 − 1 = 2 3 - 1 = 2 3 − 1 = 2 92.1667 92.1667 92.1667 92.1667 2 = 46.0833 \frac{92.1667}{2} = 46.0833 2 92.1667 = 46.0833 46.0833 2.0833 ≈ 22.12 \frac{46.0833}{2.0833} \approx 22.12 2.0833 46.0833 ≈ 22.12 4.26 4.26 4.26 Error (Residual) 12 − 3 = 9 12 - 3 = 9 12 − 3 = 9 18.7500 18.7500 18.7500 18.7500 9 = 2.0833 \frac{18.7500}{9} = 2.0833 9 18.7500 = 2.0833 - - Total 12 − 1 = 11 12 - 1 = 11 12 − 1 = 11 110.9167 110.9167 110.9167 - - -
5. Decision & Academic Conclusion
Comparison: F cal = 22.12 > F tab = 4.26 F_{\text{cal}} = 22.12 > F_{\text{tab}} = 4.26 F cal = 22.12 > F tab = 4.26 .
Decision: Reject H 0 H_0 H 0 at the 5 % 5\% 5% level of significance.
Conclusion: "Since the calculated F F F -value (22.12 22.12 22.12 ) is greater than the tabulated F F F -value (4.26 4.26 4.26 ) at the 5 % 5\% 5% level of significance, the null hypothesis (H 0 H_0 H 0 ) is rejected, indicating that there is a statistically significant difference in mean yield among the three fertilizers."
Problem 2: Randomized Block Design (RBD)
Problem Statement
An experiment measures cholesterol content (%) across 4 different diet foods (D 1 , D 2 , D 3 , D 4 D_1, D_2, D_3, D_4 D 1 , D 2 , D 3 , D 4 ) tested across 4 distinct laboratories (A , B , C , D A, B, C, D A , B , C , D ) to control for inter-laboratory testing bias (t = 4 , r = 4 , N = 16 t = 4, r = 4, N = 16 t = 4 , r = 4 , N = 16 ).
Laboratory (Block) D 1 D_1 D 1 D 2 D_2 D 2 D 3 D_3 D 3 D 4 D_4 D 4 A 2 5 4 3 B 3 7 4 2 C 5 5 6 4 D 5 5 8 5
Analyze the data using RBD and state your conclusion at α = 0.05 \alpha = 0.05 α = 0.05 . (Given: Critical F 0.05 ( 3 , 9 ) = 3.86 F_{0.05}(3, 9) = 3.86 F 0.05 ( 3 , 9 ) = 3.86 ).
View Step-by-Step Solution 1. Totals Calculation
Treatment Totals (T i T_i T i ):
T D 1 = 2 + 3 + 5 + 5 = 15 T_{D1} = 2 + 3 + 5 + 5 = 15 T D 1 = 2 + 3 + 5 + 5 = 15
T D 2 = 5 + 7 + 5 + 5 = 22 T_{D2} = 5 + 7 + 5 + 5 = 22 T D 2 = 5 + 7 + 5 + 5 = 22
T D 3 = 4 + 4 + 6 + 8 = 22 T_{D3} = 4 + 4 + 6 + 8 = 22 T D 3 = 4 + 4 + 6 + 8 = 22
T D 4 = 3 + 2 + 4 + 5 = 14 T_{D4} = 3 + 2 + 4 + 5 = 14 T D 4 = 3 + 2 + 4 + 5 = 14
Block Totals (B j B_j B j ):
B A = 2 + 5 + 4 + 3 = 14 B_A = 2 + 5 + 4 + 3 = 14 B A = 2 + 5 + 4 + 3 = 14
B B = 3 + 7 + 4 + 2 = 16 B_B = 3 + 7 + 4 + 2 = 16 B B = 3 + 7 + 4 + 2 = 16
B C = 5 + 5 + 6 + 4 = 20 B_C = 5 + 5 + 6 + 4 = 20 B C = 5 + 5 + 6 + 4 = 20
B D = 5 + 5 + 8 + 5 = 23 B_D = 5 + 5 + 8 + 5 = 23 B D = 5 + 5 + 8 + 5 = 23
Grand Total G = 15 + 22 + 22 + 14 = 73 \text{Grand Total } G = 15 + 22 + 22 + 14 = 73 Grand Total G = 15 + 22 + 22 + 14 = 73
C F = G 2 N = 73 2 16 = 5329 16 = 333.0625 CF = \frac{G^2}{N} = \frac{73^2}{16} = \frac{5329}{16} = 333.0625 C F = N G 2 = 16 7 3 2 = 16 5329 = 333.0625
2. Sum of Squares Partitioning
∑ y i j 2 = 4 + 25 + 16 + 9 + 9 + 49 + 16 + 4 + 25 + 25 + 36 + 16 + 25 + 25 + 64 + 25 = 373 \sum y_{ij}^2 = 4 + 25 + 16 + 9 + 9 + 49 + 16 + 4 + 25 + 25 + 36 + 16 + 25 + 25 + 64 + 25 = 373 ∑ y ij 2 = 4 + 25 + 16 + 9 + 9 + 49 + 16 + 4 + 25 + 25 + 36 + 16 + 25 + 25 + 64 + 25 = 373
SST = 373 − 333.0625 = 39.9375 \text{SST} = 373 - 333.0625 = 39.9375 SST = 373 − 333.0625 = 39.9375
SSTreat = 15 2 + 22 2 + 22 2 + 14 2 4 − C F = 1389 4 − 333.0625 = 347.25 − 333.0625 = 14.1875 \text{SSTreat} = \frac{15^2 + 22^2 + 22^2 + 14^2}{4} - CF = \frac{1389}{4} - 333.0625 = 347.25 - 333.0625 = 14.1875 SSTreat = 4 1 5 2 + 2 2 2 + 2 2 2 + 1 4 2 − C F = 4 1389 − 333.0625 = 347.25 − 333.0625 = 14.1875
SSBlock = 14 2 + 16 2 + 20 2 + 23 2 4 − C F = 1381 4 − 333.0625 = 345.25 − 333.0625 = 12.1875 \text{SSBlock} = \frac{14^2 + 16^2 + 20^2 + 23^2}{4} - CF = \frac{1381}{4} - 333.0625 = 345.25 - 333.0625 = 12.1875 SSBlock = 4 1 4 2 + 1 6 2 + 2 0 2 + 2 3 2 − C F = 4 1381 − 333.0625 = 345.25 − 333.0625 = 12.1875
SSE = SST − SSTreat − SSBlock = 39.9375 − 14.1875 − 12.1875 = 13.5625 \text{SSE} = \text{SST} - \text{SSTreat} - \text{SSBlock} = 39.9375 - 14.1875 - 12.1875 = 13.5625 SSE = SST − SSTreat − SSBlock = 39.9375 − 14.1875 − 12.1875 = 13.5625
3. RBD ANOVA Table Source df \text{df} df SS \text{SS} SS MS \text{MS} MS Calculated F F F Critical F 0.05 F_{0.05} F 0.05 Diets (Treatments) 4 − 1 = 3 4 - 1 = 3 4 − 1 = 3 14.1875 14.1875 14.1875 14.1875 3 ≈ 4.7292 \frac{14.1875}{3} \approx 4.7292 3 14.1875 ≈ 4.7292 4.7292 1.5069 ≈ 3.14 \frac{4.7292}{1.5069} \approx 3.14 1.5069 4.7292 ≈ 3.14 3.86 3.86 3.86 Labs (Blocks) 4 − 1 = 3 4 - 1 = 3 4 − 1 = 3 12.1875 12.1875 12.1875 12.1875 3 ≈ 4.0625 \frac{12.1875}{3} \approx 4.0625 3 12.1875 ≈ 4.0625 - - Error ( 4 − 1 ) ( 4 − 1 ) = 9 (4-1)(4-1)=9 ( 4 − 1 ) ( 4 − 1 ) = 9 13.5625 13.5625 13.5625 13.5625 9 ≈ 1.5069 \frac{13.5625}{9} \approx 1.5069 9 13.5625 ≈ 1.5069 - - Total 16 − 1 = 15 16 - 1 = 15 16 − 1 = 15 39.9375 39.9375 39.9375 - - -
4. Conclusion
F cal = 3.14 < F tab = 3.86 F_{\text{cal}} = 3.14 < F_{\text{tab}} = 3.86 F cal = 3.14 < F tab = 3.86 .
Fail to reject H 0 H_0 H 0 . At the 5 % 5\% 5% level of significance, there is no statistically significant difference in cholesterol content among the four diet foods after controlling for laboratory variation.
Problem 3: Latin Square Design (LSD)
Problem Statement
Four treatments (A , B , C , D A, B, C, D A , B , C , D ) are allocated in a 4 × 4 4 \times 4 4 × 4 Latin square to control for two independent nuisance sources (e.g., Row = Machine, Column = Operator). The experimental layout and observations are:
Row 1: A ( 18 ) , B ( 22 ) , C ( 25 ) , D ( 28 ) A(18), B(22), C(25), D(28) A ( 18 ) , B ( 22 ) , C ( 25 ) , D ( 28 )
Row 2: B ( 21 ) , C ( 24 ) , D ( 26 ) , A ( 20 ) B(21), C(24), D(26), A(20) B ( 21 ) , C ( 24 ) , D ( 26 ) , A ( 20 )
Row 3: C ( 19 ) , D ( 23 ) , A ( 22 ) , B ( 25 ) C(19), D(23), A(22), B(25) C ( 19 ) , D ( 23 ) , A ( 22 ) , B ( 25 )
Row 4: D ( 20 ) , A ( 21 ) , B ( 23 ) , C ( 27 ) D(20), A(21), B(23), C(27) D ( 20 ) , A ( 21 ) , B ( 23 ) , C ( 27 )
Test at α = 0.05 \alpha = 0.05 α = 0.05 whether treatment differences are significant. (Given: Critical F 0.05 ( 3 , 6 ) = 4.76 F_{0.05}(3, 6) = 4.76 F 0.05 ( 3 , 6 ) = 4.76 , ∑ y 2 = 8408 \sum y^2 = 8408 ∑ y 2 = 8408 ).
View Step-by-Step Solution 1. Totals Calculation (p = 4 , N = 16 p = 4, N = 16 p = 4 , N = 16 )
Treatment Totals:
T A = 18 + 20 + 22 + 21 = 81 T_A = 18 + 20 + 22 + 21 = 81 T A = 18 + 20 + 22 + 21 = 81
T B = 22 + 21 + 25 + 23 = 91 T_B = 22 + 21 + 25 + 23 = 91 T B = 22 + 21 + 25 + 23 = 91
T C = 25 + 24 + 19 + 27 = 95 T_C = 25 + 24 + 19 + 27 = 95 T C = 25 + 24 + 19 + 27 = 95
T D = 28 + 26 + 23 + 20 = 97 T_D = 28 + 26 + 23 + 20 = 97 T D = 28 + 26 + 23 + 20 = 97
Row Totals:
R 1 = 18 + 22 + 25 + 28 = 93 R_1 = 18 + 22 + 25 + 28 = 93 R 1 = 18 + 22 + 25 + 28 = 93
R 2 = 21 + 24 + 26 + 20 = 91 R_2 = 21 + 24 + 26 + 20 = 91 R 2 = 21 + 24 + 26 + 20 = 91
R 3 = 19 + 23 + 22 + 25 = 89 R_3 = 19 + 23 + 22 + 25 = 89 R 3 = 19 + 23 + 22 + 25 = 89
R 4 = 20 + 21 + 23 + 27 = 91 R_4 = 20 + 21 + 23 + 27 = 91 R 4 = 20 + 21 + 23 + 27 = 91
Column Totals:
C 1 = 18 + 21 + 19 + 20 = 78 C_1 = 18 + 21 + 19 + 20 = 78 C 1 = 18 + 21 + 19 + 20 = 78
C 2 = 22 + 24 + 23 + 21 = 90 C_2 = 22 + 24 + 23 + 21 = 90 C 2 = 22 + 24 + 23 + 21 = 90
C 3 = 25 + 26 + 22 + 23 = 96 C_3 = 25 + 26 + 22 + 23 = 96 C 3 = 25 + 26 + 22 + 23 = 96
C 4 = 28 + 20 + 25 + 27 = 100 C_4 = 28 + 20 + 25 + 27 = 100 C 4 = 28 + 20 + 25 + 27 = 100
G = 364 G = 364 G = 364
C F = 364 2 16 = 8281.0 CF = \frac{364^2}{16} = 8281.0 C F = 16 36 4 2 = 8281.0
2. Sum of Squares Partitioning
SST = 8408.0 − 8281.0 = 127.0 \text{SST} = 8408.0 - 8281.0 = 127.0 SST = 8408.0 − 8281.0 = 127.0
SSTreat = 81 2 + 91 2 + 95 2 + 97 2 4 − 8281.0 = 33276 4 − 8281.0 = 8319.0 − 8281.0 = 38.0 \text{SSTreat} = \frac{81^2 + 91^2 + 95^2 + 97^2}{4} - 8281.0 = \frac{33276}{4} - 8281.0 = 8319.0 - 8281.0 = 38.0 SSTreat = 4 8 1 2 + 9 1 2 + 9 5 2 + 9 7 2 − 8281.0 = 4 33276 − 8281.0 = 8319.0 − 8281.0 = 38.0
SSRow = 93 2 + 91 2 + 89 2 + 91 2 4 − 8281.0 = 33132 4 − 8281.0 = 8283.0 − 8281.0 = 2.0 \text{SSRow} = \frac{93^2 + 91^2 + 89^2 + 91^2}{4} - 8281.0 = \frac{33132}{4} - 8281.0 = 8283.0 - 8281.0 = 2.0 SSRow = 4 9 3 2 + 9 1 2 + 8 9 2 + 9 1 2 − 8281.0 = 4 33132 − 8281.0 = 8283.0 − 8281.0 = 2.0
SSCol = 78 2 + 90 2 + 96 2 + 100 2 4 − 8281.0 = 33400 4 − 8281.0 = 8350.0 − 8281.0 = 69.0 \text{SSCol} = \frac{78^2 + 90^2 + 96^2 + 100^2}{4} - 8281.0 = \frac{33400}{4} - 8281.0 = 8350.0 - 8281.0 = 69.0 SSCol = 4 7 8 2 + 9 0 2 + 9 6 2 + 10 0 2 − 8281.0 = 4 33400 − 8281.0 = 8350.0 − 8281.0 = 69.0
SSE = 127.0 − 38.0 − 2.0 − 69.0 = 18.0 \text{SSE} = 127.0 - 38.0 - 2.0 - 69.0 = 18.0 SSE = 127.0 − 38.0 − 2.0 − 69.0 = 18.0
3. LSD ANOVA Table Source df \text{df} df SS \text{SS} SS MS \text{MS} MS Calculated F F F Critical F 0.05 F_{0.05} F 0.05 Rows 4 − 1 = 3 4 - 1 = 3 4 − 1 = 3 2.0 2.0 2.0 0.6667 0.6667 0.6667 - - Columns 4 − 1 = 3 4 - 1 = 3 4 − 1 = 3 69.0 69.0 69.0 23.0000 23.0000 23.0000 - - Treatments 4 − 1 = 3 4 - 1 = 3 4 − 1 = 3 38.0 38.0 38.0 12.6667 12.6667 12.6667 12.6667 3.0000 ≈ 4.222 \frac{12.6667}{3.0000} \approx 4.222 3.0000 12.6667 ≈ 4.222 4.76 4.76 4.76 Error ( 4 − 1 ) ( 4 − 2 ) = 6 (4-1)(4-2)=6 ( 4 − 1 ) ( 4 − 2 ) = 6 18.0 18.0 18.0 3.0000 3.0000 3.0000 - - Total 16 − 1 = 15 16 - 1 = 15 16 − 1 = 15 127.0 127.0 127.0 - - -
4. Conclusion
F cal = 4.222 < F tab = 4.76 F_{\text{cal}} = 4.222 < F_{\text{tab}} = 4.76 F cal = 4.222 < F tab = 4.76 .
Fail to reject H 0 H_0 H 0 . After adjusting for row and column effects, there is no statistically significant treatment effect at the 5 % 5\% 5% level.
Problem 4: 2 2 2^2 2 2 Full Factorial Design
Problem Statement
A chemical engineer investigates the effect of Temperature (A A A ) and Pressure (B B B ) on product yield (%). Both factors are set at two coded levels: Low (− 1 -1 − 1 ) and High (+ 1 +1 + 1 ) with r = 2 r=2 r = 2 replicates (N = 2 2 × 2 = 8 N = 2^2 \times 2 = 8 N = 2 2 × 2 = 8 ).
Run Factor A (Temp) Factor B (Pressure) Replicate 1 Replicate 2 1 − 1 -1 − 1 − 1 -1 − 1 42 43 2 + 1 +1 + 1 − 1 -1 − 1 50 49 3 − 1 -1 − 1 + 1 +1 + 1 46 45 4 + 1 +1 + 1 + 1 +1 + 1 62 61
Perform a complete 2 2 2^2 2 2 factorial analysis: compute contrasts, main effects, interaction effect, construct the ANOVA table, and test for significance at α = 0.05 \alpha = 0.05 α = 0.05 . (Given: Critical F 0.05 ( 1 , 4 ) = 7.708 F_{0.05}(1, 4) = 7.708 F 0.05 ( 1 , 4 ) = 7.708 ).
View Step-by-Step Solution 1. Treatment Totals
Run 1 (− 1 , − 1 -1, -1 − 1 , − 1 ): y 1 = 42 + 43 = 85 y_1 = 42 + 43 = 85 y 1 = 42 + 43 = 85
Run 2 (+ 1 , − 1 +1, -1 + 1 , − 1 ): y 2 = 50 + 49 = 99 y_2 = 50 + 49 = 99 y 2 = 50 + 49 = 99
Run 3 (− 1 , + 1 -1, +1 − 1 , + 1 ): y 3 = 46 + 45 = 91 y_3 = 46 + 45 = 91 y 3 = 46 + 45 = 91
Run 4 (+ 1 , + 1 +1, +1 + 1 , + 1 ): y 4 = 62 + 61 = 123 y_4 = 62 + 61 = 123 y 4 = 62 + 61 = 123
Grand Total G = 85 + 99 + 91 + 123 = 398 \text{Grand Total } G = 85 + 99 + 91 + 123 = 398 Grand Total G = 85 + 99 + 91 + 123 = 398
2. Contrasts & Effects (2 k − 1 ⋅ r = 2 2 − 1 ⋅ 2 = 4 2^{k-1} \cdot r = 2^{2-1} \cdot 2 = 4 2 k − 1 ⋅ r = 2 2 − 1 ⋅ 2 = 4 )
C A = − 85 + 99 − 91 + 123 = 46 ⟹ A ^ = 46 4 = 11.5 C_A = -85 + 99 - 91 + 123 = 46 \implies \hat{A} = \frac{46}{4} = 11.5 C A = − 85 + 99 − 91 + 123 = 46 ⟹ A ^ = 4 46 = 11.5
C B = − 85 − 99 + 91 + 123 = 30 ⟹ B ^ = 30 4 = 7.5 C_B = -85 - 99 + 91 + 123 = 30 \implies \hat{B} = \frac{30}{4} = 7.5 C B = − 85 − 99 + 91 + 123 = 30 ⟹ B ^ = 4 30 = 7.5
C A B = + 85 − 99 − 91 + 123 = 18 ⟹ A B ^ = 18 4 = 4.5 C_{AB} = +85 - 99 - 91 + 123 = 18 \implies \widehat{AB} = \frac{18}{4} = 4.5 C A B = + 85 − 99 − 91 + 123 = 18 ⟹ A B = 4 18 = 4.5
3. Sum of Squares (2 k ⋅ r = 8 2^k \cdot r = 8 2 k ⋅ r = 8 )
SS A = 46 2 8 = 264.5 \text{SS}_A = \frac{46^2}{8} = 264.5 SS A = 8 4 6 2 = 264.5
SS B = 30 2 8 = 112.5 \text{SS}_B = \frac{30^2}{8} = 112.5 SS B = 8 3 0 2 = 112.5
SS A B = 18 2 8 = 40.5 \text{SS}_{AB} = \frac{18^2}{8} = 40.5 SS A B = 8 1 8 2 = 40.5
∑ y 2 = 42 2 + 43 2 + ⋯ + 61 2 = 20220 \sum y^2 = 42^2 + 43^2 + \dots + 61^2 = 20220 ∑ y 2 = 4 2 2 + 4 3 2 + ⋯ + 6 1 2 = 20220
C F = 398 2 8 = 19800.5 ⟹ SST = 20220 − 19800.5 = 419.5 CF = \frac{398^2}{8} = 19800.5 \implies \text{SST} = 20220 - 19800.5 = 419.5 C F = 8 39 8 2 = 19800.5 ⟹ SST = 20220 − 19800.5 = 419.5
SSE = SST − ( SS A + SS B + SS A B ) = 419.5 − ( 264.5 + 112.5 + 40.5 ) = 2.0 \text{SSE} = \text{SST} - (\text{SS}_A + \text{SS}_B + \text{SS}_{AB}) = 419.5 - (264.5 + 112.5 + 40.5) = 2.0 SSE = SST − ( SS A + SS B + SS A B ) = 419.5 − ( 264.5 + 112.5 + 40.5 ) = 2.0
4. Factorial ANOVA Table Source df \text{df} df SS \text{SS} SS MS \text{MS} MS Calculated F F F Critical F 0.05 F_{0.05} F 0.05 Decision A (Temperature) 1 264.5 264.5 264.5 264.5 264.5 264.5 264.5 0.5 = 529.0 \frac{264.5}{0.5} = 529.0 0.5 264.5 = 529.0 7.708 7.708 7.708 Reject H 0 H_0 H 0 B (Pressure) 1 112.5 112.5 112.5 112.5 112.5 112.5 112.5 0.5 = 225.0 \frac{112.5}{0.5} = 225.0 0.5 112.5 = 225.0 7.708 7.708 7.708 Reject H 0 H_0 H 0 AB (Interaction) 1 40.5 40.5 40.5 40.5 40.5 40.5 40.5 0.5 = 81.0 \frac{40.5}{0.5} = 81.0 0.5 40.5 = 81.0 7.708 7.708 7.708 Reject H 0 H_0 H 0 Error (Residual) 4 2.0 2.0 2.0 MSE = 0.5 \text{MSE} = 0.5 MSE = 0.5 - - - Total 7 419.5 419.5 419.5 - - - -
5. Conclusion All factors (A , B , A, B, A , B , and interaction A B AB A B ) are statistically significant at p < 0.001 p < 0.001 p < 0.001 . Temperature increases yield by 11.5 % 11.5\% 11.5% , Pressure increases yield by 7.5 % 7.5\% 7.5% , and positive interaction indicates the benefits of temperature are greatest when paired with high pressure.
Problem 5: Taguchi L 9 L_9 L 9 Orthogonal Array Optimization
Problem Statement
A manufacturing engineer seeks to optimize three 3-level parameters to maximize tensile strength (y y y in MPa \text{MPa} MPa ):
Factor A (Temperature): Level 1 = 200 ∘ C 200^\circ\text{C} 20 0 ∘ C , Level 2 = 220 ∘ C 220^\circ\text{C} 22 0 ∘ C , Level 3 = 240 ∘ C 240^\circ\text{C} 24 0 ∘ C
Factor B (Holding Time): Level 1 = 10 min 10\,\text{min} 10 min , Level 2 = 20 min 20\,\text{min} 20 min , Level 3 = 30 min 30\,\text{min} 30 min
Factor C (Pressure): Level 1 = 5 bar 5\,\text{bar} 5 bar , Level 2 = 7 bar 7\,\text{bar} 7 bar , Level 3 = 9 bar 9\,\text{bar} 9 bar
The measured single-run responses in an L 9 L_9 L 9 array are:
Run 1: (1, 1, 1) ⟹ y = 42 MPa \implies y = 42\,\text{MPa} ⟹ y = 42 MPa
Run 2: (1, 2, 2) ⟹ y = 48 MPa \implies y = 48\,\text{MPa} ⟹ y = 48 MPa
Run 3: (1, 3, 3) ⟹ y = 50 MPa \implies y = 50\,\text{MPa} ⟹ y = 50 MPa
Run 4: (2, 1, 2) ⟹ y = 55 MPa \implies y = 55\,\text{MPa} ⟹ y = 55 MPa
Run 5: (2, 2, 3) ⟹ y = 60 MPa \implies y = 60\,\text{MPa} ⟹ y = 60 MPa
Run 6: (2, 3, 1) ⟹ y = 58 MPa \implies y = 58\,\text{MPa} ⟹ y = 58 MPa
Run 7: (3, 1, 3) ⟹ y = 65 MPa \implies y = 65\,\text{MPa} ⟹ y = 65 MPa
Run 8: (3, 2, 1) ⟹ y = 62 MPa \implies y = 62\,\text{MPa} ⟹ y = 62 MPa
Run 9: (3, 3, 2) ⟹ y = 59 MPa \implies y = 59\,\text{MPa} ⟹ y = 59 MPa
Using the Larger-the-better S / N S/N S / N ratio formulation:
Compute the S / N S/N S / N ratio for each trial.
Build the factor-level response table and rank the factors.
Determine the optimal parameter combination.
Predict the optimum tensile strength performance.
View Step-by-Step Solution 1. S/N Computation (S / N = 20 log 10 ( y ) S/N = 20 \log_{10}(y) S / N = 20 log 10 ( y ) )
Run 1 : 20 log 10 ( 42 ) = 32.4650 dB \text{Run 1}: 20 \log_{10}(42) = 32.4650\,\text{dB} Run 1 : 20 log 10 ( 42 ) = 32.4650 dB
Run 2 : 20 log 10 ( 48 ) = 33.6248 dB \text{Run 2}: 20 \log_{10}(48) = 33.6248\,\text{dB} Run 2 : 20 log 10 ( 48 ) = 33.6248 dB
Run 3 : 20 log 10 ( 50 ) = 33.9794 dB \text{Run 3}: 20 \log_{10}(50) = 33.9794\,\text{dB} Run 3 : 20 log 10 ( 50 ) = 33.9794 dB
Run 4 : 20 log 10 ( 55 ) = 34.8073 dB \text{Run 4}: 20 \log_{10}(55) = 34.8073\,\text{dB} Run 4 : 20 log 10 ( 55 ) = 34.8073 dB
Run 5 : 20 log 10 ( 60 ) = 35.5630 dB \text{Run 5}: 20 \log_{10}(60) = 35.5630\,\text{dB} Run 5 : 20 log 10 ( 60 ) = 35.5630 dB
Run 6 : 20 log 10 ( 58 ) = 35.2686 dB \text{Run 6}: 20 \log_{10}(58) = 35.2686\,\text{dB} Run 6 : 20 log 10 ( 58 ) = 35.2686 dB
Run 7 : 20 log 10 ( 65 ) = 36.2583 dB \text{Run 7}: 20 \log_{10}(65) = 36.2583\,\text{dB} Run 7 : 20 log 10 ( 65 ) = 36.2583 dB
Run 8 : 20 log 10 ( 62 ) = 35.8478 dB \text{Run 8}: 20 \log_{10}(62) = 35.8478\,\text{dB} Run 8 : 20 log 10 ( 62 ) = 35.8478 dB
Run 9 : 20 log 10 ( 59 ) = 35.4170 dB \text{Run 9}: 20 \log_{10}(59) = 35.4170\,\text{dB} Run 9 : 20 log 10 ( 59 ) = 35.4170 dB
Overall Average: S / N ‾ = 34.8035 dB \overline{S/N} = 34.8035\,\text{dB} S / N = 34.8035 dB
2. Factor-Level S/N Response Table
Factor A (Temp):
Level 1: 32.4650 + 33.6248 + 33.9794 3 = 33.3564 dB \frac{32.4650 + 33.6248 + 33.9794}{3} = 33.3564\,\text{dB} 3 32.4650 + 33.6248 + 33.9794 = 33.3564 dB
Level 2: 34.8073 + 35.5630 + 35.2686 3 = 35.2130 dB \frac{34.8073 + 35.5630 + 35.2686}{3} = 35.2130\,\text{dB} 3 34.8073 + 35.5630 + 35.2686 = 35.2130 dB
Level 3: 36.2583 + 35.8478 + 35.4170 3 = 35.8410 dB \frac{36.2583 + 35.8478 + 35.4170}{3} = 35.8410\,\text{dB} 3 36.2583 + 35.8478 + 35.4170 = 35.8410 dB
Δ A = 35.8410 − 33.3564 = 2.4846 dB \Delta_A = 35.8410 - 33.3564 = 2.4846\,\text{dB} Δ A = 35.8410 − 33.3564 = 2.4846 dB (Rank 1)
Factor B (Time):
Level 1: 32.4650 + 34.8073 + 36.2583 3 = 34.5102 dB \frac{32.4650 + 34.8073 + 36.2583}{3} = 34.5102\,\text{dB} 3 32.4650 + 34.8073 + 36.2583 = 34.5102 dB
Level 2: 33.6248 + 35.5630 + 35.8478 3 = 35.0119 dB \frac{33.6248 + 35.5630 + 35.8478}{3} = 35.0119\,\text{dB} 3 33.6248 + 35.5630 + 35.8478 = 35.0119 dB
Level 3: 33.9794 + 35.2686 + 35.4170 3 = 34.8883 dB \frac{33.9794 + 35.2686 + 35.4170}{3} = 34.8883\,\text{dB} 3 33.9794 + 35.2686 + 35.4170 = 34.8883 dB
Δ B = 35.0119 − 34.5102 = 0.5017 dB \Delta_B = 35.0119 - 34.5102 = 0.5017\,\text{dB} Δ B = 35.0119 − 34.5102 = 0.5017 dB (Rank 3)
Factor C (Pressure):
Level 1: 32.4650 + 35.2686 + 35.8478 3 = 34.5271 dB \frac{32.4650 + 35.2686 + 35.8478}{3} = 34.5271\,\text{dB} 3 32.4650 + 35.2686 + 35.8478 = 34.5271 dB
Level 2: 33.6248 + 34.8073 + 35.4170 3 = 34.6164 dB \frac{33.6248 + 34.8073 + 35.4170}{3} = 34.6164\,\text{dB} 3 33.6248 + 34.8073 + 35.4170 = 34.6164 dB
Level 3: 33.9794 + 35.5630 + 36.2583 3 = 35.2669 dB \frac{33.9794 + 35.5630 + 36.2583}{3} = 35.2669\,\text{dB} 3 33.9794 + 35.5630 + 36.2583 = 35.2669 dB
Δ C = 35.2669 − 34.5271 = 0.7398 dB \Delta_C = 35.2669 - 34.5271 = 0.7398\,\text{dB} Δ C = 35.2669 − 34.5271 = 0.7398 dB (Rank 2)
3. Optimal Setting Selection Always pick the level with the highest mean S / N S/N S / N :
Factor A: Level 3 (240 ∘ C 240^\circ\text{C} 24 0 ∘ C )
Factor B: Level 2 (20 min 20\,\text{min} 20 min )
Factor C: Level 3 (9 bar 9\,\text{bar} 9 bar )
Optimal Setting: A 3 B 2 C 3 A_3 B_2 C_3 A 3 B 2 C 3
S / N ‾ pred = 34.8035 + ( 35.8410 − 34.8035 ) + ( 35.0119 − 34.8035 ) + ( 35.2669 − 34.8035 ) = 36.5128 dB \overline{S/N}_{\text{pred}} = 34.8035 + (35.8410 - 34.8035) + (35.0119 - 34.8035) + (35.2669 - 34.8035) = 36.5128\,\text{dB} S / N pred = 34.8035 + ( 35.8410 − 34.8035 ) + ( 35.0119 − 34.8035 ) + ( 35.2669 − 34.8035 ) = 36.5128 dB
Predicted Tensile Strength: y ^ = 10 36.5128 20 ≈ 66.93 MPa \text{Predicted Tensile Strength: } \hat{y} = 10^{\frac{36.5128}{20}} \approx 66.93\,\text{MPa} Predicted Tensile Strength: y ^ = 1 0 20 36.5128 ≈ 66.93 MPa
Problem 6: Genetic Algorithm Optimization
Problem Statement
A Genetic Algorithm is used to maximize the function f ( x ) = x 2 f(x) = x^2 f ( x ) = x 2 over the integer domain 0 ≤ x ≤ 31 0 \le x \le 31 0 ≤ x ≤ 31 encoded as 5-bit binary strings. The initial random population consists of:
C 1 = 10110 2 C_1 = 10110_2 C 1 = 1011 0 2
C 2 = 00101 2 C_2 = 00101_2 C 2 = 0010 1 2
C 3 = 11100 2 C_3 = 11100_2 C 3 = 1110 0 2
C 4 = 01011 2 C_4 = 01011_2 C 4 = 0101 1 2
Perform one full generation of the algorithm:
Decode each chromosome to decimal and evaluate its fitness.
Select the top two chromosomes as parents.
Perform single-point crossover after the 2nd bit to produce two children.
Mutate the last bit of Child 2.
Form the next generation using the children and the top two elite parents from the original generation, and compute the new average population fitness.
View Step-by-Step Solution 1. Decoding & Fitness (f ( x ) = x 2 f(x) = x^2 f ( x ) = x 2 )
C 1 = 10110 2 = 16 + 4 + 2 = 22 ⟹ f ( 22 ) = 22 2 = 484 C_1 = 10110_2 = 16 + 4 + 2 = 22 \implies f(22) = 22^2 = 484 C 1 = 1011 0 2 = 16 + 4 + 2 = 22 ⟹ f ( 22 ) = 2 2 2 = 484
C 2 = 00101 2 = 4 + 1 = 5 ⟹ f ( 5 ) = 5 2 = 25 C_2 = 00101_2 = 4 + 1 = 5 \implies f(5) = 5^2 = 25 C 2 = 0010 1 2 = 4 + 1 = 5 ⟹ f ( 5 ) = 5 2 = 25
C 3 = 11100 2 = 16 + 8 + 4 = 28 ⟹ f ( 28 ) = 28 2 = 784 C_3 = 11100_2 = 16 + 8 + 4 = 28 \implies f(28) = 28^2 = 784 C 3 = 1110 0 2 = 16 + 8 + 4 = 28 ⟹ f ( 28 ) = 2 8 2 = 784
C 4 = 01011 2 = 8 + 2 + 1 = 11 ⟹ f ( 11 ) = 11 2 = 121 C_4 = 01011_2 = 8 + 2 + 1 = 11 \implies f(11) = 11^2 = 121 C 4 = 0101 1 2 = 8 + 2 + 1 = 11 ⟹ f ( 11 ) = 1 1 2 = 121
Gen 0 Average Fitness: f ˉ 0 = 484 + 25 + 784 + 121 4 = 353.5 \bar{f}_0 = \frac{484 + 25 + 784 + 121}{4} = 353.5 f ˉ 0 = 4 484 + 25 + 784 + 121 = 353.5
2. Parent Selection The two highest fitness individuals are:
Parent 1: C 3 = 11100 C_3 = 11100 C 3 = 11100 (x = 28 , f = 784 x = 28, f = 784 x = 28 , f = 784 )
Parent 2: C 1 = 10110 C_1 = 10110 C 1 = 10110 (x = 22 , f = 484 x = 22, f = 484 x = 22 , f = 484 )
3. Crossover (Cut after Bit 2)
Parent 1: [ 11 ∣ 100 ] \text{Parent 1: } [11 \mid 100] Parent 1: [ 11 ∣ 100 ]
Parent 2: [ 10 ∣ 110 ] \text{Parent 2: } [10 \mid 110] Parent 2: [ 10 ∣ 110 ]
Child 1: [ 11 ∣ 110 ] = 11110 2 = 30 ⟹ f ( 30 ) = 30 2 = 900 \text{Child 1: } [11 \mid 110] = 11110_2 = 30 \implies f(30) = 30^2 = 900 Child 1: [ 11 ∣ 110 ] = 1111 0 2 = 30 ⟹ f ( 30 ) = 3 0 2 = 900
Child 2: [ 10 ∣ 100 ] = 10100 2 = 20 ⟹ f ( 20 ) = 20 2 = 400 \text{Child 2: } [10 \mid 100] = 10100_2 = 20 \implies f(20) = 20^2 = 400 Child 2: [ 10 ∣ 100 ] = 1010 0 2 = 20 ⟹ f ( 20 ) = 2 0 2 = 400
4. Mutation Flip the 5th bit of Child 2 (10100 → 10101 10100 \to 10101 10100 → 10101 ):
Mutated Child 2 = 10101 2 = 21 ⟹ f ( 21 ) = 21 2 = 441 \text{Mutated Child 2} = 10101_2 = 21 \implies f(21) = 21^2 = 441 Mutated Child 2 = 1010 1 2 = 21 ⟹ f ( 21 ) = 2 1 2 = 441
Child 1: 11110 2 ⟹ x = 30 , f = 900 11110_2 \implies x = 30, f = 900 1111 0 2 ⟹ x = 30 , f = 900
Mutated Child 2: 10101 2 ⟹ x = 21 , f = 441 10101_2 \implies x = 21, f = 441 1010 1 2 ⟹ x = 21 , f = 441
Elite Parent 1: 11100 2 ⟹ x = 28 , f = 784 11100_2 \implies x = 28, f = 784 1110 0 2 ⟹ x = 28 , f = 784
Elite Parent 2: 10110 2 ⟹ x = 22 , f = 484 10110_2 \implies x = 22, f = 484 1011 0 2 ⟹ x = 22 , f = 484
Gen 1 Average Fitness:
f ˉ 1 = 900 + 441 + 784 + 484 4 = 2609 4 = 652.25 \bar{f}_1 = \frac{900 + 441 + 784 + 484}{4} = \frac{2609}{4} = 652.25 f ˉ 1 = 4 900 + 441 + 784 + 484 = 4 2609 = 652.25
Improvement: Average fitness increased by 652.25 − 353.5 353.5 × 100 % ≈ 84.5 % \frac{652.25 - 353.5}{353.5} \times 100\% \approx 84.5\% 353.5 652.25 − 353.5 × 100% ≈ 84.5% .
Problem 7: Response Surface Methodology (RSM) Optimization
Problem Statement
A chemical company wants to maximize product yield (Y Y Y , in % \% % ) . Two continuous factors are investigated:
Temperature (X 1 X_1 X 1 ): Low level (− 1 -1 − 1 ) = 100 ∘ C 100^\circ\text{C} 10 0 ∘ C , High level (+ 1 +1 + 1 ) = 140 ∘ C 140^\circ\text{C} 14 0 ∘ C (Center point 0 = 120 ∘ C 0 = 120^\circ\text{C} 0 = 12 0 ∘ C ).
Pressure (X 2 X_2 X 2 ): Low level (− 1 -1 − 1 ) = 10 psi 10\,\text{psi} 10 psi , High level (+ 1 +1 + 1 ) = 20 psi 20\,\text{psi} 20 psi (Center point 0 = 15 psi 0 = 15\,\text{psi} 0 = 15 psi ).
The experimental results from an initial 2 2 2^2 2 2 factorial design matrix are:
Run X 1 X_1 X 1 (Temperature)X 2 X_2 X 2 (Pressure)Yield Y Y Y (%) 1 − 1 -1 − 1 − 1 -1 − 1 70 2 − 1 -1 − 1 + 1 +1 + 1 76 3 + 1 +1 + 1 − 1 -1 − 1 80 4 + 1 +1 + 1 + 1 +1 + 1 86
Perform the following:
Assume a first-order model Y = β 0 + β 1 X 1 + β 2 X 2 Y = \beta_0 + \beta_1 X_1 + \beta_2 X_2 Y = β 0 + β 1 X 1 + β 2 X 2 and calculate the coefficients β 0 , β 1 , β 2 \beta_0, \beta_1, \beta_2 β 0 , β 1 , β 2 .
Write the fitted mathematical model.
Verify the model at run ( X 1 = + 1 , X 2 = + 1 ) (X_1 = +1, X_2 = +1) ( X 1 = + 1 , X 2 = + 1 ) .
Predict the yield at the center point ( X 1 = 0 , X 2 = 0 ) (X_1 = 0, X_2 = 0) ( X 1 = 0 , X 2 = 0 ) .
Determine the direction of improvement and identify which factor has the stronger effect.
Execute the Method of Steepest Ascent : determine the direction vector and list the coded and physical factor coordinates for the first 3 steps using step increments Δ X 1 = 0.5 \Delta X_1 = 0.5 Δ X 1 = 0.5 and Δ X 2 = 0.3 \Delta X_2 = 0.3 Δ X 2 = 0.3 .
Suppose subsequent trials exhibit curvature and a second-order model is fitted:
Y = 90 + 4 X 1 + 3 X 2 − 2 X 1 X 2 − 5 X 1 2 − 4 X 2 2 Y = 90 + 4X_1 + 3X_2 - 2X_1 X_2 - 5X_1^2 - 4X_2^2 Y = 90 + 4 X 1 + 3 X 2 − 2 X 1 X 2 − 5 X 1 2 − 4 X 2 2
If mathematical optimization predicts an optimum at ( X 1 = 0.4 , X 2 = 0.3 ) (X_1 = 0.4, X_2 = 0.3) ( X 1 = 0.4 , X 2 = 0.3 ) yielding 92 % 92\% 92% , validate the model if an actual confirmation trial produces 91.5 % 91.5\% 91.5% .
View Step-by-Step Solution 1. Calculate First-Order Regression Coefficients Using the orthogonal properties of the 2 2 2^2 2 2 coded matrix (N = 4 N = 4 N = 4 ):
Calculate β 0 \beta_0 β 0 (Intercept / Average Response):
β 0 = ∑ Y i N = 70 + 76 + 80 + 86 4 = 312 4 = 78 \beta_0 = \frac{\sum Y_i}{N} = \frac{70 + 76 + 80 + 86}{4} = \frac{312}{4} = 78 β 0 = N ∑ Y i = 4 70 + 76 + 80 + 86 = 4 312 = 78
β 0 = 78 \boxed{\beta_0 = 78} β 0 = 78
Calculate β 1 \beta_1 β 1 (Main Effect of Temperature X 1 X_1 X 1 ):
β 1 = ∑ X 1 i Y i N = ( − 1 ) ( 70 ) + ( − 1 ) ( 76 ) + ( + 1 ) ( 80 ) + ( + 1 ) ( 86 ) 4 = − 70 − 76 + 80 + 86 4 = 20 4 = 5 \beta_1 = \frac{\sum X_{1i} Y_i}{N} = \frac{(-1)(70) + (-1)(76) + (+1)(80) + (+1)(86)}{4} = \frac{-70 - 76 + 80 + 86}{4} = \frac{20}{4} = 5 β 1 = N ∑ X 1 i Y i = 4 ( − 1 ) ( 70 ) + ( − 1 ) ( 76 ) + ( + 1 ) ( 80 ) + ( + 1 ) ( 86 ) = 4 − 70 − 76 + 80 + 86 = 4 20 = 5
β 1 = 5 \boxed{\beta_1 = 5} β 1 = 5
Calculate β 2 \beta_2 β 2 (Main Effect of Pressure X 2 X_2 X 2 ):
β 2 = ∑ X 2 i Y i N = ( − 1 ) ( 70 ) + ( + 1 ) ( 76 ) + ( − 1 ) ( 80 ) + ( + 1 ) ( 86 ) 4 = − 70 + 76 − 80 + 86 4 = 12 4 = 3 \beta_2 = \frac{\sum X_{2i} Y_i}{N} = \frac{(-1)(70) + (+1)(76) + (-1)(80) + (+1)(86)}{4} = \frac{-70 + 76 - 80 + 86}{4} = \frac{12}{4} = 3 β 2 = N ∑ X 2 i Y i = 4 ( − 1 ) ( 70 ) + ( + 1 ) ( 76 ) + ( − 1 ) ( 80 ) + ( + 1 ) ( 86 ) = 4 − 70 + 76 − 80 + 86 = 4 12 = 3
β 2 = 3 \boxed{\beta_2 = 3} β 2 = 3
2. Write the Fitted Mathematical Model Y = 78 + 5 X 1 + 3 X 2 \boxed{Y = 78 + 5X_1 + 3X_2} Y = 78 + 5 X 1 + 3 X 2
where:
78 78 78 = Intercept (baseline expected yield at the center point X 1 = 0 , X 2 = 0 X_1 = 0, X_2 = 0 X 1 = 0 , X 2 = 0 ).
5 5 5 = Main effect coefficient of Temperature (X 1 X_1 X 1 ).
3 3 3 = Main effect coefficient of Pressure (X 2 X_2 X 2 ).
3. Verify the Model Substitute X 1 = 1 , X 2 = 1 X_1 = 1, X_2 = 1 X 1 = 1 , X 2 = 1 (Run 4):
Y ^ = 78 + 5 ( 1 ) + 3 ( 1 ) = 78 + 5 + 3 = 86 \hat{Y} = 78 + 5(1) + 3(1) = 78 + 5 + 3 = 86 Y ^ = 78 + 5 ( 1 ) + 3 ( 1 ) = 78 + 5 + 3 = 86
Observed value: 86 % 86\% 86%
Predicted value: 86 % 86\% 86%
Residual: e = 86 − 86 = 0 e = 86 - 86 = 0 e = 86 − 86 = 0 .
Conclusion: The model fits the orthogonal factorial data points perfectly.
4. Predict Response at Center Point Substitute X 1 = 0 , X 2 = 0 X_1 = 0, X_2 = 0 X 1 = 0 , X 2 = 0 (Temperature = 120 ∘ C = 120^\circ\text{C} = 12 0 ∘ C , Pressure = 15 psi = 15\,\text{psi} = 15 psi ):
Y ^ = 78 + 5 ( 0 ) + 3 ( 0 ) = 78 \hat{Y} = 78 + 5(0) + 3(0) = 78 Y ^ = 78 + 5 ( 0 ) + 3 ( 0 ) = 78
Predicted yield at center point: 78 % \mathbf{78\%} 78% .
5. Direction of Improvement & Factor Significance
Coefficients: Both β 1 = + 5 \beta_1 = +5 β 1 = + 5 and β 2 = + 3 \beta_2 = +3 β 2 = + 3 are positive.
Direction: Increase both Temperature and Pressure to increase yield.
Relative Strength: Since 5 > 3 5 > 3 5 > 3 , Temperature has a stronger positive effect on yield than Pressure.
6. Method of Steepest Ascent Move in the direction of the gradient vector:
∇ Y = ( 5 , 3 ) \nabla Y = (5, 3) ∇ Y = ( 5 , 3 )
Taking step increments proportional to the ( 5 : 3 ) (5 : 3) ( 5 : 3 ) ratio with Δ X 1 = 0.5 \Delta X_1 = 0.5 Δ X 1 = 0.5 and Δ X 2 = 0.3 \Delta X_2 = 0.3 Δ X 2 = 0.3 :
Step Coded X 1 X_1 X 1 Coded X 2 X_2 X 2 Physical Temp (∘ C ^\circ\text{C} ∘ C ) Physical Pressure (psi \text{psi} psi ) Start 0.0 0.0 120.0 ∘ C 120.0^\circ\text{C} 120. 0 ∘ C 15.0 psi 15.0\,\text{psi} 15.0 psi 1 0.5 0.3 130.0 ∘ C 130.0^\circ\text{C} 130. 0 ∘ C 16.5 psi 16.5\,\text{psi} 16.5 psi 2 1.0 0.6 140.0 ∘ C 140.0^\circ\text{C} 140. 0 ∘ C 18.0 psi 18.0\,\text{psi} 18.0 psi 3 1.5 0.9 150.0 ∘ C 150.0^\circ\text{C} 150. 0 ∘ C 19.5 psi 19.5\,\text{psi} 19.5 psi
Continue running experiments along this path until observed yield ceases to increase.
7. Second-Order Model & Confirmation Experiment
Second-Order Model: When curvature exists, the quadratic model accounts for diminishing returns (− 5 X 1 2 , − 4 X 2 2 -5X_1^2, -4X_2^2 − 5 X 1 2 , − 4 X 2 2 ) and factor interaction (− 2 X 1 X 2 -2X_1 X_2 − 2 X 1 X 2 ):
Y = 90 + 4 X 1 + 3 X 2 − 2 X 1 X 2 − 5 X 1 2 − 4 X 2 2 Y = 90 + 4X_1 + 3X_2 - 2X_1 X_2 - 5X_1^2 - 4X_2^2 Y = 90 + 4 X 1 + 3 X 2 − 2 X 1 X 2 − 5 X 1 2 − 4 X 2 2
Confirmation Trial:
Predicted Yield: Y ^ = 92 % \hat{Y} = 92\% Y ^ = 92% at ( X 1 = 0.4 , X 2 = 0.3 ) (X_1 = 0.4, X_2 = 0.3) ( X 1 = 0.4 , X 2 = 0.3 ) .
Actual Experimental Yield: Y actual = 91.5 % Y_{\text{actual}} = 91.5\% Y actual = 91.5% .
Validation Decision: Since the actual experimental yield (91.5 % 91.5\% 91.5% ) closely matches the model prediction (92 % 92\% 92% ) within experimental uncertainty, the fitted model is validated.