(* Content-type: application/vnd.wolfram.mathematica *)

(*** Wolfram Notebook File ***)
(* http://www.wolfram.com/nb *)

(* CreatedBy='Mathematica 10.4' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     22107,        683]
NotebookOptionsPosition[     19379,        593]
NotebookOutlinePosition[     19926,        614]
CellTagsIndexPosition[     19883,        611]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell["Pure functions. LinearModelFit", "Title",
 CellChangeTimes->{{3.668251272826833*^9, 3.6682513045541053`*^9}, {
   3.6682516559473243`*^9, 3.66825166290457*^9}, {3.682842826258929*^9, 
   3.682842848254422*^9}, 3.7005577919607916`*^9}],

Cell[CellGroupData[{

Cell["\<\
Pure functions: when a function is needed in a very limited scope.\
\>", "Section",
 CellChangeTimes->{{3.5914514831933737`*^9, 3.5914514870374446`*^9}, {
  3.591458753632032*^9, 3.5914588062071047`*^9}, {3.6828428573982954`*^9, 
  3.682842885014832*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"Function", "[", 
   RowBox[{"x", ",", 
    RowBox[{"x", "^", "2"}]}], "]"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{
    RowBox[{
    "First", " ", "parameter", " ", "is", " ", "the", " ", "anonymous", " ", 
     "variable"}], ",", " ", 
    RowBox[{
    "which", " ", "is", " ", "transformed", " ", "according", " ", "the", " ",
      "second"}]}], " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.5914514915846996`*^9, 3.5914514996166196`*^9}, {
  3.5914516931015205`*^9, 3.591451721572647*^9}, {3.6198066124890056`*^9, 
  3.6198066166924314`*^9}, {3.682842894484006*^9, 3.682842942337006*^9}}],

Cell[BoxData[
 RowBox[{"%", "[", 
  RowBox[{"a", "+", "b"}], "]"}]], "Input",
 CellChangeTimes->{{3.5914515361821685`*^9, 3.591451544292222*^9}}],

Cell[CellGroupData[{

Cell["Now, with runes", "Subsection",
 CellChangeTimes->{{3.6682513891725593`*^9, 3.668251392417314*^9}, {
  3.668362436417982*^9, 3.6683624399810266`*^9}, {3.6828430397705708`*^9, 
  3.682843044189286*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{"#", "^", "2"}], "&"}], " ", 
  RowBox[{"(*", " ", 
   RowBox[{"#", ",", " ", 
    RowBox[{"equivalent", " ", "to", " ", "#1"}], ",", " ", 
    RowBox[{"is", " ", "the", " ", "parameter"}], ",", " ", 
    RowBox[{
     RowBox[{
      RowBox[{
       RowBox[{
       "therefore", " ", "it", " ", "is", " ", "already", " ", "in", " ", 
        "the", " ", "place", " ", "of", " ", "the", " ", "argument"}], " ", 
       "-", " ", 
       RowBox[{"the", " ", "ampersand"}]}], " ", "&"}], "  ", "tells", " ", 
     "the", " ", "program", " ", "that", " ", "it", " ", "is", " ", "a", " ", 
     "pure", " ", "function"}]}], " ", "*)"}]}]], "Input",
 CellChangeTimes->{{3.5914516440661745`*^9, 3.5914516479258747`*^9}, {
  3.5914517393553815`*^9, 3.5914518172212563`*^9}, {3.6828430488911467`*^9, 
  3.682843120064143*^9}}],

Cell[BoxData[
 RowBox[{"%", "[", 
  RowBox[{"a", "+", "b"}], "]"}]], "Input",
 CellChangeTimes->{{3.59145165533274*^9, 3.591451658176731*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Selecting values of a list", "Section",
 CellChangeTimes->{{3.591451876460576*^9, 3.591451888914741*^9}, {
  3.6682515333142056`*^9, 3.668251535131667*^9}, {3.6828431344951887`*^9, 
  3.682843140643964*^9}}],

Cell["a bigger test list for this function", "Text",
 CellChangeTimes->{{3.5914518958684464`*^9, 3.5914518975560875`*^9}, {
  3.682843144314788*^9, 3.6828431460301332`*^9}, {3.6828452103367157`*^9, 
  3.682845211690592*^9}, {3.682845301018221*^9, 3.682845315340974*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"testList", "=", 
   RowBox[{"RandomVariate", "[", 
    RowBox[{
     RowBox[{"NormalDistribution", "[", 
      RowBox[{"10", ",", "4"}], "]"}], ",", "1000"}], "]"}]}], ";"}]], "Input"],

Cell[BoxData[
 RowBox[{"Select", "[", 
  RowBox[{"testList", ",", 
   RowBox[{
    RowBox[{"#", "<", "0"}], "&"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.591451920604883*^9, 3.591451961155146*^9}, {
  3.668251477873813*^9, 3.6682514829462233`*^9}}],

Cell[BoxData[
 RowBox[{"Select", "[", 
  RowBox[{"testList", ",", " ", 
   RowBox[{
    RowBox[{"0.5", "<", "#", "<", "1"}], "&"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.5914519676088066`*^9, 3.5914519828288307`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["LinearModelFit", "Section",
 CellChangeTimes->{{3.7005530024233317`*^9, 3.700553005536691*^9}}],

Cell["An example with Angular Distribution", "Text",
 CellChangeTimes->{{3.6822243689328966`*^9, 3.682224377904173*^9}, {
  3.700552987945752*^9, 3.7005529916585455`*^9}}],

Cell[CellGroupData[{

Cell["Experimental data", "Subsection",
 CellChangeTimes->{{3.6822255722801604`*^9, 3.6822255747428355`*^9}}],

Cell[BoxData[
 RowBox[{"xws", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"0", ",", "1.17", ",", "0.12"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"0.5236", ",", "1.43", ",", "0.11"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"1.0472", ",", "1.16", ",", "0.09"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"1.5708", ",", "0.5", ",", "0.07"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"2.0944", ",", "0.95", ",", "0.09"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"2.6180", ",", "1.38", ",", "0.11"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"3.141592653589793", ",", "0.96", ",", "0.12"}], "}"}]}], 
   "}"}]}]], "Input",
 CellChangeTimes->{{3.6822248534514904`*^9, 3.6822249524825835`*^9}, {
  3.682227694716255*^9, 3.68222769536155*^9}}],

Cell[BoxData[
 RowBox[{"xw", "=", 
  RowBox[{"xws", "[", 
   RowBox[{"[", 
    RowBox[{"All", ",", 
     RowBox[{"{", 
      RowBox[{"1", ",", "2"}], "}"}]}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.6822246231512203`*^9, 3.6822246366555614`*^9}, {
  3.6822276982142553`*^9, 3.682227702232012*^9}}],

Cell[BoxData[
 RowBox[{"\[Sigma]", "=", 
  RowBox[{"xws", "[", 
   RowBox[{"[", 
    RowBox[{"All", ",", "3"}], "]"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.682224646475955*^9, 3.682224659236464*^9}, {
  3.6822277096240125`*^9, 3.682227710141057*^9}}],

Cell[BoxData[
 RowBox[{"Needs", "[", "\"\<ErrorBarPlots`\>\"", "]"}]], "Input"],

Cell[BoxData[
 RowBox[{"wDataPlot", "=", 
  RowBox[{"ErrorListPlot", "[", 
   RowBox[{"xws", ",", 
    RowBox[{"AxesLabel", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"\"\<\[Theta] (rd)\>\"", ",", "\"\<W\>\""}], "}"}]}]}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.682224767626813*^9, 3.6822248076915226`*^9}, {
  3.6822251233150535`*^9, 3.6822251315832396`*^9}, {3.6822277139147635`*^9, 
  3.682227714437364*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Fit: first unweighted, then using the variances", "Subsection",
 CellChangeTimes->{{3.6822255891074524`*^9, 3.6822256339312444`*^9}}],

Cell[BoxData[
 RowBox[{"aUnweightedFit", "=", 
  RowBox[{"LinearModelFit", "[", 
   RowBox[{"xw", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"LegendreP", "[", 
       RowBox[{"2", ",", 
        RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "]"}], ",", 
      RowBox[{"LegendreP", "[", 
       RowBox[{"4", ",", 
        RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "]"}]}], "}"}], ",", 
    RowBox[{"{", "\[Theta]", "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.682225046365136*^9, 3.682225096885393*^9}, {
  3.682227722690777*^9, 3.6822277232275505`*^9}}],

Cell[BoxData[
 RowBox[{"aUnweightedFit", "[", "\"\<ParameterTable\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.682225119459906*^9, 3.6822251494493504`*^9}}],

Cell[BoxData[
 RowBox[{"aUwFitPlot", "=", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{"aUnweightedFit", "[", "x", "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", "0", ",", "Pi"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.682225204779869*^9, 3.6822252379322968`*^9}}],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"wDataPlot", ",", "aUwFitPlot"}], "]"}]], "Input",
 CellChangeTimes->{{3.6822252508221006`*^9, 3.682225263528223*^9}}],

Cell[BoxData[
 RowBox[{"aFit", "=", 
  RowBox[{"LinearModelFit", "[", 
   RowBox[{"xw", ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"LegendreP", "[", 
       RowBox[{"2", ",", 
        RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "]"}], ",", 
      RowBox[{"LegendreP", "[", 
       RowBox[{"4", ",", 
        RowBox[{"Cos", "[", "\[Theta]", "]"}]}], "]"}]}], "}"}], ",", 
    RowBox[{"{", "\[Theta]", "}"}], ",", 
    RowBox[{"Weights", "\[Rule]", 
     RowBox[{"1", "/", 
      RowBox[{"\[Sigma]", "^", "2"}]}]}], ",", 
    RowBox[{"VarianceEstimatorFunction", "\[Rule]", 
     RowBox[{"(", 
      RowBox[{"1", "&"}], ")"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.682225283221429*^9, 3.6822253186977196`*^9}, {
  3.6822277353582864`*^9, 3.682227735999838*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Comparing weighted and unweighted parameter fits", "Subsection",
 CellChangeTimes->{{3.6822256594305396`*^9, 3.682225689903594*^9}}],

Cell[BoxData[
 RowBox[{"aFit", "[", "\"\<ParameterTable\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.682225327240012*^9, 3.6822253395098667`*^9}}],

Cell[BoxData[
 RowBox[{"aUnweightedFit", "[", "\"\<ParameterTable\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.6822253525528154`*^9, 3.6822253552861633`*^9}}],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"wDataPlot", ",", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"aUnweightedFit", "[", "x", "]"}], ",", 
       RowBox[{"aFit", "[", "x", "]"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", "Pi"}], "}"}], ",", 
     RowBox[{"PlotStyle", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{"Blue", ",", "Red"}], "}"}]}]}], "]"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.6822253884497333`*^9, 3.682225427803136*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["\<\
Mandatory tests - ALWAYS look at these graphics and numbers when using linear \
models\
\>", "Subsection",
 CellChangeTimes->{{3.6822261232024226`*^9, 3.682226197715024*^9}, {
  3.6822263336366067`*^9, 3.6822263398684473`*^9}}],

Cell[BoxData[
 RowBox[{"residuals", "=", 
  RowBox[{"Table", "[", " ", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"xw", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "1"}], "]"}], "]"}], ",", 
      RowBox[{
       RowBox[{"aFit", "[", "\"\<StandardizedResiduals\>\"", "]"}], "[", 
       RowBox[{"[", "i", "]"}], "]"}]}], " ", "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "7"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.682225698243644*^9, 3.6822257026090965`*^9}, {
  3.6822257382983055`*^9, 3.682225811499585*^9}, {3.6822277475280733`*^9, 
  3.6822277481408567`*^9}}],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"residuals", ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<\[Theta] (rd)\>\"", ",", "\"\<normalized residual\>\""}], 
     "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.6822258197328253`*^9, 3.6822258966964617`*^9}}],

Cell[BoxData[
 RowBox[{"residualsW", "=", 
  RowBox[{"Table", "[", " ", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"xw", "[", 
       RowBox[{"[", 
        RowBox[{"i", ",", "2"}], "]"}], "]"}], ",", 
      RowBox[{
       RowBox[{"aFit", "[", "\"\<StandardizedResiduals\>\"", "]"}], "[", 
       RowBox[{"[", "i", "]"}], "]"}]}], " ", "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "7"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.6822259360361433`*^9, 3.6822259397220383`*^9}, {
  3.682225989862965*^9, 3.682225996926787*^9}, {3.682228085269411*^9, 
  3.6822280854321613`*^9}}],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"residualsW", ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"\"\<w exp\>\"", ",", "\"\<normalized residual\>\""}], "}"}]}]}],
   "]"}]], "Input",
 CellChangeTimes->{{3.6822258197328253`*^9, 3.6822258966964617`*^9}, {
  3.6822259686820793`*^9, 3.6822259941095467`*^9}}],

Cell[BoxData[
 RowBox[{"chiSquare", "=", 
  RowBox[{"Total", "[", " ", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"aFit", "[", "\"\<FitResiduals\>\"", "]"}], "/", "\[Sigma]"}], 
     ")"}], "^", "2"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.6822260165366745`*^9, 3.68222603414252*^9}, {
  3.682226067513317*^9, 3.682226085412595*^9}}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["NonlinearModelFit", "Section",
 CellChangeTimes->{{3.68284329339003*^9, 3.68284329874432*^9}}],

Cell["an example with a simple sinus function", "Text",
 CellChangeTimes->{{3.7005531136839504`*^9, 3.700553157332909*^9}, {
  3.7005532074573474`*^9, 3.700553208160711*^9}}],

Cell[BoxData[
 RowBox[{"data", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{"0.00", ",", 
      RowBox[{"-", "0.8"}], ",", "0.5"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"0.25", ",", "9.5", ",", "0.8"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"0.5", ",", 
      RowBox[{"-", "7.3"}], ",", "0.7"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"0.75", ",", 
      RowBox[{"-", "2.3"}], ",", "0.5"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"1.00", ",", "11.0", ",", "0.9"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"1.25", ",", 
      RowBox[{"-", "5.5"}], ",", "0.6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"1.5", ",", 
      RowBox[{"-", "5.7"}], ",", "0.6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"1.75", ",", "10.4", ",", "0.9"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"2.00", ",", 
      RowBox[{"-", "2.3"}], ",", "0.5"}], "}"}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.6828443365478234`*^9, 3.682844409107017*^9}}],

Cell[BoxData[
 RowBox[{"fit", "=", 
  RowBox[{"NonlinearModelFit", "[", 
   RowBox[{
    RowBox[{"data", "[", 
     RowBox[{"[", 
      RowBox[{"All", ",", 
       RowBox[{"{", 
        RowBox[{"1", ",", "2"}], "}"}]}], "]"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"a", " ", 
      RowBox[{"Sin", "[", 
       RowBox[{"w", " ", "t"}], "]"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"{", 
       RowBox[{"a", ",", "10"}], "}"}], ",", 
      RowBox[{"{", 
       RowBox[{"w", ",", "8"}], "}"}]}], "}"}], ",", 
    RowBox[{"{", "t", "}"}], ",", 
    RowBox[{"Weights", "\[Rule]", 
     RowBox[{"1", "/", 
      RowBox[{
       RowBox[{"data", "[", 
        RowBox[{"[", 
         RowBox[{"All", ",", "3"}], "]"}], "]"}], "^", "2"}]}]}], ",", 
    RowBox[{"VarianceEstimatorFunction", "\[Rule]", 
     RowBox[{"(", 
      RowBox[{"1", "&"}], ")"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.682844414078909*^9, 3.6828445128216915`*^9}, {
  3.682844732526126*^9, 3.6828447828130245`*^9}}],

Cell[BoxData[
 RowBox[{"fit", "[", "\"\<ParameterTable\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.6828447134528666`*^9, 3.682844727024871*^9}}],

Cell[BoxData[
 RowBox[{"\[Chi]2", "=", 
  RowBox[{"Total", "[", 
   RowBox[{
    RowBox[{"(", 
     RowBox[{
      RowBox[{"fit", "[", "\"\<FitResiduals\>\"", "]"}], "/", 
      RowBox[{"data", "[", 
       RowBox[{"[", 
        RowBox[{"All", ",", "3"}], "]"}], "]"}]}], ")"}], "^", "2"}], 
   "]"}]}]], "Input",
 CellChangeTimes->{{3.6828449106866007`*^9, 3.6828449510604887`*^9}}],

Cell[BoxData[
 RowBox[{"p\[Chi]2", "=", 
  RowBox[{"100", 
   RowBox[{"(", 
    RowBox[{"1", "-", 
     RowBox[{"CDF", "[", 
      RowBox[{
       RowBox[{"ChiSquareDistribution", "[", 
        RowBox[{
         RowBox[{"Length", "[", "data", "]"}], "-", "2"}], "]"}], ",", 
       "\[Chi]2"}], "]"}]}], ")"}]}]}]], "Input",
 CellChangeTimes->{{3.682844954892821*^9, 3.68284499317365*^9}}],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"Transpose", "[", 
   RowBox[{"{", " ", 
    RowBox[{
     RowBox[{"data", "[", 
      RowBox[{"[", 
       RowBox[{"All", ",", "1"}], "]"}], "]"}], ",", 
     RowBox[{"fit", "[", "\"\<FitResiduals\>\"", "]"}]}], "}"}], "]"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.682845004131*^9, 3.6828450640035625`*^9}, {
  3.6828451027077007`*^9, 3.682845111078694*^9}}],

Cell[BoxData[
 RowBox[{"ListPlot", "[", 
  RowBox[{"Transpose", "[", 
   RowBox[{"{", " ", 
    RowBox[{
     RowBox[{"data", "[", 
      RowBox[{"[", 
       RowBox[{"All", ",", "2"}], "]"}], "]"}], ",", 
     RowBox[{"fit", "[", "\"\<FitResiduals\>\"", "]"}]}], "}"}], "]"}], 
  "]"}]], "Input",
 CellChangeTimes->{{3.6828451256169167`*^9, 3.68284512583237*^9}}],

Cell[BoxData[
 RowBox[{"Needs", "[", "\"\<ErrorBarPlots`\>\"", "]"}]], "Input",
 CellChangeTimes->{{3.6828451341677027`*^9, 3.6828451395569067`*^9}}],

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{
   RowBox[{"ErrorListPlot", "[", 
    RowBox[{"data", ",", 
     RowBox[{"PlotRange", "\[Rule]", 
      RowBox[{"{", 
       RowBox[{
        RowBox[{"-", "12"}], ",", "12"}], "}"}]}]}], "]"}], ",", 
   RowBox[{"Plot", "[", 
    RowBox[{
     RowBox[{"fit", "[", "x", "]"}], ",", 
     RowBox[{"{", 
      RowBox[{"x", ",", "0", ",", "2"}], "}"}]}], "]"}]}], "]"}]], "Input",
 CellChangeTimes->{{3.68284508676406*^9, 3.6828450936666965`*^9}, {
  3.682845129319728*^9, 3.6828451961435003`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell["Exercise", "Section",
 CellChangeTimes->{{3.6828445553919277`*^9, 3.6828445572606373`*^9}}],

Cell[TextData[{
 "Ajuste os par\[AHat]metros do polin\[OHat]mio ",
 Cell[BoxData[
  FormBox[
   RowBox[{"y", "=", 
    RowBox[{"a", "+", 
     RowBox[{"b", " ", "x"}], " ", "+", " ", 
     RowBox[{"c", " ", 
      SuperscriptBox["x", "3"], " "}]}]}], TraditionalForm]],
  FormatType->"TraditionalForm"],
 "ao conjunto de dados z= {x,y} abaixo. Os desvios-padr~eos s\[ATilde]o todos \
iguais a 0.5. Calcule ",
 Cell[BoxData[
  FormBox[
   RowBox[{
    SuperscriptBox["\[Chi]", "2"], " "}], TraditionalForm]],
  FormatType->"TraditionalForm"],
 "e fa\[CCedilla]a os gr\[AAcute]ficos necess\[AAcute]rios \[AGrave] avalia\
\[CCedilla]\[ATilde]o da qualidade do ajuste."
}], "Text",
 CellChangeTimes->{{3.682844566459696*^9, 3.6828446106409965`*^9}, {
  3.7005532917722836`*^9, 3.7005533625321827`*^9}, {3.7005567960877156`*^9, 
  3.700556809192773*^9}, {3.7005571832984114`*^9, 3.700557255166252*^9}}],

Cell[BoxData[
 RowBox[{"z", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "3.0"}], ",", "8.7"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "2.5"}], ",", "5.6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "2.0"}], ",", "4.6"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1.5"}], ",", "4.4"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "1.0"}], ",", "3.3"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "0.5"}], ",", "2.8"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"0", ",", "1.7"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"0.5", ",", "2.8"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"1.0", ",", "1.8"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"1.5", ",", "0.9"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"2.0", ",", "0.8"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"2.5", ",", 
      RowBox[{"-", "0.4"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"3.0", ",", 
      RowBox[{"-", "2.7"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"3.5", ",", 
      RowBox[{"-", "5.3"}]}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"4.0", ",", 
      RowBox[{"-", "6.4"}]}], "}"}]}], "}"}]}]], "Input"]
}, Open  ]]
}, Open  ]]
},
WindowSize->{1040, 709},
WindowMargins->{{420, Automatic}, {90, Automatic}},
PrintingCopies->1,
PrintingPageRange->{32000, 32000},
PrintingOptions->{"Magnification"->1.,
"PaperOrientation"->"Portrait",
"PaperSize"->{595.1999999999999, 841.68}},
Magnification:>1.3 Inherited,
FrontEndVersion->"11.0 for Microsoft Windows (64-bit) (September 21, 2016)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 240, 3, 131, "Title"],
Cell[CellGroupData[{
Cell[845, 29, 265, 5, 91, "Section"],
Cell[1113, 36, 629, 15, 65, "Input"],
Cell[1745, 53, 145, 3, 39, "Input"],
Cell[CellGroupData[{
Cell[1915, 60, 207, 3, 61, "Subsection"],
Cell[2125, 65, 863, 20, 90, "Input"],
Cell[2991, 87, 142, 3, 39, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[3182, 96, 213, 3, 91, "Section"],
Cell[3398, 101, 270, 3, 39, "Text"],
Cell[3671, 106, 219, 6, 39, "Input"],
Cell[3893, 114, 250, 6, 39, "Input"],
Cell[4146, 122, 220, 5, 39, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[4403, 132, 101, 1, 91, "Section"],
Cell[4507, 135, 171, 2, 39, "Text"],
Cell[CellGroupData[{
Cell[4703, 141, 109, 1, 61, "Subsection"],
Cell[4815, 144, 779, 20, 65, "Input"],
Cell[5597, 166, 305, 8, 39, "Input"],
Cell[5905, 176, 254, 6, 39, "Input"],
Cell[6162, 184, 79, 1, 39, "Input"],
Cell[6244, 187, 419, 10, 39, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[6700, 202, 139, 1, 61, "Subsection"],
Cell[6842, 205, 560, 14, 39, "Input"],
Cell[7405, 221, 156, 2, 39, "Input"],
Cell[7564, 225, 286, 7, 39, "Input"],
Cell[7853, 234, 165, 3, 39, "Input"],
Cell[8021, 239, 772, 20, 65, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[8830, 264, 138, 1, 61, "Subsection"],
Cell[8971, 267, 146, 2, 39, "Input"],
Cell[9120, 271, 158, 2, 39, "Input"],
Cell[9281, 275, 522, 14, 39, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9840, 294, 237, 5, 98, "Subsection"],
Cell[10080, 301, 611, 16, 39, "Input"],
Cell[10694, 319, 303, 7, 39, "Input"],
Cell[11000, 328, 610, 16, 39, "Input"],
Cell[11613, 346, 346, 8, 39, "Input"],
Cell[11962, 356, 352, 9, 39, "Input"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[12363, 371, 100, 1, 91, "Section"],
Cell[12466, 374, 174, 2, 39, "Text"],
Cell[12643, 378, 981, 28, 65, "Input"],
Cell[13627, 408, 1013, 30, 65, "Input"],
Cell[14643, 440, 145, 2, 39, "Input"],
Cell[14791, 444, 383, 11, 39, "Input"],
Cell[15177, 457, 389, 11, 39, "Input"],
Cell[15569, 470, 413, 11, 39, "Input"],
Cell[15985, 483, 364, 10, 39, "Input"],
Cell[16352, 495, 149, 2, 39, "Input"],
Cell[16504, 499, 541, 15, 39, "Input"]
}, Open  ]],
Cell[CellGroupData[{
Cell[17082, 519, 97, 1, 91, "Section"],
Cell[17182, 522, 897, 22, 68, "Text"],
Cell[18082, 546, 1269, 43, 65, "Input"]
}, Open  ]]
}, Open  ]]
}
]
*)