(***********************************************************************

                    Mathematica-Compatible Notebook

This notebook can be used on any computer system with Mathematica 4.0,
MathReader 4.0, or any compatible application. The data for the notebook 
starts with the line containing stars above.

To get the notebook into a Mathematica-compatible application, do one of 
the following:

* Save the data starting with the line of stars above into a file
  with a name ending in .nb, then open the file inside the application;

* Copy the data starting with the line of stars above to the
  clipboard, then use the Paste menu command inside the application.

Data for notebooks contains only printable 7-bit ASCII and can be
sent directly in email or through ftp in text mode.  Newlines can be
CR, LF or CRLF (Unix, Macintosh or MS-DOS style).

NOTE: If you modify the data for this notebook not in a Mathematica-
compatible application, you must delete the line below containing the 
word CacheID, otherwise Mathematica-compatible applications may try to 
use invalid cache data.

For more information on notebooks and Mathematica-compatible 
applications, contact Wolfram Research:
  web: http://www.wolfram.com
  email: info@wolfram.com
  phone: +1-217-398-0700 (U.S.)

Notebook reader applications are available free of charge from 
Wolfram Research.
***********************************************************************)

(*CacheID: 232*)


(*NotebookFileLineBreakTest
NotebookFileLineBreakTest*)
(*NotebookOptionsPosition[     52606,       1400]*)
(*NotebookOutlinePosition[     59990,       1626]*)
(*  CellTagsIndexPosition[     59251,       1593]*)
(*WindowFrame->Normal*)



Notebook[{
Cell[TextData[{
  "Statistics",
  StyleBox["`", "MB"],
  "NonlinearFit",
  StyleBox["`", "MB"]
}], "Subsection",
  CellTags->{"S6.5.1", "1.1"}],

Cell[TextData[{
  "The built\[Hyphen]in function ",
  StyleBox["Fit", "MR"],
  " finds a least\[Hyphen]squares fit to a list of data for a model that is a \
linear combination of the given basis functions; the coefficients of the \
linear combination are the parameters of the fit. Often a more sophisticated \
model is desired, where the model is not linear in the parameters. The \
function ",
  StyleBox["NonlinearFit", "MR"],
  " allows you to perform this kind of least\[Hyphen]squares fit. ",
  StyleBox["NonlinearFit", "MR"],
  " gives an expression for the model, similar to the output of the built\
\[Hyphen]in function ",
  StyleBox["Fit", "MR"],
  ". The function ",
  StyleBox["NonlinearRegress", "MR"],
  " gives a number of regression diagnostics and allows you to specify \
exactly what will be included in the output, as does the package function ",
  StyleBox["Regress", "MR"],
  ". "
}], "Text",
  CellTags->{"S6.5.1", "1.2"}],

Cell[BoxData[
    FormBox[
      StyleBox[
        FrameBox[
          StyleBox[GridBox[{
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"NonlinearFit[\\\\\\\"\\\",\
 \\\"MR\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\"data\\\\\\\"\\\", \\\"TI\\\"]\\)\
\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \
\\!\\(\\*StyleBox[\\\"\\\\\\\"model\\\\\\\"\\\", \
\\\"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \\!\
\\(\\*StyleBox[\\\"\\\\\\\"variables\\\\\\\"\\\", \
\\\"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \\!\
\\(\\*StyleBox[\\\"\\\\\\\"parameters\\\\\\\"\\\", \
\\\"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\"]\\\\\\\"\\\", \\\"MR\\\"]\\) \
\"\>", \("fit the data to the model with the named variables and parameters, \
returning the model evaluated at the parameter estimates achieving the least\
\[Hyphen]squares fit"\)},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"NonlinearRegress[\\\\\\\"\
\\\", \\\"MR\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\"data\\\\\\\"\\\", \
\\\"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \\!\
\\(\\*StyleBox[\\\"\\\\\\\"model\\\\\\\"\\\", \
\\\"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \\!\
\\(\\*StyleBox[\\\"\\\\\\\"variables\\\\\\\"\\\", \
\\\"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \\!\
\\(\\*StyleBox[\\\"\\\\\\\"parameters\\\\\\\"\\\", \
\\\"TI\\\"]\\)\\!\\(\\*StyleBox[\\\"\\\\\\\"]\\\\\\\"\\\", \\\"MR\\\"]\\) \
\"\>", \("fit the data to the model with the named variables and parameters, \
returning a list of rules expressing the fit parameter estimates and fit \
diagnostics"\)}
                },
              ColumnAlignments->{Right, Left}],
            GridBoxOptions->{RowLines->False}]],
        "2ColumnBox"], TraditionalForm]], "Text",
  CellTags->{"S6.5.1", "1.3"}],

Cell[TextData[{
  "The ",
  StyleBox["NonlinearFit", "MR"],
  " and ",
  StyleBox["NonlinearRegress", "MR"],
  " functions. "
}], "Caption",
  CellTags->"S6.5.1"],

Cell[TextData[{
  "The model argument of ",
  StyleBox["NonlinearFit", "MR"],
  " and ",
  StyleBox["NonlinearRegress", "MR"],
  " must be completely specified by the symbols in the variables argument and \
the symbols in the parameters argument. The variables argument specifies the \
independent variables represented in the data. The parameters argument \
specifies the model parameters for which you would like estimates. \nThe data \
argument can be a list of vectors, each vector consisting of the values of \
the independent variables ",
  Cell[BoxData[
      \(TraditionalForm\`x\_1\)], "InlineFormula"],
  ", ",
  Cell[BoxData[
      \(TraditionalForm\`x\_2\)], "InlineFormula"],
  ", \[Ellipsis] , followed by the observed value of the associated response \
",
  Cell[BoxData[
      \(TraditionalForm\`y\)], "InlineFormula"],
  ". The data argument can also be a vector, in which case it is assumed that \
the vector represents the observed values of the response variable, with the \
independent variable(s) for the ",
  Cell[BoxData[
      \(TraditionalForm\`i\^th\)], "InlineFormula"],
  " response equal to ",
  Cell[BoxData[
      \(TraditionalForm\`i\)], "InlineFormula"],
  ". "
}], "Text",
  CellTags->{"S6.5.1", "1.4"}],

Cell[BoxData[
    FormBox[
      StyleBox[
        FrameBox[
          StyleBox[GridBox[{
                {
                  RowBox[{
                    StyleBox["{",
                      "InlineFormula"], 
                    StyleBox[
                      
                      SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"y\\\\\\\
\"\\\", \\\"TI\\\"]\\)\"\>", "1"],
                      
                      "InlineFormula"], \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\)\"\>", 
                    StyleBox[
                      
                      SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"y\\\\\\\
\"\\\", \\\"TI\\\"]\\)\"\>", "2"],
                      
                      "InlineFormula"], \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \[Ellipsis] \
\"\>", 
                    StyleBox["}",
                      "InlineFormula"], "\<\" \"\>"}], 
                  
                  RowBox[{"\<\"data points specified by a list of response \
values, where \"\>", 
                    StyleBox["i",
                      
                      "InlineFormula"], "\<\" is the value of the independent \
variable associated with the \"\>", 
                    StyleBox["i",
                      "InlineFormula"], 
                    StyleBox[\(\[Null]\^"th"\),
                      "InlineFormula"], "\<\" response, \"\>", 
                    StyleBox[\(y\_i\),
                      "InlineFormula"], "\<\" \"\>"}]},
                {
                  RowBox[{
                    StyleBox["{",
                      "InlineFormula"], 
                    StyleBox["{",
                      "InlineFormula"], 
                    StyleBox[
                      
                      SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"x\\\\\\\
\"\\\", \\\"TI\\\"]\\)\"\>", "11"],
                      
                      "InlineFormula"], \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \"\>", 
                    StyleBox[
                      
                      SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"x\\\\\\\
\"\\\", \\\"TI\\\"]\\)\"\>", "12"],
                      
                      "InlineFormula"], \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \[Ellipsis] \
\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \"\>", 
                    StyleBox[
                      
                      SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"y\\\\\\\
\"\\\", \\\"TI\\\"]\\)\"\>", "1"],
                      "InlineFormula"], 
                    StyleBox["}",
                      
                      "InlineFormula"], \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\)\"\>", 
                    StyleBox["{",
                      "InlineFormula"], 
                    StyleBox[
                      
                      SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"x\\\\\\\
\"\\\", \\\"TI\\\"]\\)\"\>", "21"],
                      
                      "InlineFormula"], \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \"\>", 
                    StyleBox[
                      
                      SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"x\\\\\\\
\"\\\", \\\"TI\\\"]\\)\"\>", "22"],
                      
                      "InlineFormula"], \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \[Ellipsis] \
\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \"\>", 
                    StyleBox[
                      
                      SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"y\\\\\\\
\"\\\", \\\"TI\\\"]\\)\"\>", "2"],
                      "InlineFormula"], 
                    StyleBox["}",
                      
                      "InlineFormula"], \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \[Ellipsis] \
\"\>", 
                    StyleBox["}",
                      "InlineFormula"], "\<\" \"\>"}], 
                  
                  RowBox[{"\<\"data points specified by a matrix, where \"\>",
                     
                    StyleBox[\(x\_\(i  k\)\),
                      "InlineFormula"], "\<\" is the value of the \"\>", 
                    StyleBox["i",
                      "InlineFormula"], 
                    StyleBox[\(\[Null]\^"th"\),
                      "InlineFormula"], "\<\" case of the \"\>", 
                    StyleBox["k",
                      "InlineFormula"], 
                    StyleBox[\(\[Null]\^"th"\),
                      "InlineFormula"], "\<\" independent variable, and \"\>",
                     
                    StyleBox[\(y\_i\),
                      "InlineFormula"], "\<\" is the \"\>", 
                    StyleBox["i",
                      "InlineFormula"], 
                    StyleBox[\(\[Null]\^"th"\),
                      "InlineFormula"], "\<\" response \"\>"}]}
                },
              ColumnAlignments->{Right, Left}],
            GridBoxOptions->{RowLines->False}]],
        "2ColumnBox"], TraditionalForm]], "Text",
  CellTags->{"S6.5.1", "1.5"}],

Cell[TextData[{
  "Ways of specifying data in ",
  StyleBox["NonlinearFit", "MR"],
  " and ",
  StyleBox["NonlinearRegress", "MR"],
  "."
}], "Caption",
  CellTags->"S6.5.1"],

Cell[TextData[{
  "The estimates of the model parameters are chosen to minimize the ",
  Cell[BoxData[
      \(TraditionalForm\`\[Chi]\^2\)], "InlineFormula"],
  " merit function given by the sum of squared residuals ",
  Cell[BoxData[
      FormBox[
        RowBox[{
          UnderscriptBox["\[Sum]", "i",
            LimitsPositioning->True], \(e\_i\%2\)}], TraditionalForm]], 
    "InlineFormula"],
  ". The optimization methods used by ",
  StyleBox["NonlinearFit", "MR"],
  " and ",
  StyleBox["NonlinearRegress", "MR"],
  " are iterative so starting values are required for the parameter estimate \
search. Careful choice of starting values may be necessary, as the parameter \
estimates found by ",
  StyleBox["NonlinearFit", "MR"],
  " and ",
  StyleBox["NonlinearRegress", "MR"],
  " may represent a local minimum in the ",
  Cell[BoxData[
      \(TraditionalForm\`\[Chi]\^2\)], "InlineFormula"],
  " merit function. \nThe simplest way to specify a parameter is as a symbol, \
assuming a starting point of 1.0, or as a ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  StyleBox["symbol", "TI"],
  StyleBox[",", "MR"],
  " ",
  StyleBox["start", "TI"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  " pair. A parameter can be specified to lie within a prescribed range using \
",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  StyleBox["symbol", "TI"],
  StyleBox[",", "MR"],
  " ",
  StyleBox["start", "TI"],
  StyleBox[",", "MR"],
  " ",
  StyleBox["min", "TI"],
  StyleBox[",", "MR"],
  " ",
  StyleBox["max", "TI"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  " or ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  StyleBox["symbol", "TI"],
  StyleBox[",", "MR"],
  " ",
  StyleBox["min", "TI"],
  StyleBox[",", "MR"],
  " ",
  StyleBox["max", "TI"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  ". When elements in the parameter are specified as ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  StyleBox["symbol", "TI"],
  StyleBox[",", "MR"],
  " ",
  StyleBox["min", "TI"],
  StyleBox[",", "MR"],
  " ",
  StyleBox["max", "TI"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  ", the starting parameter values are taken to be those minimizing the ",
  Cell[BoxData[
      \(TraditionalForm\`\[Chi]\^2\)], "InlineFormula"],
  " merit function out of the set forming a ",
  Cell[BoxData[
      \(TraditionalForm\`2\^p\)], "InlineFormula"],
  " factorial design based on the parameter ranges, where ",
  Cell[BoxData[
      \(TraditionalForm\`p\)], "InlineFormula"],
  " is the number of parameters. For example, if the parameter list is \
specified as ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  StyleBox["a", "TI"],
  StyleBox[",", "MR"],
  " ",
  StyleBox["0,", "MR"],
  " ",
  StyleBox["1", "MR"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  StyleBox[",", "MR"],
  " ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  StyleBox["b", "TI"],
  StyleBox[",", "MR"],
  " ",
  StyleBox["0,", "MR"],
  " ",
  StyleBox["3", "MR"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  ", then the value of ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  StyleBox["a", "TI"],
  StyleBox[",", "MR"],
  " ",
  StyleBox["b", "TI"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  " in ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  StyleBox["1/3,", "MR"],
  " ",
  StyleBox["1", "MR"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  StyleBox[",", "MR"],
  " ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  StyleBox["1/3,", "MR"],
  " ",
  StyleBox["2", "MR"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  StyleBox[",", "MR"],
  " ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  StyleBox["2/3,", "MR"],
  " ",
  StyleBox["1", "MR"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  StyleBox[",", "MR"],
  " ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  StyleBox["2/3,", "MR"],
  " ",
  StyleBox["2", "MR"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  " that yields the minimum ",
  Cell[BoxData[
      \(TraditionalForm\`\[Chi]\^2\)], "InlineFormula"],
  " gives starting values for parameters ",
  Cell[BoxData[
      \(TraditionalForm\`a\)], "InlineFormula"],
  " and ",
  Cell[BoxData[
      \(TraditionalForm\`b\)], "InlineFormula"],
  ". When the parameter range is given, the search for the parameter \
terminates if the search goes outside the specified range. \nIf a parameter \
is specified as ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  StyleBox["symbol", "TI"],
  StyleBox[",", "MR"],
  " ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  Cell[BoxData[
      FormBox[
        SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"start\\\\\\\"\\\", \\\
\"TI\\\"]\\)\"\>", "0"], TraditionalForm]], "InlineFormula"],
  StyleBox[",", "MR"],
  " ",
  Cell[BoxData[
      FormBox[
        SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"start\\\\\\\"\\\", \\\
\"TI\\\"]\\)\"\>", "1"], TraditionalForm]], "InlineFormula"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  " or ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  StyleBox["symbol", "TI"],
  StyleBox[",", "MR"],
  " ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  Cell[BoxData[
      FormBox[
        SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"start\\\\\\\"\\\", \\\
\"TI\\\"]\\)\"\>", "0"], TraditionalForm]], "InlineFormula"],
  StyleBox[",", "MR"],
  " ",
  Cell[BoxData[
      FormBox[
        SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"start\\\\\\\"\\\", \\\
\"TI\\\"]\\)\"\>", "1"], TraditionalForm]], "InlineFormula"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  StyleBox[",", "MR"],
  " ",
  StyleBox["min", "TI"],
  StyleBox[",", "MR"],
  " ",
  StyleBox["max", "TI"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  ", then the search for parameter estimates uses ",
  Cell[BoxData[
      \(TraditionalForm\`s  t  a  r  t\_0\)], "InlineFormula"],
  " and ",
  Cell[BoxData[
      \(TraditionalForm\`s  t  a  r  t\_1\)], "InlineFormula"],
  " as the first two values of ",
  Cell[BoxData[
      \(TraditionalForm\`s  y  m  b  o  l\)], "InlineFormula"],
  ". This form must be used if symbolic derivatives of ",
  Cell[BoxData[
      \(TraditionalForm\`\[Chi]\^2\)], "InlineFormula"],
  " with respect to the parameters cannot be found. "
}], "Text",
  CellTags->{"S6.5.1", "1.6"}],

Cell["This loads the package. ", "MathCaption",
  CellTags->{"S6.5.1", "1.7"}],

Cell["<< Statistics`NonlinearFit`", "Input",
  CellLabel->"In[1]:= ",
  CellTags->"S6.5.1"],

Cell[TextData[{
  "This data set [Meyer and Roth, 1972] gives five values for ",
  Cell[BoxData[
      \(TraditionalForm\`x\_1\)], "InlineFormula"],
  ", ",
  Cell[BoxData[
      \(TraditionalForm\`x\_2\)], "InlineFormula"],
  ", and ",
  Cell[BoxData[
      \(TraditionalForm\`y\)], "InlineFormula"],
  ", describing a reaction involving the catalytic dehydration of \
n\[Hyphen]hexyl alcohol. Here ",
  Cell[BoxData[
      \(TraditionalForm\`x\_1\)], "InlineFormula"],
  " is partial pressure of alchohol, ",
  Cell[BoxData[
      \(TraditionalForm\`x\_2\)], "InlineFormula"],
  " is olefin, and ",
  Cell[BoxData[
      \(TraditionalForm\`y\)], "InlineFormula"],
  " is the rate of reaction. "
}], "MathCaption",
  CellTags->{"S6.5.1", "1.8"}],

Cell["\<\
data = {{1.0, 1.0, .126}, {2.0, 1.0, .219},
{1.0, 2.0, .076}, {2.0, 2.0, .126}, {.1, .0, .186}};\
\>", "Input",
  CellLabel->"In[2]:= ",
  CellTags->"S6.5.1"],

Cell[TextData[{
  "This finds the values of the parameters ",
  Cell[BoxData[
      \(TraditionalForm\`\[Theta]\_1\)], "InlineFormula"],
  ", ",
  Cell[BoxData[
      \(TraditionalForm\`\[Theta]\_2\)], "InlineFormula"],
  ", and ",
  Cell[BoxData[
      \(TraditionalForm\`\[Theta]\_3\)], "InlineFormula"],
  ", such that the model ",
  Cell[BoxData[
      \(TraditionalForm\`y = \(\[Theta]\_1\) \(\[Theta]\_3\) 
        x\_1/\((1 + \(\[Theta]\_1\) x\_1 + \(\[Theta]\_2\) x\_2)\)\)], 
    "InlineFormula"],
  " best fits the ",
  Cell[BoxData[
      \(TraditionalForm\`\((x\_1, x\_2, y)\)\)], "InlineFormula"],
  " data triples in a least\[Hyphen]squares sense. The least\[Hyphen]squares \
parameters are substituted into the model, just as in ",
  StyleBox["Fit", "MR"],
  ". "
}], "MathCaption",
  CellTags->{"S6.5.1", "1.9"}],

Cell[CellGroupData[{

Cell["\<\
NonlinearFit[data,
theta1 theta3 x1 / (1 + theta1 x1 + theta2 x2),
{x1, x2}, {theta1, theta2, theta3}]\
\>", "Input",
  CellLabel->"In[3]:= ",
  CellTags->"S6.5.1"],

Cell[BoxData[
    \(\(2.442770447552637`\ x1\)\/\(1 + 3.1315068139503124`\ x1 + \
15.159362262745494`\ x2\)\)], "Output",
  CellLabel->"Out[3]= ",
  CellTags->"S6.5.1"]
}, Open  ]],

Cell[TextData[{
  StyleBox["NonlinearRegress", "MR"],
  " can be used to produce a list of replacement rules for the parameters. "
}], "MathCaption",
  CellTags->{"S6.5.1", "1.10"}],

Cell[CellGroupData[{

Cell["\<\
BestFitParameters /. NonlinearRegress[data,
theta1 theta3 x1 / (1 + theta1 x1 + theta2 x2),
{x1, x2}, {theta1, theta2, theta3},
RegressionReport -> BestFitParameters]\
\>", "Input",
  CellLabel->"In[4]:= ",
  CellTags->"S6.5.1"],

Cell[BoxData[
    \({theta1 \[Rule] 3.1315068139503124`, 
      theta2 \[Rule] 15.159362262745494`, 
      theta3 \[Rule] 0.7800623127085414`}\)], "Output",
  CellLabel->"Out[4]= ",
  CellTags->"S6.5.1"]
}, Open  ]],

Cell[BoxData[
    FormBox[
      StyleBox[
        FrameBox[
          StyleBox[GridBox[{
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"option name\\\\\\\"\\\", \
\\\"TI\\\"]\\) \"\>", "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"default \
value\\\\\\\"\\\", \\\"TI\\\"]\\) \"\>", "\[Null]"},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"AccuracyGoal\\\\\\\"\\\", \
\\\"MR\\\"]\\) \"\>", \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Automatic\\\\\\\"\\\", \\\"MR\\\"]\\) \
\"\>", 
                  RowBox[{"\<\"accuracy desired for the \"\>", 
                    StyleBox[\(\[Chi]\^2\),
                      "InlineFormula"], "\<\" merit function\"\>"}]},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Compiled\\\\\\\"\\\", \
\\\"MR\\\"]\\) \"\>", "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"True\\\\\\\"\\\", \\\
\"MR\\\"]\\) \"\>", \("whether the merit function should be compiled "\)},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Gradient\\\\\\\"\\\", \
\\\"MR\\\"]\\) \"\>", \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Automatic\\\\\\\"\\\", \\\"MR\\\"]\\) \
\"\>", "\<\"list of gradient functions\"\>"},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"MaxIterations\\\\\\\"\\\",\
 \\\"MR\\\"]\\) \"\>", "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"30\\\\\\\"\\\", \
\\\"MR\\\"]\\) \"\>", \("maximum number of iterations in search"\)},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Method\\\\\\\"\\\", \
\\\"MR\\\"]\\) \"\>", "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"LevenbergMarquardt\\\
\\\\\"\\\", \\\"MR\\\"]\\) \"\>", "\<\"method used in fit\"\>"},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"PrecisionGoal\\\\\\\"\\\",\
 \\\"MR\\\"]\\) \"\>", "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Automatic\\\\\\\"\\\
\", \\\"MR\\\"]\\) \"\>", 
                  RowBox[{"\<\"precision desired for the \"\>", 
                    StyleBox[\(\[Chi]\^2\),
                      "InlineFormula"], "\<\" merit function\"\>"}]},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Tolerance\\\\\\\"\\\", \\\
\"MR\\\"]\\) \"\>", "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Automatic\\\\\\\"\\\", \
\\\"MR\\\"]\\) \"\>", \("numerical tolerance for matrix operations"\)},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Weights\\\\\\\"\\\", \
\\\"MR\\\"]\\) \"\>", \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"Automatic\\\\\\\"\\\", \\\"MR\\\"]\\) \
\"\>", \("list of weights for each data point or pure function"\)},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"WorkingPrecision\\\\\\\"\\\
\", \\\"MR\\\"]\\) \"\>", \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"$MachinePrecision\\\\\\\"\\\", \\\"MR\\\"]\
\\) \"\>", \("precision of the arithmetic to use in calculations"\)}
                },
              ColumnAlignments->{Left}],
            GridBoxOptions->{RowLines->{True, False}}]],
        "3ColumnBox"], TraditionalForm]], "Text",
  CellTags->{"S6.5.1", "1.11"}],

Cell[TextData[{
  "Options for ",
  StyleBox["NonlinearFit", "MR"],
  " and ",
  StyleBox["NonlinearRegress", "MR"],
  "."
}], "Caption",
  CellTags->"S6.5.1"],

Cell[TextData[{
  StyleBox["NonlinearFit", "MR"],
  " and ",
  StyleBox["NonlinearRegress", "MR"],
  " use ",
  StyleBox["FindMinimum", "MR"],
  " to find parameter estimates, thus both functions accept the ",
  StyleBox["FindMinimum", "MR"],
  " options ",
  StyleBox["AccuracyGoal", "MR"],
  ", ",
  StyleBox["Compiled", "MR"],
  ", ",
  StyleBox["Gradient", "MR"],
  ", ",
  StyleBox["MaxIterations", "MR"],
  ", ",
  StyleBox["Method", "MR"],
  ", ",
  StyleBox["PrecisionGoal", "MR"],
  ", and ",
  StyleBox["WorkingPrecision", "MR"],
  ". The ",
  StyleBox["Method", "MR"],
  " option allows you to choose between algorithms for performing the \
minimization of the ",
  Cell[BoxData[
      \(TraditionalForm\`\[Chi]\^2\)], "InlineFormula"],
  " merit function. The ",
  StyleBox["LevenbergMarquardt", "MR"],
  " method gradually shifts the search for the minimum of ",
  Cell[BoxData[
      \(TraditionalForm\`\[Chi]\^2\)], "InlineFormula"],
  " from steepest descent to quadratic minimization. Other possible settings \
for ",
  StyleBox["Method", "MR"],
  " are ",
  StyleBox["Gradient", "MR"],
  " (steepest descent), ",
  StyleBox["Newton", "MR"],
  ", ",
  StyleBox["QuasiNewton", "MR"],
  ", and ",
  StyleBox["Automatic", "MR"],
  ". The ",
  StyleBox["Automatic", "MR"],
  " method does linear fitting for linear models and ",
  StyleBox["LevenbergMarquardt", "MR"],
  " search for nonlinear models, thus the meaning of the automatic method \
differs between ",
  StyleBox["NonlinearFit", "MR"],
  " (or ",
  StyleBox["NonlinearRegress", "MR"],
  ") and ",
  StyleBox["FindMinimum", "MR"],
  ". \nThe ",
  StyleBox["Weights", "MR"],
  " option allows you to implement weighted least squares by specifying a \
list of weights, one for each data point; the default ",
  StyleBox["Weights", "MR"],
  " ",
  StyleBox["->", "MR"],
  " ",
  StyleBox["Automatic", "MR"],
  " implies a weight of unity for each data point. When ",
  StyleBox["Weights", "MR"],
  " ",
  StyleBox["->", "MR"],
  " ",
  Cell[BoxData[
      \(TraditionalForm\`{\)], "InlineFormula"],
  Cell[BoxData[
      FormBox[
        SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"w\\\\\\\"\\\", \\\"TI\
\\\"]\\)\"\>", "1"], TraditionalForm]], "InlineFormula"],
  StyleBox[",", "MR"],
  " \[Ellipsis] ",
  StyleBox[",", "MR"],
  " ",
  Cell[BoxData[
      FormBox[
        SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"w\\\\\\\"\\\", \\\"TI\
\\\"]\\)\"\>", "n"], TraditionalForm]], "InlineFormula"],
  Cell[BoxData[
      \(TraditionalForm\`}\)], "InlineFormula"],
  ", the parameter estimates are chosen to minimize the weighted sum of \
squared residuals ",
  Cell[BoxData[
      FormBox[
        RowBox[{
          UnderoverscriptBox["\[Sum]", \(i = 1\), "n",
            LimitsPositioning->True], \(w\_i\), " ", \(e\_i\%2\)}], 
        TraditionalForm]], "InlineFormula"],
  ". \n",
  StyleBox["Weights", "MR"],
  " can also specify a pure function of the response. For example, to choose \
parameter estimates to minimize ",
  Cell[BoxData[
      FormBox[
        RowBox[{
          UnderoverscriptBox["\[Sum]", \(i = 1\), "n",
            LimitsPositioning->True], \(\@y\_i\), " ", \(e\_i\%2\)}], 
        TraditionalForm]], "InlineFormula"],
  ", set ",
  StyleBox["Weights", "MR"],
  " ",
  StyleBox["->", "MR"],
  " ",
  StyleBox["(Sqrt[#]", "MR"],
  " ",
  StyleBox["&)", "MR"],
  ". "
}], "Text",
  CellTags->{"S6.5.1", "1.12"}],

Cell[BoxData[
    FormBox[
      StyleBox[
        FrameBox[
          StyleBox[GridBox[{
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"option name\\\\\\\"\\\", \
\\\"TI\\\"]\\) \"\>", "\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"default \
value\\\\\\\"\\\", \\\"TI\\\"]\\) \"\>", "\[Null]"},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"RegressionReport\\\\\\\"\\\
\", \\\"MR\\\"]\\) \"\>", \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"SummaryReport\\\\\\\"\\\", \\\"MR\\\"]\\) \
\"\>", \("objects to be included in output"\)}
                },
              ColumnAlignments->{Left}],
            GridBoxOptions->{RowLines->True}]],
        "3ColumnBox"], TraditionalForm]], "Text",
  CellTags->{"S6.5.1", "1.13"}],

Cell[TextData[{
  "Options unique to ",
  StyleBox["NonlinearRegress", "MR"],
  "."
}], "Caption",
  CellTags->"S6.5.1"],

Cell[TextData[{
  "The option ",
  StyleBox["RegressionReport", "MR"],
  " is accepted by ",
  StyleBox["NonlinearRegress", "MR"],
  ", but not ",
  StyleBox["NonlinearFit", "MR"],
  ". If ",
  StyleBox["RegressionReport", "MR"],
  " is not specified, ",
  StyleBox["NonlinearRegress", "MR"],
  " automatically gives a list including values for ",
  StyleBox["BestFitParameters", "MR"],
  ", ",
  StyleBox["ParameterCITable", "MR"],
  ", ",
  StyleBox["EstimatedVariance", "MR"],
  ", ",
  StyleBox["ANOVATable", "MR"],
  ", ",
  StyleBox["AsymptoticCorrelationMatrix", "MR"],
  ", and ",
  StyleBox["FitCurvatureTable", "MR"],
  ". This set of objects comprises the default ",
  StyleBox["SummaryReport", "MR"],
  ". The option ",
  StyleBox["RegressionReport", "MR"],
  " can be used to specify a single object or a list of objects so that more \
(or less) than the default set of results is included in the output. ",
  StyleBox["RegressionReportValues[NonlinearRegress]", "MR"],
  " gives the objects that may be included in the ",
  StyleBox["RegressionReport", "MR"],
  " list for the ",
  StyleBox["NonlinearRegress", "MR"],
  " function. \n",
  StyleBox["NonlinearRegress", "MR"],
  " will also accept any option that can be specified for ",
  StyleBox["StudentTCI", "MR"],
  ". In particular, ",
  StyleBox["ConfidenceLevel", "MR"],
  " specifies the reference parameter confidence region (the relative \
curvature of which is given in ",
  StyleBox["FitCurvatureTable", "MR"],
  ") and the individual parameter confidence intervals given in ",
  StyleBox["ParameterCITable", "MR"],
  ". "
}], "Text",
  CellTags->{"S6.5.1", "1.14"}],

Cell[TextData[{
  "If the ",
  StyleBox["RegressionReport", "MR"],
  " option is not specified, the default regression output is produced. Note \
that the ",
  StyleBox["AsymptoticCorrelationMatrix", "MR"],
  " indicates that the parameters ",
  Cell[BoxData[
      \(TraditionalForm\`\[Theta]\_1\)], "InlineFormula"],
  " and ",
  Cell[BoxData[
      \(TraditionalForm\`\[Theta]\_3\)], "InlineFormula"],
  " are highly correlated, a negative aspect of this model. "
}], "MathCaption",
  CellTags->{"S6.5.1", "1.15"}],

Cell[CellGroupData[{

Cell["\<\
NonlinearRegress[data,
theta1 theta3 x1 / (1 + theta1 x1 + theta2 x2),
{x1, x2}, {theta1, theta2, theta3}]\
\>", "Input",
  CellLabel->"In[5]:= ",
  CellTags->"S6.5.1"],

Cell[BoxData[
    RowBox[{"{", 
      RowBox[{\(BestFitParameters \[Rule] {theta1 \[Rule] 
              3.1315068139503124`, theta2 \[Rule] 15.159362262745494`, 
            theta3 \[Rule] 0.7800623127085414`}\), ",", 
        RowBox[{"ParameterCITable", "\[Rule]", 
          TagBox[GridBox[{
                {"\<\"\"\>", "\<\"Estimate\"\>", "\<\"Asymptotic SE\"\>", "\<\
\"CI\"\>"},
                {"theta1", "3.1315068139503124`", 
                  "0.8084171841304373`", \({\(-0.3468315901248884`\), 
                    6.609845218025513`}\)},
                {"theta2", "15.159362262745494`", 
                  "0.6312241558637214`", \({12.443423925434653`, 
                    17.875300600056335`}\)},
                {"theta3", "0.7800623127085414`", 
                  "0.15179702167157247`", \({0.12693244304551166`, 
                    1.4331921823715712`}\)}
                },
              RowSpacings->1,
              ColumnSpacings->3,
              RowAlignments->Baseline,
              ColumnAlignments->{Left}],
            (
            TableForm[ #, TableDepth -> 2, TableHeadings -> {{theta1, theta2, 
              theta3}, {"Estimate", "Asymptotic SE", "CI"}}]&)]}], 
        ",", \(EstimatedVariance \[Rule] 0.000021776330970991434`\), ",", 
        RowBox[{"ANOVATable", "\[Rule]", 
          InterpretationBox[GridBox[{
                {"\<\"\"\>", "\<\"DF\"\>", "\<\"SumOfSq\"\>", \
"\<\"MeanSq\"\>"},
                {"\<\"Model\"\>", "3", "0.12004144733805802`", 
                  "0.040013815779352674`"},
                {"\<\"Error\"\>", "2", "0.00004355266194198287`", 
                  "0.000021776330970991434`"},
                {"\<\"Uncorrected Total\"\>", "5", 
                  "0.12008500000000001`", "\<\"\"\>"},
                {"\<\"Corrected Total\"\>", "4", 
                  "0.012627199999999998`", "\<\"\"\>"}
                },
              RowSpacings->1,
              ColumnSpacings->3,
              RowAlignments->Baseline,
              ColumnAlignments->{Left}],
            
            TableForm[ {{3, 0.12004144733805802, 0.040013815779352674}, {2, 
              0.000043552661941982868`17, 0.000021776330970991434`17}, {5, 
              0.12008500000000001}, {4, 0.012627199999999998}}, 
              TableHeadings -> {{"Model", "Error", "Uncorrected Total", 
              "Corrected Total"}, {"DF", "SumOfSq", "MeanSq"}}]]}], ",", 
        RowBox[{"AsymptoticCorrelationMatrix", "\[Rule]", 
          TagBox[
            RowBox[{"(", "\[NoBreak]", GridBox[{
                  {"0.9999999999999998`", 
                    "0.2870147850116305`", \(-0.9918299862926232`\)},
                  {"0.2870147850116305`", 
                    "1.0000000000000002`", \(-0.18202091524481895`\)},
                  {\(-0.9918299862926235`\), \(-0.18202091524481895`\), 
                    "1.`"}
                  }], "\[NoBreak]", ")"}],
            (MatrixForm[ #]&)]}], ",", 
        RowBox[{"FitCurvatureTable", "\[Rule]", 
          TagBox[GridBox[{
                {"\<\"\"\>", "\<\"Curvature\"\>"},
                {"\<\"Max Intrinsic\"\>", "0.037466954505343`"},
                {"\<\"Max Parameter-Effects\"\>", "12.816079819177288`"},
                {"\<\"95. % Confidence Region\"\>", "0.22843024672369935`"}
                },
              RowSpacings->1,
              ColumnSpacings->3,
              RowAlignments->Baseline,
              ColumnAlignments->{Left}],
            (
            TableForm[ #, TableDepth -> 2, TableHeadings -> {{"Max Intrinsic",
               "Max Parameter-Effects", "95. % Confidence Region"}, {
              "Curvature"}}]&)]}]}], "}"}]], "Output",
  CellLabel->"Out[5]= ",
  CellTags->"S6.5.1"]
}, Open  ]],

Cell[BoxData[
    FormBox[
      StyleBox[
        FrameBox[
          StyleBox[GridBox[{
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"BestFit\\\\\\\"\\\", \
\\\"MR\\\"]\\) \"\>", "\<\"best fit function \"\>"},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"BestFitParameters\\\\\\\"\
\\\", \\\"MR\\\"]\\) \"\>", "\<\"best fit parameter estimates\"\>"},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"ANOVATable\\\\\\\"\\\", \
\\\"MR\\\"]\\) \"\>", "\<\"analysis of variance table\"\>"},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"EstimatedVariance\\\\\\\"\
\\\", \\\"MR\\\"]\\) \"\>", "\<\"estimated error variance \"\>"},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"ParameterTable\\\\\\\"\\\"\
, \\\"MR\\\"]\\) \"\>", \("table of parameter information including \
asymptotic standard errors and test statistics "\)},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"ParameterCITable\\\\\\\"\\\
\", \\\"MR\\\"]\\) \"\>", \("table of asymptotic confidence intervals for the \
parameters"\)},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"ParameterConfidenceRegion\
\\\\\\\"\\\", \\\"MR\\\"]\\) \"\>", \("asymptotic ellipsoidal joint \
confidence region for the parameters"\)},
                {
                  
                  RowBox[{"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"\
ParameterConfidenceRegion[\\\\\\\"\\\", \\\"MR\\\"]\\)\"\>", 
                    StyleBox["{",
                      "InlineFormula"], 
                    StyleBox[
                      
                      SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"parm\\\\\
\\\"\\\", \\\"TI\\\"]\\)\"\>", \(i  1\)],
                      
                      "InlineFormula"], \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \"\>", 
                    StyleBox[
                      
                      SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"parm\\\\\
\\\"\\\", \\\"TI\\\"]\\)\"\>", \(i  2\)],
                      
                      "InlineFormula"], \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \[Ellipsis] \
\"\>", 
                    StyleBox["}",
                      
                      "InlineFormula"], \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"]\\\\\\\"\\\", \\\"MR\\\"]\\) \"\>"}], 
                  
                  RowBox[{"\<\"asymptotic ellipsoidal conditional joint \
confidence region for the parameters \"\>", 
                    StyleBox["{",
                      "InlineFormula"], 
                    StyleBox[
                      
                      SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"parm\\\\\
\\\"\\\", \\\"TI\\\"]\\)\"\>", \(i  1\)],
                      
                      "InlineFormula"], \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \"\>", 
                    StyleBox[
                      
                      SubscriptBox["\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"parm\\\\\
\\\"\\\", \\\"TI\\\"]\\)\"\>", \(i  2\)],
                      
                      "InlineFormula"], \
"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\",\\\\\\\"\\\", \\\"MR\\\"]\\) \[Ellipsis] \
\"\>", 
                    StyleBox["}",
                      "InlineFormula"]}]},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"FitResiduals\\\\\\\"\\\", \
\\\"MR\\\"]\\) \"\>", \("differences between the observed responses and the \
predicted responses "\)},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"PredictedResponse\\\\\\\"\
\\\", \\\"MR\\\"]\\) \"\>", \("fitted values obtained by evaluating the best \
fit function at the observed values of the independent variables"\)},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"SinglePredictionCITable\\\
\\\\\"\\\", \\\"MR\\\"]\\) \"\>", \("table of asymptotic confidence intervals \
for predicting a single observation of the response variable "\)},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"MeanPredictionCITable\\\\\
\\\"\\\", \\\"MR\\\"]\\) \"\>", \("table of asymptotic confidence intervals \
for predicting the expected value of the response variable"\)},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"\
AsymptoticCovarianceMatrix\\\\\\\"\\\", \\\"MR\\\"]\\) \"\>", "\<\"asymptotic \
covariance matrix of the parameters \"\>"},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"\
AsymptoticCorrelationMatrix\\\\\\\"\\\", \\\"MR\\\"]\\) \"\>", \
"\<\"asymptotic correlation matrix of the parameters \"\>"}
                },
              ColumnAlignments->{Right, Left}],
            GridBoxOptions->{RowLines->False}]],
        "2ColumnBox"], TraditionalForm]], "Text",
  CellTags->{"S6.5.1", "1.16"}],

Cell[TextData[{
  "Some option settings for ",
  StyleBox["RegressionReport", "MR"],
  " or objects that may be included in a list specified by ",
  StyleBox["RegressionReport", "MR"],
  "."
}], "Caption",
  CellTags->"S6.5.1"],

Cell[TextData[{
  StyleBox["ANOVATable", "MR"],
  " provides a partitioning of the sum of squares. Both the uncorrected total \
sum of squares (without intercept term) and the corrected total sum of \
squares (with intercept term) are included. The corrected total sum of \
squares is useful for calculating a value for ",
  Cell[BoxData[
      \(TraditionalForm\`R\^2\)], "InlineFormula"],
  " comparable to that produced by ",
  StyleBox["Regress", "MR"],
  ", so that the fit of a nonlinear model can be compared to the fit of a \
linear model. For example, one possible definition of ",
  Cell[BoxData[
      \(TraditionalForm\`R\^2\)], "InlineFormula"],
  " for a nonlinear regression is the ratio of the difference between the \
corrected total sum of squares and the residual sum of squares to the \
corrected total sum of squares. The other objects in this table have meanings \
similar to what they would have if specified by the ",
  StyleBox["RegressionReport", "MR"],
  " option of the linear regression function ",
  StyleBox["Regress", "MR"],
  ". "
}], "Text",
  CellTags->{"S6.5.1", "1.17"}],

Cell[BoxData[
    FormBox[
      StyleBox[
        FrameBox[
          StyleBox[GridBox[{
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"StartingParameters\\\\\\\"\
\\\", \\\"MR\\\"]\\) \"\>", "\<\"starting values for the parameters\"\>"}
                },
              ColumnAlignments->{Right, Left}],
            GridBoxOptions->{RowLines->False}]],
        "2ColumnBox"], TraditionalForm]], "Text",
  CellTags->{"S6.5.1", "1.18"}],

Cell["Objects describing the search for the least\[Hyphen]squares fit.", \
"Caption",
  CellTags->"S6.5.1"],

Cell[TextData[{
  "Including ",
  StyleBox["StartingParameters", "MR"],
  " in the list of objects requested by ",
  StyleBox["RegressionReport", "MR"],
  " is helpful if the search for a least\[Hyphen]squares fit is unsuccessful \
and it is necessary to restart the search at a new point. "
}], "Text",
  CellTags->{"S6.5.1", "1.19"}],

Cell[BoxData[
    FormBox[
      StyleBox[
        FrameBox[
          StyleBox[GridBox[{
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"HatDiagonal\\\\\\\"\\\", \
\\\"MR\\\"]\\) \"\>", 
                  RowBox[{"\<\"diagonal of the hat matrix \"\>", 
                    StyleBox[\(X \(\((\(X\^T\) X)\)\^\(-1\)\) X\^T\),
                      "InlineFormula"], "\<\", where \"\>", 
                    StyleBox["X",
                      
                      "InlineFormula"], "\<\" is the design matrix for the \
approximate linear model at the least\[Hyphen]squares parameter estimates \
\"\>"}]},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"StandardizedResiduals\\\\\
\\\"\\\", \\\"MR\\\"]\\) \"\>", \("fit residuals scaled by their asymptotic \
standard errors, computed using the estimated error variance"\)}
                },
              ColumnAlignments->{Right, Left}],
            GridBoxOptions->{RowLines->False}]],
        "2ColumnBox"], TraditionalForm]], "Text",
  CellTags->{"S6.5.1", "1.20"}],

Cell["Diagnostics for detecting outliers.", "Caption",
  CellTags->"S6.5.1"],

Cell[TextData[{
  "If the nonlinear model is approximately linear at the \
least\[Hyphen]squares fit, ",
  StyleBox["HatDiagonal", "MR"],
  " and ",
  StyleBox["StandardizedResiduals", "MR"],
  " are useful for detecting outliers in the data. ",
  StyleBox["HatDiagonal", "MR"],
  " gives the leverage of each point, measuring whether each observation of \
the predictor variables is unusual; a leverage of zero indicates no \
influence, while a leverage of one indicates that a degree of freedom has \
been lost to fitting that point. "
}], "Text",
  CellTags->{"S6.5.1", "1.21"}],

Cell[BoxData[
    FormBox[
      StyleBox[
        FrameBox[
          StyleBox[GridBox[{
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"FitCurvatureTable\\\\\\\"\
\\\", \\\"MR\\\"]\\) \"\>", "\<\"maximum relative intrinsic curvature, \
maximum relative parameter\[Hyphen]effects curvature, and relative curvature \
of the least\[Hyphen]squares solution confidence region (where confidence \
level is given by \\!\\(\\*StyleBox[\\\"\\\\\\\"ConfidenceLevel\\\\\\\"\\\", \
\\\"MR\\\"]\\))\"\>"},
                {"\<\"\\!\\(\\*StyleBox[\\\"\\\\\\\"ParameterBias\\\\\\\"\\\",\
 \\\"MR\\\"]\\) \"\>", \("bias in the least\[Hyphen]squares estimates of the \
parameters in a nonlinear model"\)}
                },
              ColumnAlignments->{Right, Left}],
            GridBoxOptions->{RowLines->False}]],
        "2ColumnBox"], TraditionalForm]], "Text",
  CellTags->{"S6.5.1", "1.22"}],

Cell["\<\
Diagnostics for evaluating the validity of a linear model \
approximation to the nonlinear model.\
\>", "Caption",
  CellTags->"S6.5.1"],

Cell[TextData[{
  StyleBox["FitCurvatureTable", "MR"],
  " expresses the relative curvature of the solution locus at the least\
\[Hyphen]squares estimate in terms of two components, intrinsic and parameter\
\[Hyphen]effects. For the combination of data and model given above, ",
  Cell[BoxData[
      \(TraditionalForm\`p = 3\)], "InlineFormula"],
  " and ",
  Cell[BoxData[
      \(TraditionalForm\`n = 5\)], "InlineFormula"],
  ". The solution locus is a ",
  Cell[BoxData[
      \(TraditionalForm\`3\)], "InlineFormula"],
  "\[Hyphen]dimensional subspace of ",
  Cell[BoxData[
      \(TraditionalForm\`R\_5\)], "InlineFormula"],
  " given by ",
  Cell[BoxData[
      \(TraditionalForm\`{\(\(1.0\)\(\ \)\) \(\[Theta]\_1\) \[Theta]\_3/\((1 \
+ \(\(1.0\)\(\ \)\) \[Theta]\_1 + \(\(1.0\)\(\ \)\) \[Theta]\_2)\), \)], 
    "InlineFormula"],
  " ",
  Cell[BoxData[
      \(TraditionalForm\`\(\(2.0\)\(\ \)\) \(\[Theta]\_1\) \[Theta]\_3/\((1 + \
\(\(2.0\)\(\ \)\) \[Theta]\_1 + \(\(1.0\)\(\ \)\) \[Theta]\_2)\), \)], 
    "InlineFormula"],
  " ",
  Cell[BoxData[
      \(TraditionalForm\`\(\(1.0\)\(\ \)\) \(\[Theta]\_1\) \[Theta]\_3/\((1 + \
\(\(1.0\)\(\ \)\) \[Theta]\_1 + \(\(2.0\)\(\ \)\) \[Theta]\_2)\), \)], 
    "InlineFormula"],
  " ",
  Cell[BoxData[
      \(TraditionalForm\`\(\(2.0\)\(\ \)\) \(\[Theta]\_1\) \[Theta]\_3/\((1 + \
\(\(2.0\)\(\ \)\) \[Theta]\_1 + \(\(2.0\)\(\ \)\) \[Theta]\_2)\), \)], 
    "InlineFormula"],
  " ",
  Cell[BoxData[
      \(TraditionalForm\`\(\(0.1\)\(\ \)\) \(\[Theta]\_1\) \[Theta]\_3/\((1 + \
\(\(0.1\)\(\ \)\) \[Theta]\_1)\)}\)], "InlineFormula"],
  ", ",
  Cell[BoxData[
      \(TraditionalForm\`{\[Theta]\_1, \[Theta]\_2, \[Theta]\_3} \[Element] 
        R\_3\)], "InlineFormula"],
  ". \nIntrinsic curvature describes the normal component of the solution \
locus curvature at the least\[Hyphen]squares estimate. \
Parameter\[Hyphen]effects curvature describes the tangential component of the \
solution locus curvature at the least\[Hyphen]squares estimate. Standardizing \
curvature to be response\[Hyphen]invariant gives relative curvature. \
Maximizing the relative intrinsic curvature over the ",
  Cell[BoxData[
      \(TraditionalForm\`\((n - p)\)\)], "InlineFormula"],
  "\[Hyphen]dimensional subspace normal to the locus gives the maximum \
relative intrinsic curvature. Similarly, maximizing the relative parameter\
\[Hyphen]effects curvature over the ",
  Cell[BoxData[
      \(TraditionalForm\`p\)], "InlineFormula"],
  "\[Hyphen]dimensional subspace tangential to the locus gives the maximum \
relative parameter\[Hyphen]effects curvature. \nBoth of these quantities can \
be compared to the relative curvature of the confidence region centered on \
the least\[Hyphen]squares parameter estimates. If the maximum relative \
intrinsic curvature is small compared to the confidence region relative \
curvature, the solution locus is approximately planar over the confidence \
region. If the maximum relative parameter\[Hyphen]effects curvature is small \
compared to the confidence region relative curvature, the parameter \
coordinates projected onto the tangential plane are approximately parallel \
and uniformly spaced over the confidence region. \n",
  StyleBox["ParameterBias", "MR"],
  " is based on the average curvature of the solution locus tangential to the \
least\[Hyphen]squares estimate. "
}], "Text",
  CellTags->{"S6.5.1", "1.23"}],

Cell[TextData[{
  "Here the same data is fit to a reparametrized model, ",
  Cell[BoxData[
      \(TraditionalForm\`y = 
        x\_1/\((\[Phi]\_0 + \(\[Phi]\_1\) x\_1 + \(\[Phi]\_2\) x\_2)\) . \)], 
    "InlineFormula"],
  " Only the ",
  StyleBox["BestFit", "MR"],
  ", ",
  StyleBox["AsymptoticCorrelationMatrix", "MR"],
  ", and ",
  StyleBox["FitCurvatureTable", "MR"],
  " are calculated. "
}], "MathCaption",
  CellTags->{"S6.5.1", "1.24"}],

Cell[CellGroupData[{

Cell["\<\
NonlinearRegress[data,
x1 / (phi0 + phi1 x1 + phi2 x2),
{x1, x2}, {phi0, phi1, phi2},
RegressionReport -> {BestFit,
 AsymptoticCorrelationMatrix, FitCurvatureTable}]\
\>", "Input",
  CellLabel->"In[6]:= ",
  CellTags->"S6.5.1"],

Cell[BoxData[
    RowBox[{"{", 
      RowBox[{\(BestFit \[Rule] 
          x1\/\(\(\(0.40937121136690974`\)\(\[InvisibleSpace]\)\) + \
1.2819492885451063`\ x1 + 6.205806546962624`\ x2\)\), ",", 
        RowBox[{"AsymptoticCorrelationMatrix", "\[Rule]", 
          TagBox[
            RowBox[{"(", "\[NoBreak]", GridBox[{
                  {"1.0000000000000002`", \(-0.881408058953822`\), 
                    "0.8026386059232623`"},
                  {\(-0.8814080589538219`\), 
                    "1.`", \(-0.9274654763586443`\)},
                  {"0.8026386059232626`", \(-0.9274654763586443`\), "1.`"}
                  }], "\[NoBreak]", ")"}],
            (MatrixForm[ #]&)]}], ",", 
        RowBox[{"FitCurvatureTable", "\[Rule]", 
          TagBox[GridBox[{
                {"\<\"\"\>", "\<\"Curvature\"\>"},
                {"\<\"Max Intrinsic\"\>", "0.0374669152092768`"},
                {"\<\"Max Parameter-Effects\"\>", "0.089320523274931`"},
                {"\<\"95. % Confidence Region\"\>", "0.22843024672369935`"}
                },
              RowSpacings->1,
              ColumnSpacings->3,
              RowAlignments->Baseline,
              ColumnAlignments->{Left}],
            (
            TableForm[ #, TableDepth -> 2, TableHeadings -> {{"Max Intrinsic",
               "Max Parameter-Effects", "95. % Confidence Region"}, {
              "Curvature"}}]&)]}]}], "}"}]], "Output",
  CellLabel->"Out[6]= ",
  CellTags->"S6.5.1"]
}, Open  ]],

Cell[TextData[{
  "Note that reparametrizing the model does not affect the intrinsic \
curvature. However, here reparametrization has reduced the \
parameter\[Hyphen]effects curvature below the critical value of 0.22843. This \
indicates that the least\[Hyphen]squares estimates of the parameters ",
  Cell[BoxData[
      \(TraditionalForm\`\[Phi]\_0\)], "InlineFormula"],
  ", ",
  Cell[BoxData[
      \(TraditionalForm\`\[Phi]\_1\)], "InlineFormula"],
  ", and ",
  Cell[BoxData[
      \(TraditionalForm\`\[Phi]\_2\)], "InlineFormula"],
  " have nearly linear behavior. This could not be said about the least\
\[Hyphen]squares estimates of the parameters ",
  Cell[BoxData[
      \(TraditionalForm\`\[Theta]\_1\)], "InlineFormula"],
  ", ",
  Cell[BoxData[
      \(TraditionalForm\`\[Theta]\_2\)], "InlineFormula"],
  ", and ",
  Cell[BoxData[
      \(TraditionalForm\`\[Theta]\_3\)], "InlineFormula"],
  " in the original model [Ratkowsky, 1983]. Unfortunately, now the \
asymptotic correlation between parameters is high for all pairs; a model with \
fewer parameters should be considered. \nThe maximum relative curvature \
diagnostics are useful if you wish to make inferences based on the linear \
approximation to the nonlinear model. For example, diagnostics based on the \
asymptotic standard errors of the parameters, such as ",
  StyleBox["ParameterCITable", "MR"],
  " and ",
  StyleBox["ParameterConfidenceRegion", "MR"],
  ", are of questionable validity when the linear approximation to the model \
is poor at the least\[Hyphen]squares parameter estimates. The linear \
approximation is based on the planar assumption, tested by maximum relative \
intrinsic curvature, and the uniform coordinate assumption, tested by maximum \
relative parameter\[Hyphen]effects curvature. It is usually the uniform \
coordinate assumption that is invalid [Bates and Watts, 1988]. "
}], "Text",
  CellTags->{"S6.5.1", "1.25"}],

Cell["References", "Subsubsection",
  CellTags->{"S6.5.1", "1.26"}],

Cell[TextData[{
  "Bates, D. M. and Watts, D. G. (1988). ",
  StyleBox["Nonlinear Regression Analysis and Its Applications", "TI"],
  ", John Wiley & Sons, New York. \nMeyer, R. R. and Roth, P. M. (1972). \
\[OpenCurlyDoubleQuote]Modified Damped Least Squares: An Algorithm for \
Nonlinear Estimation,\[CloseCurlyDoubleQuote] ",
  StyleBox["J. Inst. Math. Appl.", "TI"],
  " 9, 218\[Dash]233.\nRatkowsky, David A. (1983). ",
  StyleBox["Nonlinear Regression Modeling, A Unified Practical Approach", 
    "TI"],
  ", Marcel Dekker, New York. "
}], "Text",
  CellTags->{"S6.5.1", "1.27"}]
},
FrontEndVersion->"4.0 for X",
ScreenRectangle->{{0, 800}, {0, 600}},
WindowSize->{520, 455},
WindowMargins->{{140, Automatic}, {Automatic, 72}},
StyleDefinitions -> "HelpBrowser.nb"
]


(***********************************************************************
Cached data follows.  If you edit this Notebook file directly, not using
Mathematica, you must remove the line containing CacheID at the top of 
the file.  The cache data will then be recreated when you save this file 
from within Mathematica.
***********************************************************************)

(*CellTagsOutline
CellTagsIndex->{
  "S6.5.1"->{
    Cell[1717, 49, 143, 6, 70, "Subsection",
      CellTags->{"S6.5.1", "1.1"}],
    Cell[1863, 57, 944, 21, 70, "Text",
      CellTags->{"S6.5.1", "1.2"}],
    Cell[2810, 80, 1865, 33, 70, "Text",
      CellTags->{"S6.5.1", "1.3"}],
    Cell[4678, 115, 162, 7, 70, "Caption",
      CellTags->"S6.5.1"],
    Cell[4843, 124, 1238, 30, 70, "Text",
      CellTags->{"S6.5.1", "1.4"}],
    Cell[6084, 156, 5138, 122, 70, "Text",
      CellTags->{"S6.5.1", "1.5"}],
    Cell[11225, 280, 174, 7, 70, "Caption",
      CellTags->"S6.5.1"],
    Cell[11402, 289, 7101, 234, 70, "Text",
      CellTags->{"S6.5.1", "1.6"}],
    Cell[18506, 525, 78, 1, 70, "MathCaption",
      CellTags->{"S6.5.1", "1.7"}],
    Cell[18587, 528, 91, 2, 70, "Input",
      CellTags->"S6.5.1"],
    Cell[18681, 532, 746, 22, 70, "MathCaption",
      CellTags->{"S6.5.1", "1.8"}],
    Cell[19430, 556, 168, 5, 70, "Input",
      CellTags->"S6.5.1"],
    Cell[19601, 563, 828, 23, 70, "MathCaption",
      CellTags->{"S6.5.1", "1.9"}],
    Cell[20454, 590, 174, 6, 70, "Input",
      CellTags->"S6.5.1"],
    Cell[20631, 598, 168, 4, 70, "Output",
      CellTags->"S6.5.1"],
    Cell[20814, 605, 181, 4, 70, "MathCaption",
      CellTags->{"S6.5.1", "1.10"}],
    Cell[21020, 613, 238, 7, 70, "Input",
      CellTags->"S6.5.1"],
    Cell[21261, 622, 203, 5, 70, "Output",
      CellTags->"S6.5.1"],
    Cell[21479, 630, 2840, 49, 70, "Text",
      CellTags->{"S6.5.1", "1.11"}],
    Cell[24322, 681, 159, 7, 70, "Caption",
      CellTags->"S6.5.1"],
    Cell[24484, 690, 3417, 113, 70, "Text",
      CellTags->{"S6.5.1", "1.12"}],
    Cell[27904, 805, 712, 16, 70, "Text",
      CellTags->{"S6.5.1", "1.13"}],
    Cell[28619, 823, 120, 5, 70, "Caption",
      CellTags->"S6.5.1"],
    Cell[28742, 830, 1642, 47, 70, "Text",
      CellTags->{"S6.5.1", "1.14"}],
    Cell[30387, 879, 517, 14, 70, "MathCaption",
      CellTags->{"S6.5.1", "1.15"}],
    Cell[30929, 897, 178, 6, 70, "Input",
      CellTags->"S6.5.1"],
    Cell[31110, 905, 3716, 77, 70, "Output",
      CellTags->"S6.5.1"],
    Cell[34841, 985, 4614, 91, 70, "Text",
      CellTags->{"S6.5.1", "1.16"}],
    Cell[39458, 1078, 227, 7, 70, "Caption",
      CellTags->"S6.5.1"],
    Cell[39688, 1087, 1107, 23, 70, "Text",
      CellTags->{"S6.5.1", "1.17"}],
    Cell[40798, 1112, 441, 11, 70, "Text",
      CellTags->{"S6.5.1", "1.18"}],
    Cell[41242, 1125, 107, 2, 70, "Caption",
      CellTags->"S6.5.1"],
    Cell[41352, 1129, 335, 8, 70, "Text",
      CellTags->{"S6.5.1", "1.19"}],
    Cell[41690, 1139, 1028, 22, 70, "Text",
      CellTags->{"S6.5.1", "1.20"}],
    Cell[42721, 1163, 76, 1, 70, "Caption",
      CellTags->"S6.5.1"],
    Cell[42800, 1166, 581, 13, 70, "Text",
      CellTags->{"S6.5.1", "1.21"}],
    Cell[43384, 1181, 891, 18, 70, "Text",
      CellTags->{"S6.5.1", "1.22"}],
    Cell[44278, 1201, 146, 4, 70, "Caption",
      CellTags->"S6.5.1"],
    Cell[44427, 1207, 3401, 71, 70, "Text",
      CellTags->{"S6.5.1", "1.23"}],
    Cell[47831, 1280, 447, 14, 70, "MathCaption",
      CellTags->{"S6.5.1", "1.24"}],
    Cell[48303, 1298, 237, 8, 70, "Input",
      CellTags->"S6.5.1"],
    Cell[48543, 1308, 1459, 31, 70, "Output",
      CellTags->"S6.5.1"],
    Cell[50017, 1342, 1926, 39, 70, "Text",
      CellTags->{"S6.5.1", "1.25"}],
    Cell[51946, 1383, 67, 1, 70, "Subsubsection",
      CellTags->{"S6.5.1", "1.26"}],
    Cell[52016, 1386, 586, 12, 70, "Text",
      CellTags->{"S6.5.1", "1.27"}]},
  "1.1"->{
    Cell[1717, 49, 143, 6, 70, "Subsection",
      CellTags->{"S6.5.1", "1.1"}]},
  "1.2"->{
    Cell[1863, 57, 944, 21, 70, "Text",
      CellTags->{"S6.5.1", "1.2"}]},
  "1.3"->{
    Cell[2810, 80, 1865, 33, 70, "Text",
      CellTags->{"S6.5.1", "1.3"}]},
  "1.4"->{
    Cell[4843, 124, 1238, 30, 70, "Text",
      CellTags->{"S6.5.1", "1.4"}]},
  "1.5"->{
    Cell[6084, 156, 5138, 122, 70, "Text",
      CellTags->{"S6.5.1", "1.5"}]},
  "1.6"->{
    Cell[11402, 289, 7101, 234, 70, "Text",
      CellTags->{"S6.5.1", "1.6"}]},
  "1.7"->{
    Cell[18506, 525, 78, 1, 70, "MathCaption",
      CellTags->{"S6.5.1", "1.7"}]},
  "1.8"->{
    Cell[18681, 532, 746, 22, 70, "MathCaption",
      CellTags->{"S6.5.1", "1.8"}]},
  "1.9"->{
    Cell[19601, 563, 828, 23, 70, "MathCaption",
      CellTags->{"S6.5.1", "1.9"}]},
  "1.10"->{
    Cell[20814, 605, 181, 4, 70, "MathCaption",
      CellTags->{"S6.5.1", "1.10"}]},
  "1.11"->{
    Cell[21479, 630, 2840, 49, 70, "Text",
      CellTags->{"S6.5.1", "1.11"}]},
  "1.12"->{
    Cell[24484, 690, 3417, 113, 70, "Text",
      CellTags->{"S6.5.1", "1.12"}]},
  "1.13"->{
    Cell[27904, 805, 712, 16, 70, "Text",
      CellTags->{"S6.5.1", "1.13"}]},
  "1.14"->{
    Cell[28742, 830, 1642, 47, 70, "Text",
      CellTags->{"S6.5.1", "1.14"}]},
  "1.15"->{
    Cell[30387, 879, 517, 14, 70, "MathCaption",
      CellTags->{"S6.5.1", "1.15"}]},
  "1.16"->{
    Cell[34841, 985, 4614, 91, 70, "Text",
      CellTags->{"S6.5.1", "1.16"}]},
  "1.17"->{
    Cell[39688, 1087, 1107, 23, 70, "Text",
      CellTags->{"S6.5.1", "1.17"}]},
  "1.18"->{
    Cell[40798, 1112, 441, 11, 70, "Text",
      CellTags->{"S6.5.1", "1.18"}]},
  "1.19"->{
    Cell[41352, 1129, 335, 8, 70, "Text",
      CellTags->{"S6.5.1", "1.19"}]},
  "1.20"->{
    Cell[41690, 1139, 1028, 22, 70, "Text",
      CellTags->{"S6.5.1", "1.20"}]},
  "1.21"->{
    Cell[42800, 1166, 581, 13, 70, "Text",
      CellTags->{"S6.5.1", "1.21"}]},
  "1.22"->{
    Cell[43384, 1181, 891, 18, 70, "Text",
      CellTags->{"S6.5.1", "1.22"}]},
  "1.23"->{
    Cell[44427, 1207, 3401, 71, 70, "Text",
      CellTags->{"S6.5.1", "1.23"}]},
  "1.24"->{
    Cell[47831, 1280, 447, 14, 70, "MathCaption",
      CellTags->{"S6.5.1", "1.24"}]},
  "1.25"->{
    Cell[50017, 1342, 1926, 39, 70, "Text",
      CellTags->{"S6.5.1", "1.25"}]},
  "1.26"->{
    Cell[51946, 1383, 67, 1, 70, "Subsubsection",
      CellTags->{"S6.5.1", "1.26"}]},
  "1.27"->{
    Cell[52016, 1386, 586, 12, 70, "Text",
      CellTags->{"S6.5.1", "1.27"}]}
  }
*)

(*CellTagsIndex
CellTagsIndex->{
  {"S6.5.1", 53230, 1417},
  {"1.1", 56703, 1508},
  {"1.2", 56796, 1511},
  {"1.3", 56884, 1514},
  {"1.4", 56973, 1517},
  {"1.5", 57063, 1520},
  {"1.6", 57154, 1523},
  {"1.7", 57246, 1526},
  {"1.8", 57341, 1529},
  {"1.9", 57438, 1532},
  {"1.10", 57536, 1535},
  {"1.11", 57634, 1538},
  {"1.12", 57727, 1541},
  {"1.13", 57821, 1544},
  {"1.14", 57913, 1547},
  {"1.15", 58006, 1550},
  {"1.16", 58105, 1553},
  {"1.17", 58198, 1556},
  {"1.18", 58292, 1559},
  {"1.19", 58385, 1562},
  {"1.20", 58477, 1565},
  {"1.21", 58571, 1568},
  {"1.22", 58664, 1571},
  {"1.23", 58757, 1574},
  {"1.24", 58851, 1577},
  {"1.25", 58951, 1580},
  {"1.26", 59045, 1583},
  {"1.27", 59145, 1586}
  }
*)

(*NotebookFileOutline
Notebook[{
Cell[1717, 49, 143, 6, 70, "Subsection",
  CellTags->{"S6.5.1", "1.1"}],
Cell[1863, 57, 944, 21, 70, "Text",
  CellTags->{"S6.5.1", "1.2"}],
Cell[2810, 80, 1865, 33, 70, "Text",
  CellTags->{"S6.5.1", "1.3"}],
Cell[4678, 115, 162, 7, 70, "Caption",
  CellTags->"S6.5.1"],
Cell[4843, 124, 1238, 30, 70, "Text",
  CellTags->{"S6.5.1", "1.4"}],
Cell[6084, 156, 5138, 122, 70, "Text",
  CellTags->{"S6.5.1", "1.5"}],
Cell[11225, 280, 174, 7, 70, "Caption",
  CellTags->"S6.5.1"],
Cell[11402, 289, 7101, 234, 70, "Text",
  CellTags->{"S6.5.1", "1.6"}],
Cell[18506, 525, 78, 1, 70, "MathCaption",
  CellTags->{"S6.5.1", "1.7"}],
Cell[18587, 528, 91, 2, 70, "Input",
  CellTags->"S6.5.1"],
Cell[18681, 532, 746, 22, 70, "MathCaption",
  CellTags->{"S6.5.1", "1.8"}],
Cell[19430, 556, 168, 5, 70, "Input",
  CellTags->"S6.5.1"],
Cell[19601, 563, 828, 23, 70, "MathCaption",
  CellTags->{"S6.5.1", "1.9"}],

Cell[CellGroupData[{
Cell[20454, 590, 174, 6, 70, "Input",
  CellTags->"S6.5.1"],
Cell[20631, 598, 168, 4, 70, "Output",
  CellTags->"S6.5.1"]
}, Open  ]],
Cell[20814, 605, 181, 4, 70, "MathCaption",
  CellTags->{"S6.5.1", "1.10"}],

Cell[CellGroupData[{
Cell[21020, 613, 238, 7, 70, "Input",
  CellTags->"S6.5.1"],
Cell[21261, 622, 203, 5, 70, "Output",
  CellTags->"S6.5.1"]
}, Open  ]],
Cell[21479, 630, 2840, 49, 70, "Text",
  CellTags->{"S6.5.1", "1.11"}],
Cell[24322, 681, 159, 7, 70, "Caption",
  CellTags->"S6.5.1"],
Cell[24484, 690, 3417, 113, 70, "Text",
  CellTags->{"S6.5.1", "1.12"}],
Cell[27904, 805, 712, 16, 70, "Text",
  CellTags->{"S6.5.1", "1.13"}],
Cell[28619, 823, 120, 5, 70, "Caption",
  CellTags->"S6.5.1"],
Cell[28742, 830, 1642, 47, 70, "Text",
  CellTags->{"S6.5.1", "1.14"}],
Cell[30387, 879, 517, 14, 70, "MathCaption",
  CellTags->{"S6.5.1", "1.15"}],

Cell[CellGroupData[{
Cell[30929, 897, 178, 6, 70, "Input",
  CellTags->"S6.5.1"],
Cell[31110, 905, 3716, 77, 70, "Output",
  CellTags->"S6.5.1"]
}, Open  ]],
Cell[34841, 985, 4614, 91, 70, "Text",
  CellTags->{"S6.5.1", "1.16"}],
Cell[39458, 1078, 227, 7, 70, "Caption",
  CellTags->"S6.5.1"],
Cell[39688, 1087, 1107, 23, 70, "Text",
  CellTags->{"S6.5.1", "1.17"}],
Cell[40798, 1112, 441, 11, 70, "Text",
  CellTags->{"S6.5.1", "1.18"}],
Cell[41242, 1125, 107, 2, 70, "Caption",
  CellTags->"S6.5.1"],
Cell[41352, 1129, 335, 8, 70, "Text",
  CellTags->{"S6.5.1", "1.19"}],
Cell[41690, 1139, 1028, 22, 70, "Text",
  CellTags->{"S6.5.1", "1.20"}],
Cell[42721, 1163, 76, 1, 70, "Caption",
  CellTags->"S6.5.1"],
Cell[42800, 1166, 581, 13, 70, "Text",
  CellTags->{"S6.5.1", "1.21"}],
Cell[43384, 1181, 891, 18, 70, "Text",
  CellTags->{"S6.5.1", "1.22"}],
Cell[44278, 1201, 146, 4, 70, "Caption",
  CellTags->"S6.5.1"],
Cell[44427, 1207, 3401, 71, 70, "Text",
  CellTags->{"S6.5.1", "1.23"}],
Cell[47831, 1280, 447, 14, 70, "MathCaption",
  CellTags->{"S6.5.1", "1.24"}],

Cell[CellGroupData[{
Cell[48303, 1298, 237, 8, 70, "Input",
  CellTags->"S6.5.1"],
Cell[48543, 1308, 1459, 31, 70, "Output",
  CellTags->"S6.5.1"]
}, Open  ]],
Cell[50017, 1342, 1926, 39, 70, "Text",
  CellTags->{"S6.5.1", "1.25"}],
Cell[51946, 1383, 67, 1, 70, "Subsubsection",
  CellTags->{"S6.5.1", "1.26"}],
Cell[52016, 1386, 586, 12, 70, "Text",
  CellTags->{"S6.5.1", "1.27"}]
}
]
*)




(***********************************************************************
End of Mathematica Notebook file.
***********************************************************************)

