Next Page Previous Page Six Sigma Home Tools & Aids Search Handbook
4. Process Modeling
4.6. Case Studies in Process Modeling
4.6.4. Thermal Expansion of Copper Case Study

4.6.4.5.

Cubic/Cubic Rational Function Model

C/C Rational Function Model Since the Q/Q model was inadequate, we next attempt to fit a cubic/cubic (C/C) rational function model.

We used Dataplot to fit the C/C rational function model with the following 7 subset points to generate the starting values.

 TEMP        THERMEXP
 ----        --------
  10             0
  30             2
  40             3
  50             5
 120            12
 200            15
 800            20
       
Exact Rational Fit Output Dataplot generated the following output from the exact rational fit command. The output has been edited for display.
EXACT RATIONAL FUNCTION FIT
NUMBER OF POINTS IN FIRST SET    =        7
DEGREE OF NUMERATOR              =        3
DEGREE OF DENOMINATOR            =        3
  
NUMERATOR  --A0  A1  A2  A3      =
   -0.2322993E+01  0.3528976E+00 -0.1382551E-01  0.1765684E-03
DENOMINATOR--B0  B1  B2  B3      = 
    0.1000000E+01 -0.3394208E-01  0.1099545E-03  0.7905308E-05
  
APPLICATION OF EXACT-FIT COEFFICIENTS
TO SECOND PAIR OF VARIABLES--
  
NUMBER OF POINTS IN SECOND SET           =      236
NUMBER OF ESTIMATED COEFFICIENTS         =        7
RESIDUAL DEGREES OF FREEDOM              =      229
  
RESIDUAL SUM OF SQUARES                  =  0.78246452E+02
RESIDUAL STANDARD DEVIATION (DENOM=N-P)  =  0.58454049E+00
AVERAGE ABSOLUTE RESIDUAL   (DENOM=N)    =  0.46998626E+00
LARGEST (IN MAGNITUDE) POSITIVE RESIDUAL =  0.95733070E+00
LARGEST (IN MAGNITUDE) NEGATIVE RESIDUAL = -0.13497944E+01
LARGEST (IN MAGNITUDE) ABSOLUTE RESIDUAL =  0.13497944E+01
  
  
       
The important information in this output is the estimates for A0, A1, A2, A3, B1, B2, and B3 (B0 is always set to 1). These values are used as the starting values for the fit in the next section.
Nonlinear Fit Output Dataplot generated the following output for the nonlinear fit. The output has been edited for display.

LEAST SQUARES NON-LINEAR FIT
SAMPLE SIZE N =      236
MODEL--THERMEXP =(A0+A1*TEMP+A2*TEMP**2+A3*TEMP**3)/
                 (1+B1*TEMP+B2*TEMP**2+B3*TEMP**3)
REPLICATION CASE
REPLICATION STANDARD DEVIATION =     0.8131711930D-01
REPLICATION DEGREES OF FREEDOM =           1
NUMBER OF DISTINCT SUBSETS     =         235
  
FINAL PARAMETER ESTIMATES           (APPROX. ST. DEV.)    T VALUE
1  A0                   1.07913       (0.1710    )          6.3
2  A1                 -0.122801       (0.1203E-01)         -10.
3  A2                  0.408837E-02   (0.2252E-03)          18.
4  A3                 -0.142848E-05   (0.2610E-06)         -5.5
5  B1                 -0.576111E-02   (0.2468E-03)         -23.
6  B2                  0.240629E-03   (0.1060E-04)          23.
7  B3                 -0.123254E-06   (0.1217E-07)         -10.

RESIDUAL    STANDARD DEVIATION =         0.0818038210
RESIDUAL    DEGREES OF FREEDOM =         229
REPLICATION STANDARD DEVIATION =         0.0813171193
REPLICATION DEGREES OF FREEDOM =           1
LACK OF FIT F RATIO =       1.0121 = THE  32.1265% POINT OF THE
F DISTRIBUTION WITH    228 AND      1 DEGREES OF FREEDOM

       
The above output yields the following estimated model.
Plot of C/C Rational Function Fit We generate a plot of the fitted rational function model with the raw data.

Looking at the fitted function with the raw data appears to show a reasonable fit.

6-Plot for Model Validation Although the plot of the fitted function with the raw data appears to show a reasonable fit, we need to validate the model assumptions. The 6-plot is an effective tool for this purpose.

The 6-plot indicates no significant violation of the model assumptions. That is, the residuals have constant location and scale (from the residual plot in row 1, column 2), seem to be random (from the lag plot in row 2, column 1), and approximated well by a normal distribution (from the histogram and normal probability plots in row 2, columns 2 and 3).

Residual Plot We generate a full size residual plot in order to show more detail.

The full size residual plot shows that the assumptions of constant location and scale for the residuals are valid. No distinguishing pattern is evident in the residuals.

Conclusion We conclude that the cubic/cubic rational function model does in fact provide a satisfactory model for this data set.
Six Sigma Home Tools & Aids Search Handbook Previous Page Next Page