Skip to main content

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,CA, B, C) are applied to 12 homogeneous experimental plots (t=3,r=4,N=12t = 3, r = 4, N = 12). The measured crop yields (in kg) are:

  • Fertilizer A: 18,20,19,2218, 20, 19, 22
  • Fertilizer B: 24,22,21,2324, 22, 21, 23
  • Fertilizer C: 26,25,27,2826, 25, 27, 28

Test at the α=0.05\alpha = 0.05 significance level whether there is a statistically significant difference in mean yield among the three fertilizers. (Given: Critical F0.05(2,9)≈4.26F_{0.05}(2, 9) \approx 4.26).

View Step-by-Step Solution

1. Hypothesis Formulation​

  • Null Hypothesis (H0H_0): μA=μB=μC\mu_A = \mu_B = \mu_C (The mean yields of all three fertilizers are equal).
  • Alternative Hypothesis (H1H_1): At least one fertilizer mean yield is different.

2. Totals & Correction Factor (CFCF)​

  • TA=18+20+19+22=79T_A = 18 + 20 + 19 + 22 = 79
  • TB=24+22+21+23=90T_B = 24 + 22 + 21 + 23 = 90
  • TC=26+25+27+28=106T_C = 26 + 25 + 27 + 28 = 106
  • Grand Total G=79+90+106=275\text{Grand Total } G = 79 + 90 + 106 = 275
  • CF=G2N=275212=7562512≈6302.0833CF = \frac{G^2}{N} = \frac{275^2}{12} = \frac{75625}{12} \approx 6302.0833

3. Sum of Squares Computations​

  • ∑yij2=182+202+192+222+242+222+212+232+262+252+272+282=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
  • SST=∑yij2−CF=6413−6302.0833=110.9167\text{SST} = \sum y_{ij}^2 - CF = 6413 - 6302.0833 = 110.9167
  • SSTreat=792+902+10624−CF=6241+8100+112364−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
  • SSE=SST−SSTreat=110.9167−92.1667=18.7500\text{SSE} = \text{SST} - \text{SSTreat} = 110.9167 - 92.1667 = 18.7500

4. CRD ANOVA Table​

Sourcedf\text{df}SS\text{SS}MS\text{MS}Calculated FFTabulated F0.05F_{0.05}
Fertilizers (Treatments)3−1=23 - 1 = 292.166792.166792.16672=46.0833\frac{92.1667}{2} = 46.083346.08332.0833≈22.12\frac{46.0833}{2.0833} \approx 22.124.264.26
Error (Residual)12−3=912 - 3 = 918.750018.750018.75009=2.0833\frac{18.7500}{9} = 2.0833--
Total12−1=1112 - 1 = 11110.9167110.9167---

5. Decision & Academic Conclusion​

  • Comparison: Fcal=22.12>Ftab=4.26F_{\text{cal}} = 22.12 > F_{\text{tab}} = 4.26.
  • Decision: Reject H0H_0 at the 5%5\% level of significance.
  • Conclusion: "Since the calculated FF-value (22.1222.12) is greater than the tabulated FF-value (4.264.26) at the 5%5\% level of significance, the null hypothesis (H0H_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 (D1,D2,D3,D4D_1, D_2, D_3, D_4) tested across 4 distinct laboratories (A,B,C,DA, B, C, D) to control for inter-laboratory testing bias (t=4,r=4,N=16t = 4, r = 4, N = 16).

Laboratory (Block)D1D_1D2D_2D3D_3D4D_4
A2543
B3742
C5564
D5585

Analyze the data using RBD and state your conclusion at α=0.05\alpha = 0.05. (Given: Critical F0.05(3,9)=3.86F_{0.05}(3, 9) = 3.86).

View Step-by-Step Solution

1. Totals Calculation​

  • Treatment Totals (TiT_i):
    • TD1=2+3+5+5=15T_{D1} = 2 + 3 + 5 + 5 = 15
    • TD2=5+7+5+5=22T_{D2} = 5 + 7 + 5 + 5 = 22
    • TD3=4+4+6+8=22T_{D3} = 4 + 4 + 6 + 8 = 22
    • TD4=3+2+4+5=14T_{D4} = 3 + 2 + 4 + 5 = 14
  • Block Totals (BjB_j):
    • BA=2+5+4+3=14B_A = 2 + 5 + 4 + 3 = 14
    • BB=3+7+4+2=16B_B = 3 + 7 + 4 + 2 = 16
    • BC=5+5+6+4=20B_C = 5 + 5 + 6 + 4 = 20
    • BD=5+5+8+5=23B_D = 5 + 5 + 8 + 5 = 23
  • Grand Total G=15+22+22+14=73\text{Grand Total } G = 15 + 22 + 22 + 14 = 73
  • CF=G2N=73216=532916=333.0625CF = \frac{G^2}{N} = \frac{73^2}{16} = \frac{5329}{16} = 333.0625

2. Sum of Squares Partitioning​

  • ∑yij2=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
  • SST=373−333.0625=39.9375\text{SST} = 373 - 333.0625 = 39.9375
  • SSTreat=152+222+222+1424−CF=13894−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
  • SSBlock=142+162+202+2324−CF=13814−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
  • 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

3. RBD ANOVA Table​

Sourcedf\text{df}SS\text{SS}MS\text{MS}Calculated FFCritical F0.05F_{0.05}
Diets (Treatments)4−1=34 - 1 = 314.187514.187514.18753≈4.7292\frac{14.1875}{3} \approx 4.72924.72921.5069≈3.14\frac{4.7292}{1.5069} \approx 3.143.863.86
Labs (Blocks)4−1=34 - 1 = 312.187512.187512.18753≈4.0625\frac{12.1875}{3} \approx 4.0625--
Error(4−1)(4−1)=9(4-1)(4-1)=913.562513.562513.56259≈1.5069\frac{13.5625}{9} \approx 1.5069--
Total16−1=1516 - 1 = 1539.937539.9375---

4. Conclusion​

  • Fcal=3.14<Ftab=3.86F_{\text{cal}} = 3.14 < F_{\text{tab}} = 3.86.
  • Fail to reject H0H_0. At the 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,DA, B, C, D) are allocated in a 4×44 \times 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)
  • Row 2: 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)
  • Row 4: D(20),A(21),B(23),C(27)D(20), A(21), B(23), C(27)

Test at α=0.05\alpha = 0.05 whether treatment differences are significant. (Given: Critical F0.05(3,6)=4.76F_{0.05}(3, 6) = 4.76, ∑y2=8408\sum y^2 = 8408).

View Step-by-Step Solution

1. Totals Calculation (p=4,N=16p = 4, N = 16)​

  • Treatment Totals:
    • TA=18+20+22+21=81T_A = 18 + 20 + 22 + 21 = 81
    • TB=22+21+25+23=91T_B = 22 + 21 + 25 + 23 = 91
    • TC=25+24+19+27=95T_C = 25 + 24 + 19 + 27 = 95
    • TD=28+26+23+20=97T_D = 28 + 26 + 23 + 20 = 97
  • Row Totals:
    • R1=18+22+25+28=93R_1 = 18 + 22 + 25 + 28 = 93
    • R2=21+24+26+20=91R_2 = 21 + 24 + 26 + 20 = 91
    • R3=19+23+22+25=89R_3 = 19 + 23 + 22 + 25 = 89
    • R4=20+21+23+27=91R_4 = 20 + 21 + 23 + 27 = 91
  • Column Totals:
    • C1=18+21+19+20=78C_1 = 18 + 21 + 19 + 20 = 78
    • C2=22+24+23+21=90C_2 = 22 + 24 + 23 + 21 = 90
    • C3=25+26+22+23=96C_3 = 25 + 26 + 22 + 23 = 96
    • C4=28+20+25+27=100C_4 = 28 + 20 + 25 + 27 = 100
  • G=364G = 364
  • CF=364216=8281.0CF = \frac{364^2}{16} = 8281.0

2. Sum of Squares Partitioning​

  • SST=8408.0−8281.0=127.0\text{SST} = 8408.0 - 8281.0 = 127.0
  • SSTreat=812+912+952+9724−8281.0=332764−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
  • SSRow=932+912+892+9124−8281.0=331324−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
  • SSCol=782+902+962+10024−8281.0=334004−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
  • SSE=127.0−38.0−2.0−69.0=18.0\text{SSE} = 127.0 - 38.0 - 2.0 - 69.0 = 18.0

3. LSD ANOVA Table​

Sourcedf\text{df}SS\text{SS}MS\text{MS}Calculated FFCritical F0.05F_{0.05}
Rows4−1=34 - 1 = 32.02.00.66670.6667--
Columns4−1=34 - 1 = 369.069.023.000023.0000--
Treatments4−1=34 - 1 = 338.038.012.666712.666712.66673.0000≈4.222\frac{12.6667}{3.0000} \approx 4.2224.764.76
Error(4−1)(4−2)=6(4-1)(4-2)=618.018.03.00003.0000--
Total16−1=1516 - 1 = 15127.0127.0---

4. Conclusion​

  • Fcal=4.222<Ftab=4.76F_{\text{cal}} = 4.222 < F_{\text{tab}} = 4.76.
  • Fail to reject H0H_0. After adjusting for row and column effects, there is no statistically significant treatment effect at the 5%5\% level.

Problem 4: 222^2 Full Factorial Design​

Problem Statement​

A chemical engineer investigates the effect of Temperature (AA) and Pressure (BB) on product yield (%). Both factors are set at two coded levels: Low (−1-1) and High (+1+1) with r=2r=2 replicates (N=22×2=8N = 2^2 \times 2 = 8).

RunFactor A (Temp)Factor B (Pressure)Replicate 1Replicate 2
1−1-1−1-14243
2+1+1−1-15049
3−1-1+1+14645
4+1+1+1+16261

Perform a complete 222^2 factorial analysis: compute contrasts, main effects, interaction effect, construct the ANOVA table, and test for significance at α=0.05\alpha = 0.05. (Given: Critical F0.05(1,4)=7.708F_{0.05}(1, 4) = 7.708).

View Step-by-Step Solution

1. Treatment Totals​

  • Run 1 (−1,−1-1, -1): y1=42+43=85y_1 = 42 + 43 = 85
  • Run 2 (+1,−1+1, -1): y2=50+49=99y_2 = 50 + 49 = 99
  • Run 3 (−1,+1-1, +1): y3=46+45=91y_3 = 46 + 45 = 91
  • Run 4 (+1,+1+1, +1): y4=62+61=123y_4 = 62 + 61 = 123
  • Grand Total G=85+99+91+123=398\text{Grand Total } G = 85 + 99 + 91 + 123 = 398

2. Contrasts & Effects (2k−1⋅r=22−1⋅2=42^{k-1} \cdot r = 2^{2-1} \cdot 2 = 4)​

  • CA=−85+99−91+123=46  ⟹  A^=464=11.5C_A = -85 + 99 - 91 + 123 = 46 \implies \hat{A} = \frac{46}{4} = 11.5
  • CB=−85−99+91+123=30  ⟹  B^=304=7.5C_B = -85 - 99 + 91 + 123 = 30 \implies \hat{B} = \frac{30}{4} = 7.5
  • CAB=+85−99−91+123=18  ⟹  AB^=184=4.5C_{AB} = +85 - 99 - 91 + 123 = 18 \implies \widehat{AB} = \frac{18}{4} = 4.5

3. Sum of Squares (2k⋅r=82^k \cdot r = 8)​

  • SSA=4628=264.5\text{SS}_A = \frac{46^2}{8} = 264.5
  • SSB=3028=112.5\text{SS}_B = \frac{30^2}{8} = 112.5
  • SSAB=1828=40.5\text{SS}_{AB} = \frac{18^2}{8} = 40.5
  • ∑y2=422+432+⋯+612=20220\sum y^2 = 42^2 + 43^2 + \dots + 61^2 = 20220
  • CF=39828=19800.5  ⟹  SST=20220−19800.5=419.5CF = \frac{398^2}{8} = 19800.5 \implies \text{SST} = 20220 - 19800.5 = 419.5
  • SSE=SST−(SSA+SSB+SSAB)=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

4. Factorial ANOVA Table​

Sourcedf\text{df}SS\text{SS}MS\text{MS}Calculated FFCritical F0.05F_{0.05}Decision
A (Temperature)1264.5264.5264.5264.5264.50.5=529.0\frac{264.5}{0.5} = 529.07.7087.708Reject H0H_0
B (Pressure)1112.5112.5112.5112.5112.50.5=225.0\frac{112.5}{0.5} = 225.07.7087.708Reject H0H_0
AB (Interaction)140.540.540.540.540.50.5=81.0\frac{40.5}{0.5} = 81.07.7087.708Reject H0H_0
Error (Residual)42.02.0MSE=0.5\text{MSE} = 0.5---
Total7419.5419.5----

5. Conclusion​

All factors (A,B,A, B, and interaction ABAB) are statistically significant at p<0.001p < 0.001. Temperature increases yield by 11.5%11.5\%, Pressure increases yield by 7.5%7.5\%, and positive interaction indicates the benefits of temperature are greatest when paired with high pressure.


Problem 5: Taguchi L9L_9 Orthogonal Array Optimization​

Problem Statement​

A manufacturing engineer seeks to optimize three 3-level parameters to maximize tensile strength (yy in MPa\text{MPa}):

  • Factor A (Temperature): Level 1 = 200∘C200^\circ\text{C}, Level 2 = 220∘C220^\circ\text{C}, Level 3 = 240∘C240^\circ\text{C}
  • Factor B (Holding Time): Level 1 = 10 min10\,\text{min}, Level 2 = 20 min20\,\text{min}, Level 3 = 30 min30\,\text{min}
  • Factor C (Pressure): Level 1 = 5 bar5\,\text{bar}, Level 2 = 7 bar7\,\text{bar}, Level 3 = 9 bar9\,\text{bar}

The measured single-run responses in an L9L_9 array are:

  • Run 1: (1, 1, 1)   ⟹  y=42 MPa\implies y = 42\,\text{MPa}
  • Run 2: (1, 2, 2)   ⟹  y=48 MPa\implies y = 48\,\text{MPa}
  • Run 3: (1, 3, 3)   ⟹  y=50 MPa\implies y = 50\,\text{MPa}
  • Run 4: (2, 1, 2)   ⟹  y=55 MPa\implies y = 55\,\text{MPa}
  • Run 5: (2, 2, 3)   ⟹  y=60 MPa\implies y = 60\,\text{MPa}
  • Run 6: (2, 3, 1)   ⟹  y=58 MPa\implies y = 58\,\text{MPa}
  • Run 7: (3, 1, 3)   ⟹  y=65 MPa\implies y = 65\,\text{MPa}
  • Run 8: (3, 2, 1)   ⟹  y=62 MPa\implies y = 62\,\text{MPa}
  • Run 9: (3, 3, 2)   ⟹  y=59 MPa\implies y = 59\,\text{MPa}

Using the Larger-the-better S/NS/N ratio formulation:

  1. Compute the S/NS/N ratio for each trial.
  2. Build the factor-level response table and rank the factors.
  3. Determine the optimal parameter combination.
  4. Predict the optimum tensile strength performance.
View Step-by-Step Solution

1. S/N Computation (S/N=20log⁡10(y)S/N = 20 \log_{10}(y))​

  • Run 1:20log⁡10(42)=32.4650 dB\text{Run 1}: 20 \log_{10}(42) = 32.4650\,\text{dB}
  • Run 2:20log⁡10(48)=33.6248 dB\text{Run 2}: 20 \log_{10}(48) = 33.6248\,\text{dB}
  • Run 3:20log⁡10(50)=33.9794 dB\text{Run 3}: 20 \log_{10}(50) = 33.9794\,\text{dB}
  • Run 4:20log⁡10(55)=34.8073 dB\text{Run 4}: 20 \log_{10}(55) = 34.8073\,\text{dB}
  • Run 5:20log⁡10(60)=35.5630 dB\text{Run 5}: 20 \log_{10}(60) = 35.5630\,\text{dB}
  • Run 6:20log⁡10(58)=35.2686 dB\text{Run 6}: 20 \log_{10}(58) = 35.2686\,\text{dB}
  • Run 7:20log⁡10(65)=36.2583 dB\text{Run 7}: 20 \log_{10}(65) = 36.2583\,\text{dB}
  • Run 8:20log⁡10(62)=35.8478 dB\text{Run 8}: 20 \log_{10}(62) = 35.8478\,\text{dB}
  • Run 9:20log⁡10(59)=35.4170 dB\text{Run 9}: 20 \log_{10}(59) = 35.4170\,\text{dB}
  • Overall Average: S/N‾=34.8035 dB\overline{S/N} = 34.8035\,\text{dB}

2. Factor-Level S/N Response Table​

  • Factor A (Temp):
    • Level 1: 32.4650+33.6248+33.97943=33.3564 dB\frac{32.4650 + 33.6248 + 33.9794}{3} = 33.3564\,\text{dB}
    • Level 2: 34.8073+35.5630+35.26863=35.2130 dB\frac{34.8073 + 35.5630 + 35.2686}{3} = 35.2130\,\text{dB}
    • Level 3: 36.2583+35.8478+35.41703=35.8410 dB\frac{36.2583 + 35.8478 + 35.4170}{3} = 35.8410\,\text{dB}
    • ΔA=35.8410−33.3564=2.4846 dB\Delta_A = 35.8410 - 33.3564 = 2.4846\,\text{dB} (Rank 1)
  • Factor B (Time):
    • Level 1: 32.4650+34.8073+36.25833=34.5102 dB\frac{32.4650 + 34.8073 + 36.2583}{3} = 34.5102\,\text{dB}
    • Level 2: 33.6248+35.5630+35.84783=35.0119 dB\frac{33.6248 + 35.5630 + 35.8478}{3} = 35.0119\,\text{dB}
    • Level 3: 33.9794+35.2686+35.41703=34.8883 dB\frac{33.9794 + 35.2686 + 35.4170}{3} = 34.8883\,\text{dB}
    • ΔB=35.0119−34.5102=0.5017 dB\Delta_B = 35.0119 - 34.5102 = 0.5017\,\text{dB} (Rank 3)
  • Factor C (Pressure):
    • Level 1: 32.4650+35.2686+35.84783=34.5271 dB\frac{32.4650 + 35.2686 + 35.8478}{3} = 34.5271\,\text{dB}
    • Level 2: 33.6248+34.8073+35.41703=34.6164 dB\frac{33.6248 + 34.8073 + 35.4170}{3} = 34.6164\,\text{dB}
    • Level 3: 33.9794+35.5630+36.25833=35.2669 dB\frac{33.9794 + 35.5630 + 36.2583}{3} = 35.2669\,\text{dB}
    • ΔC=35.2669−34.5271=0.7398 dB\Delta_C = 35.2669 - 34.5271 = 0.7398\,\text{dB} (Rank 2)

3. Optimal Setting Selection​

Always pick the level with the highest mean S/NS/N:

  • Factor A: Level 3 (240∘C240^\circ\text{C})
  • Factor B: Level 2 (20 min20\,\text{min})
  • Factor C: Level 3 (9 bar9\,\text{bar})
  • Optimal Setting: A3B2C3A_3 B_2 C_3

4. Predicted Optimum Performance​

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} Predicted Tensile Strength: y^=1036.512820≈66.93 MPa\text{Predicted Tensile Strength: } \hat{y} = 10^{\frac{36.5128}{20}} \approx 66.93\,\text{MPa}


Problem 6: Genetic Algorithm Optimization​

Problem Statement​

A Genetic Algorithm is used to maximize the function f(x)=x2f(x) = x^2 over the integer domain 0≤x≤310 \le x \le 31 encoded as 5-bit binary strings. The initial random population consists of:

  • C1=101102C_1 = 10110_2
  • C2=001012C_2 = 00101_2
  • C3=111002C_3 = 11100_2
  • C4=010112C_4 = 01011_2

Perform one full generation of the algorithm:

  1. Decode each chromosome to decimal and evaluate its fitness.
  2. Select the top two chromosomes as parents.
  3. Perform single-point crossover after the 2nd bit to produce two children.
  4. Mutate the last bit of Child 2.
  5. 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)=x2f(x) = x^2)​

  • C1=101102=16+4+2=22  ⟹  f(22)=222=484C_1 = 10110_2 = 16 + 4 + 2 = 22 \implies f(22) = 22^2 = 484
  • C2=001012=4+1=5  ⟹  f(5)=52=25C_2 = 00101_2 = 4 + 1 = 5 \implies f(5) = 5^2 = 25
  • C3=111002=16+8+4=28  ⟹  f(28)=282=784C_3 = 11100_2 = 16 + 8 + 4 = 28 \implies f(28) = 28^2 = 784
  • C4=010112=8+2+1=11  ⟹  f(11)=112=121C_4 = 01011_2 = 8 + 2 + 1 = 11 \implies f(11) = 11^2 = 121
  • Gen 0 Average Fitness: fˉ0=484+25+784+1214=353.5\bar{f}_0 = \frac{484 + 25 + 784 + 121}{4} = 353.5

2. Parent Selection​

The two highest fitness individuals are:

  • Parent 1: C3=11100C_3 = 11100 (x=28,f=784x = 28, f = 784)
  • Parent 2: C1=10110C_1 = 10110 (x=22,f=484x = 22, f = 484)

3. Crossover (Cut after Bit 2)​

  • Parent 1: [11∣100]\text{Parent 1: } [11 \mid 100]
  • Parent 2: [10∣110]\text{Parent 2: } [10 \mid 110]
  • Child 1: [11∣110]=111102=30  ⟹  f(30)=302=900\text{Child 1: } [11 \mid 110] = 11110_2 = 30 \implies f(30) = 30^2 = 900
  • Child 2: [10∣100]=101002=20  ⟹  f(20)=202=400\text{Child 2: } [10 \mid 100] = 10100_2 = 20 \implies f(20) = 20^2 = 400

4. Mutation​

Flip the 5th bit of Child 2 (10100→1010110100 \to 10101):

  • Mutated Child 2=101012=21  ⟹  f(21)=212=441\text{Mutated Child 2} = 10101_2 = 21 \implies f(21) = 21^2 = 441

5. New Generation Formation (Elitism)​

  1. Child 1: 111102  ⟹  x=30,f=90011110_2 \implies x = 30, f = 900
  2. Mutated Child 2: 101012  ⟹  x=21,f=44110101_2 \implies x = 21, f = 441
  3. Elite Parent 1: 111002  ⟹  x=28,f=78411100_2 \implies x = 28, f = 784
  4. Elite Parent 2: 101102  ⟹  x=22,f=48410110_2 \implies x = 22, f = 484
  • Gen 1 Average Fitness: fˉ1=900+441+784+4844=26094=652.25\bar{f}_1 = \frac{900 + 441 + 784 + 484}{4} = \frac{2609}{4} = 652.25
  • Improvement: Average fitness increased by 652.25−353.5353.5×100%≈84.5%\frac{652.25 - 353.5}{353.5} \times 100\% \approx 84.5\%.

Problem 7: Response Surface Methodology (RSM) Optimization​

Problem Statement​

A chemical company wants to maximize product yield (YY, in %\%). Two continuous factors are investigated:

  • Temperature (X1X_1): Low level (−1-1) = 100∘C100^\circ\text{C}, High level (+1+1) = 140∘C140^\circ\text{C} (Center point 0=120∘C0 = 120^\circ\text{C}).
  • Pressure (X2X_2): Low level (−1-1) = 10 psi10\,\text{psi}, High level (+1+1) = 20 psi20\,\text{psi} (Center point 0=15 psi0 = 15\,\text{psi}).

The experimental results from an initial 222^2 factorial design matrix are:

RunX1X_1 (Temperature)X2X_2 (Pressure)Yield YY (%)
1−1-1−1-170
2−1-1+1+176
3+1+1−1-180
4+1+1+1+186

Perform the following:

  1. Assume a first-order model Y=β0+β1X1+β2X2Y = \beta_0 + \beta_1 X_1 + \beta_2 X_2 and calculate the coefficients β0,β1,β2\beta_0, \beta_1, \beta_2.
  2. Write the fitted mathematical model.
  3. Verify the model at run (X1=+1,X2=+1)(X_1 = +1, X_2 = +1).
  4. Predict the yield at the center point (X1=0,X2=0)(X_1 = 0, X_2 = 0).
  5. Determine the direction of improvement and identify which factor has the stronger effect.
  6. 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 ΔX1=0.5\Delta X_1 = 0.5 and ΔX2=0.3\Delta X_2 = 0.3.
  7. Suppose subsequent trials exhibit curvature and a second-order model is fitted: Y=90+4X1+3X2−2X1X2−5X12−4X22Y = 90 + 4X_1 + 3X_2 - 2X_1 X_2 - 5X_1^2 - 4X_2^2 If mathematical optimization predicts an optimum at (X1=0.4,X2=0.3)(X_1 = 0.4, X_2 = 0.3) yielding 92%92\%, validate the model if an actual confirmation trial produces 91.5%91.5\%.
View Step-by-Step Solution

1. Calculate First-Order Regression Coefficients​

Using the orthogonal properties of the 222^2 coded matrix (N=4N = 4):

  • Calculate β0\beta_0 (Intercept / Average Response): β0=∑YiN=70+76+80+864=3124=78\beta_0 = \frac{\sum Y_i}{N} = \frac{70 + 76 + 80 + 86}{4} = \frac{312}{4} = 78 β0=78\boxed{\beta_0 = 78}

  • Calculate β1\beta_1 (Main Effect of Temperature X1X_1): β1=∑X1iYiN=(−1)(70)+(−1)(76)+(+1)(80)+(+1)(86)4=−70−76+80+864=204=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=5\boxed{\beta_1 = 5}

  • Calculate β2\beta_2 (Main Effect of Pressure X2X_2): β2=∑X2iYiN=(−1)(70)+(+1)(76)+(−1)(80)+(+1)(86)4=−70+76−80+864=124=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=3\boxed{\beta_2 = 3}

2. Write the Fitted Mathematical Model​

Y=78+5X1+3X2\boxed{Y = 78 + 5X_1 + 3X_2} where:

  • 7878 = Intercept (baseline expected yield at the center point X1=0,X2=0X_1 = 0, X_2 = 0).
  • 55 = Main effect coefficient of Temperature (X1X_1).
  • 33 = Main effect coefficient of Pressure (X2X_2).

3. Verify the Model​

Substitute X1=1,X2=1X_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

  • Observed value: 86%86\%
  • Predicted value: 86%86\%
  • Residual: e=86−86=0e = 86 - 86 = 0.
  • Conclusion: The model fits the orthogonal factorial data points perfectly.

4. Predict Response at Center Point​

Substitute X1=0,X2=0X_1 = 0, X_2 = 0 (Temperature =120∘C= 120^\circ\text{C}, Pressure =15 psi= 15\,\text{psi}): Y^=78+5(0)+3(0)=78\hat{Y} = 78 + 5(0) + 3(0) = 78

  • Predicted yield at center point: 78%\mathbf{78\%}.

5. Direction of Improvement & Factor Significance​

  • Coefficients: Both β1=+5\beta_1 = +5 and β2=+3\beta_2 = +3 are positive.
  • Direction: Increase both Temperature and Pressure to increase yield.
  • Relative Strength: Since 5>35 > 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) Taking step increments proportional to the (5:3)(5 : 3) ratio with ΔX1=0.5\Delta X_1 = 0.5 and ΔX2=0.3\Delta X_2 = 0.3:

StepCoded X1X_1Coded X2X_2Physical Temp (∘C^\circ\text{C})Physical Pressure (psi\text{psi})
Start0.00.0120.0∘C120.0^\circ\text{C}15.0 psi15.0\,\text{psi}
10.50.3130.0∘C130.0^\circ\text{C}16.5 psi16.5\,\text{psi}
21.00.6140.0∘C140.0^\circ\text{C}18.0 psi18.0\,\text{psi}
31.50.9150.0∘C150.0^\circ\text{C}19.5 psi19.5\,\text{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 (−5X12,−4X22-5X_1^2, -4X_2^2) and factor interaction (−2X1X2-2X_1 X_2): Y=90+4X1+3X2−2X1X2−5X12−4X22Y = 90 + 4X_1 + 3X_2 - 2X_1 X_2 - 5X_1^2 - 4X_2^2
  • Confirmation Trial:
    • Predicted Yield: Y^=92%\hat{Y} = 92\% at (X1=0.4,X2=0.3)(X_1 = 0.4, X_2 = 0.3).
    • Actual Experimental Yield: Yactual=91.5%Y_{\text{actual}} = 91.5\%.
  • Validation Decision: Since the actual experimental yield (91.5%91.5\%) closely matches the model prediction (92%92\%) within experimental uncertainty, the fitted model is validated.