
=== Main Effects: yield ===
Factor                   Effect    Std Error   % Contribution
--------------------------------------------------------------
temperature             15.1567       3.0390            46.8%
pressure                12.6667       3.0390            39.1%
catalyst                 4.5550       3.0390            14.1%

=== ANOVA Table: yield ===
Source                      DF           SS           MS          F    p-value
-----------------------------------------------------------------------------
temperature                  1     689.1736     689.1736   2055.445     0.0000
pressure                     1     481.3333     481.3333   1435.566     0.0000
catalyst                     2      46.4647      23.2323     69.290     0.0001
temperature*pressure         1       0.0800       0.0800      0.239     0.6425
Error                        6       2.0118       0.3353
Total                       11    1219.0634     110.8239

=== Interaction Effects: yield ===
Factor A             Factor B              Interaction   % Contribution
------------------------------------------------------------------------
temperature          pressure                  -0.1633           100.0%

=== Summary Statistics: yield ===

temperature:
  Level               N       Mean        Std        Min        Max
  ------------------------------------------------------------
  100                 6    64.8717     7.3484    56.3300    74.2800
  200                 6    80.0283     7.2096    72.2600    89.2600

pressure:
  Level               N       Mean        Std        Min        Max
  ------------------------------------------------------------
  1                   6    66.1167     8.6358    56.3300    76.5400
  5                   6    78.7833     8.5422    69.4000    89.2600

catalyst:
  Level               N       Mean        Std        Min        Max
  ------------------------------------------------------------
  A                   4    71.5400    11.0371    58.3900    85.3300
  B                   4    75.1825    11.7193    60.6500    89.2600
  C                   4    70.6275    11.4763    56.3300    84.2500

=== Model Adequacy (quadratic): yield ===
  R²              = 0.9988
  Adj R²          = 0.9933
  Predicted R²    = 0.9845   (PRESS = 18.8851)
  Shapiro-Wilk    = W=0.916, p=0.252
  Durbin-Watson   = 2.907
  Run-order drift = slope=-0.01686, p=0.609
  Max leverage    = 0.750   (threshold 1.667)
  Max Cook's D    = 0.574   (threshold 1.345)

=== Stationary Point: yield ===
  Nature          : saddle
  Predicted value : 24.1011
  Location        :
    temperature        = 100            (coded -0.151)
    pressure           = 1              (coded -0.127)
    catalyst           = B              (coded -0.060)
  Eigenvalues     : -8.1986, +50.0405, +50.2039

=== Achieved Power: yield ===
  Runs: 12, df_error: 6
  Sigma (sqrt residual MS): 0.5790   (residual MS = 0.3353)
  Alpha: 0.05   |   delta: 1.1581   |   target power: 0.8

  Factor                   Levels  Power @ delta   MDE @ target
  ------------------------ ------ -------------- --------------
  temperature                   2          0.981         0.7959
  pressure                      2          0.981         0.7959
  catalyst                      3          0.786         1.1789
  Note: power < 0.80 for some factors at the chosen delta.

=== Cross-Validation (quadratic, k=5): yield ===
  Note: Skipped: training fold size 9 <= parameters 10. Use --no-rsm or supply more runs.

=== Main Effects: cost ===
Factor                   Effect    Std Error   % Contribution
--------------------------------------------------------------
pressure                 4.4100       0.7181            67.7%
temperature              1.7400       0.7181            26.7%
catalyst                 0.3675       0.7181             5.6%

=== ANOVA Table: cost ===
Source                      DF           SS           MS          F    p-value
-----------------------------------------------------------------------------
temperature                  1       9.0828       9.0828    153.794     0.0000
pressure                     1      58.3443      58.3443    987.910     0.0000
catalyst                     2       0.2793       0.1397      2.365     0.1749
temperature*pressure         1       0.0003       0.0003      0.005     0.9455
Error                        6       0.3544       0.0591
Total                       11      68.0611       6.1874

=== Interaction Effects: cost ===
Factor A             Factor B              Interaction   % Contribution
------------------------------------------------------------------------
temperature          pressure                   0.0100           100.0%

=== Summary Statistics: cost ===

temperature:
  Level               N       Mean        Std        Min        Max
  ------------------------------------------------------------
  100                 6    15.5567     2.4157    13.2400    17.9500
  200                 6    17.2967     2.4413    14.7400    19.7400

pressure:
  Level               N       Mean        Std        Min        Max
  ------------------------------------------------------------
  1                   6    14.2217     0.9930    13.2400    15.5900
  5                   6    18.6317     0.9785    17.5200    19.7400

catalyst:
  Level               N       Mean        Std        Min        Max
  ------------------------------------------------------------
  A                   4    16.6300     2.8277    13.2400    19.7400
  B                   4    16.3875     2.6347    13.5200    19.3000
  C                   4    16.2625     2.7670    13.3100    19.4800

=== Model Adequacy (quadratic): cost ===
  R²              = 0.9958
  Adj R²          = 0.9771
  Predicted R²    = 0.9493   (PRESS = 3.4523)
  Shapiro-Wilk    = W=0.930, p=0.380
  Durbin-Watson   = 2.176
  Run-order drift = slope=0.00271, p=0.851
  Max leverage    = 0.750   (threshold 1.667)
  Max Cook's D    = 0.478   (threshold 1.345)

=== Stationary Point: cost ===
  Nature          : ridge  (outside design region)
  Predicted value : 5.0688
  Location        :
    temperature        = 100            (coded -0.066)
    pressure           = 1              (coded -0.203)
    catalyst           = C              (coded +1.530)
  Eigenvalues     : +0.1167, +10.9205, +10.9303
  Ridge axis      : temperature=-0.009, pressure=+0.001, catalyst=-1.000

=== Achieved Power: cost ===
  Runs: 12, df_error: 6
  Sigma (sqrt residual MS): 0.2430   (residual MS = 0.0591)
  Alpha: 0.05   |   delta: 0.4860   |   target power: 0.8

  Factor                   Levels  Power @ delta   MDE @ target
  ------------------------ ------ -------------- --------------
  temperature                   2          0.981         0.3340
  pressure                      2          0.981         0.3340
  catalyst                      3          0.786         0.4948
  Note: power < 0.80 for some factors at the chosen delta.

=== Cross-Validation (quadratic, k=5): cost ===
  Note: Skipped: training fold size 9 <= parameters 10. Use --no-rsm or supply more runs.

Pareto chart (yield): results/pareto_yield.png

Pareto chart (cost): results/pareto_cost.png
Main effects (yield): results/main_effects_yield.png
Main effects (cost): results/main_effects_cost.png
HTML report written to: results/report.html

HTML report: results/report.html
