Exercise28.2.nb

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

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

(* CreatedBy='Mathematica 10.1' *)

(*CacheID: 234*)
(* Internal cache information:
NotebookFileLineBreakTest
NotebookFileLineBreakTest
NotebookDataPosition[       158,          7]
NotebookDataLength[     30527,        971]
NotebookOptionsPosition[     27335,        855]
NotebookOutlinePosition[     27744,        873]
CellTagsIndexPosition[     27701,        870]
WindowFrame->Normal*)

(* Beginning of Notebook Content *)
Notebook[{

Cell[CellGroupData[{
Cell[TextData[{
 "Evidence of ",
 Cell[BoxData[
  FormBox[
   RowBox[{
    SubscriptBox["H", "1"], ",", " ", 
    RowBox[{"P", "(", 
     RowBox[{"D", "|", 
      SubscriptBox["H", "1"]}], ")"}]}], TraditionalForm]],
  FormatType->"TraditionalForm"]
}], "Subchapter",
 CellChangeTimes->{{3.7009312250855007`*^9, 3.700931237388666*^9}, {
  3.7009394389791937`*^9, 3.7009394482944527`*^9}, {3.7009394791270742`*^9, 
  3.700939501365468*^9}}],

Cell[BoxData[
 RowBox[{"Clear", "[", 
  RowBox[{"x", ",", "y", ",", "t", ",", "w1", ",", "w0"}], "]"}]], "Input",
 CellChangeTimes->{{3.7009520866010027`*^9, 3.7009521001164207`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"data", "=", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "8"}], ",", "8"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"-", "2"}], ",", "10"}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{"6", ",", "11"}], "}"}]}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{3.700951333023745*^9}],

Cell[BoxData[
 RowBox[{
  RowBox[{"y", "=", 
   RowBox[{"w0", "+", 
    RowBox[{"w1", " ", "x"}]}]}], ";"}]], "Input",
 CellChangeTimes->{{3.7009313575196657`*^9, 3.7009313895665007`*^9}, {
  3.700937636531591*^9, 3.70093763663034*^9}, {3.700952079169611*^9, 
  3.700952080099824*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"w1", "=", "0"}], ";", 
  RowBox[{
   SubsuperscriptBox["\[Integral]", 
    RowBox[{"-", "\[Infinity]"}], "\[Infinity]"], 
   RowBox[{
    SqrtBox[
     FractionBox["1", 
      RowBox[{"2", "\[Pi]"}]]], 
    RowBox[{"Exp", "[", 
     FractionBox[
      RowBox[{"-", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"y", "-", "t"}], ")"}], "2"]}], "2"], "]"}], 
    SqrtBox[
     FractionBox["1", 
      RowBox[{"2", "\[Pi]"}]]], 
    RowBox[{"Exp", "[", 
     FractionBox[
      RowBox[{"-", 
       SuperscriptBox[
        RowBox[{"(", "w0", ")"}], "2"]}], "2"], "]"}], 
    RowBox[{"\[DifferentialD]", "w0"}]}]}]}]], "Input",
 CellChangeTimes->{{3.700931239508768*^9, 3.700931350656206*^9}, {
  3.700931395231771*^9, 3.7009314550521193`*^9}, {3.700937544994338*^9, 
  3.700937616115135*^9}, {3.7009390737539873`*^9, 3.700939091343781*^9}, {
  3.700951958989813*^9, 3.70095206640495*^9}, {3.70095211882808*^9, 
  3.7009521193230257`*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   FractionBox[
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{"-", 
      FractionBox[
       SuperscriptBox["t", "2"], "4"]}]], 
    RowBox[{"2", " ", 
     SqrtBox["\[Pi]"]}]], "/.", 
   RowBox[{"t", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"8", ",", "10", ",", "11"}], "}"}]}]}], "//", "N"}]], "Input",
 CellChangeTimes->{{3.7009521305108023`*^9, 3.700952153408086*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "3.1745586679666396`*^-8", ",", "3.917716632754334`*^-12", ",", 
   "2.0558290113157305`*^-14"}], "}"}]], "Output",
 CellChangeTimes->{{3.700952150603806*^9, 3.700952153849834*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"L1", "=", 
   RowBox[{"{", 
    RowBox[{
    "3.1745586679666396`*^-8", ",", "3.917716632754334`*^-12", ",", 
     "2.0558290113157305`*^-14"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.700940747599984*^9, 3.700940752873045*^9}, {
   3.7009408354899673`*^9, 3.700940835619882*^9}, 3.700952160984796*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"evidence1", "=", 
  RowBox[{"Product", "[", 
   RowBox[{
    RowBox[{"L1", "[", 
     RowBox[{"[", "i", "]"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.700940737181302*^9, 3.700940769140851*^9}, {
  3.700940838500054*^9, 3.7009408386279488`*^9}, {3.7009408741652393`*^9, 
  3.7009408794024754`*^9}, {3.700952022326285*^9, 3.70095202263698*^9}}],

Cell[BoxData["2.5568389192915458`*^-33"], "Output",
 CellChangeTimes->{3.700940769412107*^9, 3.7009408391717854`*^9, 
  3.700940879740674*^9, 3.700952162683248*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell[TextData[{
 "Evidence of ",
 Cell[BoxData[
  FormBox[
   RowBox[{
    SubscriptBox["H", "2"], ",", " ", 
    RowBox[{"P", "(", 
     RowBox[{"D", "|", 
      SubscriptBox["H", "2"]}], ")"}]}], TraditionalForm]],
  FormatType->"TraditionalForm"]
}], "Subchapter",
 CellChangeTimes->{{3.7009312250855007`*^9, 3.700931237388666*^9}, {
  3.7009394389791937`*^9, 3.7009394482944527`*^9}, {3.7009394791270742`*^9, 
  3.700939523379356*^9}}],

Cell[BoxData[
 RowBox[{"Clear", "[", 
  RowBox[{"w0", ",", "w1", ",", "x", ",", "y", ",", "t"}], "]"}]], "Input",
 CellChangeTimes->{{3.700931477121502*^9, 3.700931479532435*^9}, {
  3.7009377501790037`*^9, 3.700937750465947*^9}, {3.700952233866604*^9, 
  3.7009522437475758`*^9}}],

Cell[BoxData[
 RowBox[{
  RowBox[{"y", "=", 
   RowBox[{"w0", "+", 
    RowBox[{"w1", " ", "x"}]}]}], ";"}]], "Input"],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  SqrtBox[
   FractionBox["1", 
    RowBox[{"2", "\[Pi]"}]]], 
  RowBox[{"Exp", "[", 
   FractionBox[
    RowBox[{"-", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"t", "-", "y"}], ")"}], "2"]}], "2"], "]"}], 
  SqrtBox[
   FractionBox["1", 
    RowBox[{"2", "\[Pi]"}]]], 
  RowBox[{"Exp", "[", 
   FractionBox[
    RowBox[{"-", 
     SuperscriptBox["w1", "2"]}], "2"], "]"}], 
  SqrtBox[
   FractionBox["1", 
    RowBox[{"2", "\[Pi]"}]]], 
  RowBox[{"Exp", "[", 
   FractionBox[
    RowBox[{"-", 
     SuperscriptBox[
      RowBox[{"(", "w0", ")"}], "2"]}], "2"], "]"}]}]], "Input",
 CellChangeTimes->{{3.700931492158019*^9, 3.70093149601506*^9}, {
   3.700931585943928*^9, 3.7009316258443127`*^9}, {3.7009376457347393`*^9, 
   3.700937647911757*^9}, 3.700937739134253*^9, {3.700952200261134*^9, 
   3.7009522265716677`*^9}, {3.7009523061481447`*^9, 3.700952306893054*^9}}],

Cell[BoxData[
 FractionBox[
  SuperscriptBox["\[ExponentialE]", 
   RowBox[{
    RowBox[{"-", 
     FractionBox[
      SuperscriptBox["w0", "2"], "2"]}], "-", 
    FractionBox[
     SuperscriptBox["w1", "2"], "2"], "-", 
    RowBox[{
     FractionBox["1", "2"], " ", 
     SuperscriptBox[
      RowBox[{"(", 
       RowBox[{"t", "-", "w0", "-", 
        RowBox[{"w1", " ", "x"}]}], ")"}], "2"]}]}]], 
  RowBox[{"2", " ", 
   SqrtBox["2"], " ", 
   SuperscriptBox["\[Pi]", 
    RowBox[{"3", "/", "2"}]]}]]], "Output",
 CellChangeTimes->{{3.70093765419534*^9, 3.7009376645777063`*^9}, {
   3.700937739847349*^9, 3.7009377525204268`*^9}, 3.700952249633843*^9, 
   3.700952319052369*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  SubsuperscriptBox["\[Integral]", 
   RowBox[{"-", "\[Infinity]"}], "\[Infinity]"], 
  RowBox[{
   SubsuperscriptBox["\[Integral]", 
    RowBox[{"-", "\[Infinity]"}], "\[Infinity]"], 
   RowBox[{
    SqrtBox[
     FractionBox["1", 
      RowBox[{"2", "\[Pi]"}]]], 
    RowBox[{"Exp", "[", 
     FractionBox[
      RowBox[{"-", 
       SuperscriptBox[
        RowBox[{"(", 
         RowBox[{"t", "-", "y"}], ")"}], "2"]}], "2"], "]"}], 
    SqrtBox[
     FractionBox["1", 
      RowBox[{"2", "\[Pi]"}]]], 
    RowBox[{"Exp", "[", 
     FractionBox[
      RowBox[{"-", 
       SuperscriptBox["w1", "2"]}], "2"], "]"}], 
    SqrtBox[
     FractionBox["1", 
      RowBox[{"2", "\[Pi]"}]]], 
    RowBox[{"Exp", "[", 
     FractionBox[
      RowBox[{"-", 
       SuperscriptBox[
        RowBox[{"(", "w0", ")"}], "2"]}], "2"], "]"}], 
    RowBox[{"\[DifferentialD]", "w1"}], 
    RowBox[{"\[DifferentialD]", "w0"}]}]}]}]], "Input",
 CellChangeTimes->{{3.700931653114694*^9, 3.700931670782487*^9}, {
  3.700937692196286*^9, 3.700937724082734*^9}, {3.700952259443089*^9, 
  3.700952270753003*^9}, {3.700952323214621*^9, 3.7009523241509647`*^9}}],

Cell[BoxData[
 RowBox[{"ConditionalExpression", "[", 
  RowBox[{
   FractionBox[
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{"-", 
      FractionBox[
       SuperscriptBox["t", "2"], 
       RowBox[{"4", "+", 
        RowBox[{"2", " ", 
         SuperscriptBox["x", "2"]}]}]]}]], 
    RowBox[{
     SqrtBox[
      RowBox[{"2", " ", "\[Pi]"}]], " ", 
     SqrtBox[
      RowBox[{"1", "+", 
       SuperscriptBox["x", "2"]}]], " ", 
     SqrtBox[
      RowBox[{"1", "+", 
       FractionBox["1", 
        RowBox[{"1", "+", 
         SuperscriptBox["x", "2"]}]]}]]}]], ",", 
   RowBox[{
    RowBox[{"Re", "[", 
     FractionBox["1", 
      RowBox[{"1", "+", 
       SuperscriptBox["x", "2"]}]], "]"}], "\[GreaterEqual]", 
    RowBox[{"-", "1"}]}]}], "]"}]], "Output",
 CellChangeTimes->{3.700952271694881*^9, 3.700952315896832*^9, 
  3.700952366859705*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   FractionBox[
    SuperscriptBox["\[ExponentialE]", 
     RowBox[{"-", 
      FractionBox[
       SuperscriptBox["t", "2"], 
       RowBox[{"4", "+", 
        RowBox[{"2", " ", 
         SuperscriptBox["x", "2"]}]}]]}]], 
    RowBox[{
     SqrtBox[
      RowBox[{"2", " ", "\[Pi]"}]], " ", 
     SqrtBox[
      RowBox[{"1", "+", 
       SuperscriptBox["x", "2"]}]], " ", 
     SqrtBox[
      RowBox[{"1", "+", 
       FractionBox["1", 
        RowBox[{"1", "+", 
         SuperscriptBox["x", "2"]}]]}]]}]], "/.", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x", "\[Rule]", 
        RowBox[{"-", "8"}]}], ",", 
       RowBox[{"t", "\[Rule]", "8"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x", "\[Rule]", 
        RowBox[{"-", "2"}]}], ",", 
       RowBox[{"t", "\[Rule]", "10"}]}], "}"}], ",", 
     RowBox[{"{", 
      RowBox[{
       RowBox[{"x", "\[Rule]", "6"}], ",", 
       RowBox[{"t", "\[Rule]", "11"}]}], "}"}]}], "}"}]}], "//", 
  "N"}]], "Input",
 CellChangeTimes->{{3.7009523818768*^9, 3.700952398387801*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.03023925424859377`", ",", "0.00003914837665445475`", ",", 
   "0.013169695605323387`"}], "}"}]], "Output",
 CellChangeTimes->{{3.70095239470648*^9, 3.7009523988158627`*^9}}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"L2", "=", 
   RowBox[{"{", 
    RowBox[{
    "0.03023925424859377`", ",", "0.00003914837665445475`", ",", 
     "0.013169695605323387`"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.70094079196052*^9, 3.7009407936144753`*^9}, {
   3.7009408432548103`*^9, 3.7009408433546753`*^9}, 3.700952416593451*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"evidence2", "=", 
  RowBox[{"Product", "[", 
   RowBox[{
    RowBox[{"L2", "[", 
     RowBox[{"[", "i", "]"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.700940737181302*^9, 3.7009407878375273`*^9}, {
  3.7009408450998573`*^9, 3.700940845179205*^9}, {3.700940883360998*^9, 
  3.700940886645372*^9}}],

Cell[BoxData["1.5590518959711003`*^-8"], "Output",
 CellChangeTimes->{{3.700940769412107*^9, 3.7009407962889557`*^9}, {
   3.7009408455983877`*^9, 3.700940847828725*^9}, 3.700940887039351*^9, 
   3.7009524182458963`*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Linear Regression and evidence from the best fit", "Subchapter",
 CellChangeTimes->{{3.7009312250855007`*^9, 3.700931237388666*^9}, {
  3.7009394389791937`*^9, 3.7009394482944527`*^9}, {3.7009394791270742`*^9, 
  3.700939501365468*^9}, {3.700939535545888*^9, 3.700939555958934*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"data", "=", 
  RowBox[{"{", 
   RowBox[{
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "8"}], ",", "8"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"-", "2"}], ",", "10"}], "}"}], ",", 
    RowBox[{"{", 
     RowBox[{"6", ",", "11"}], "}"}]}], "}"}]}]], "Input",
 CellChangeTimes->{{3.700937906708914*^9, 3.70093792047473*^9}, {
  3.700938372512368*^9, 3.700938375251837*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "8"}], ",", "8"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "2"}], ",", "10"}], "}"}], ",", 
   RowBox[{"{", 
    RowBox[{"6", ",", "11"}], "}"}]}], "}"}]], "Output",
 CellChangeTimes->{3.700938375943356*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"P1", "=", 
  RowBox[{"ListPlot", "[", 
   RowBox[{"data", ",", " ", 
    RowBox[{"PlotRange", "\[Rule]", 
     RowBox[{"{", 
      RowBox[{"0", ",", "15"}], "}"}]}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.7009386205581923`*^9, 3.7009386813097963`*^9}, {
  3.700944772211248*^9, 3.700944792481885*^9}, {3.700944826003543*^9, 
  3.7009448383668528`*^9}}],

Cell[BoxData[
 GraphicsBox[{{}, {{}, 
    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.012833333333333334`],
      AbsoluteThickness[1.6], 
     PointBox[{{-8., 8.}, {-2., 10.}, {6., 11.}}]}, {}}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{},
  PlotRange->{{-8., 6.}, {0, 15}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {0, 0}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.700938624671417*^9, 3.700938681614414*^9}, 
   3.700944183983465*^9, {3.7009447825801067`*^9, 3.700944792928329*^9}, {
   3.7009448341670094`*^9, 3.700944838873356*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"P1", ",", 
   RowBox[{"AxesLabel", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{
      RowBox[{"HoldForm", "[", "x", "]"}], ",", 
      RowBox[{"HoldForm", "[", "y", "]"}]}], "}"}]}], ",", 
   RowBox[{"PlotRange", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"0", ",", "15"}], "}"}]}], ",", 
   RowBox[{"PlotLabel", "\[Rule]", "None"}], ",", 
   RowBox[{"LabelStyle", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"GrayLevel", "[", "0", "]"}], "}"}]}]}], "]"}]], "Input",
 CellChangeTimes->{{3.700944856750125*^9, 3.700944858058929*^9}},
 NumberMarks->False],

Cell[BoxData[
 GraphicsBox[{{}, {{}, 
    {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.012833333333333334`],
      AbsoluteThickness[1.6], 
     PointBox[{{-8., 8.}, {-2., 10.}, {6., 11.}}]}, {}}, {}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{
    FormBox[
     TagBox["x", HoldForm], TraditionalForm], 
    FormBox[
     TagBox["y", HoldForm], TraditionalForm]},
  AxesOrigin->{0, 0},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  LabelStyle->{
    GrayLevel[0]},
  Method->{},
  PlotLabel->None,
  PlotRange->{0, 15},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {0, 0}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{3.7009441991713963`*^9, 3.700944860420968*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"model", "=", 
  RowBox[{"LinearModelFit", "[", 
   RowBox[{"data", ",", "x", ",", "x"}], "]"}]}]], "Input",
 CellChangeTimes->{{3.70093849789102*^9, 3.7009385744216423`*^9}}],

Cell[BoxData[
 TagBox[
  RowBox[{"FittedModel", "[", 
   TagBox[
    PanelBox[
     TagBox[
      RowBox[{"9.945945945945946`", "\[VeryThinSpace]", "+", 
       RowBox[{"0.20945945945945946`", " ", "x"}]}],
      Short[#, 2]& ],
     FrameMargins->5],
    Editable -> False], "]"}],
  InterpretTemplate[
  FittedModel[{
    "Linear", {9.945945945945946, 
     0.20945945945945946`}, {{$CellContext`x}, {1, $CellContext`x}}, {0, 
     0}}, {{1., 1., 1.}}, {{-8, 8}, {-2, 10}, {6, 11}}, {{1., -8.}, {
    1., -2.}, {1., 6.}}, 
    Function[Null, 
     Internal`LocalizedBlock[{$CellContext`x}, #], {HoldAll}]]& ],
  Editable->False,
  SelectWithContents->True,
  Selectable->True]], "Output",
 CellChangeTimes->{
  3.700938377192505*^9, {3.700938501126392*^9, 3.7009385748453417`*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"model", "[", "\"\<BestFit\>\"", "]"}]], "Input"],

Cell[BoxData[
 RowBox[{"9.945945945945946`", "\[VeryThinSpace]", "+", 
  RowBox[{"0.20945945945945946`", " ", "x"}]}]], "Output",
 CellChangeTimes->{3.700938390730682*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"P2", "=", 
  RowBox[{"Plot", "[", 
   RowBox[{
    RowBox[{"9.945945945945946`", "\[VeryThinSpace]", "+", 
     RowBox[{"0.20945945945945946`", " ", "x"}]}], ",", 
    RowBox[{"{", 
     RowBox[{"x", ",", 
      RowBox[{"-", "9"}], ",", "12"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.700938683477749*^9, 3.700938781685828*^9}, 
   3.700938837912992*^9}],

Cell[BoxData[
 GraphicsBox[{{}, {}, 
   {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
    1.], LineBox[CompressedData["
1:eJwllHs4lXcAx1VHsS7mRJLb+77n8nsnkkhIfEOLlLCIbaKe7CJlEUlLpFxy
XzEUSRcprXJWKtZ5k9GFjnUqmZykMSexFKHEzp798X2+z+efz/N8//nSm8J9
QiarqamtUeW/ji4oHpiYoLmiHfZHGBMKDe0BNyVjNHfgxthWuYrnMro5347S
3FYNKfZTFK6fTVsoe0Nzy466dHfRFMarY8KLu2hOcWu1ZbmQQrLCp8/hHs0Z
zw5qsDClkC+YptyVR3MllxLfOiymcO18eMdbM5o72d/osd6Dwp/v6pyJG8Vt
2xLR6bKbgrTc21Ojy5gT5vrtmHKGwouTFjUCZyNulhW3tlNO4ZbmHO0buw24
GrmD9q4pNMLy9UPPyfU5X3qkqdOCRnO26KyexlxOsr3+B14wDc/bK7O6XOdw
rjbzAldk0Tj8WXLM1BId7mZZhWtANY28wftup1v43CRDH76ol0aw0sqowFyb
W7DMpu2WHoMU+w3mbXpaXOaDY30JLgwaokL6EzKmc5Whht3DEQx+5pVV3XfS
4GovtLTiGAPFN6Oe+ud4XN7wKclyGQP/lTELCZnEFUlGH28cZ1D/bpP1FekH
aY5h+tA6UwHuqg+f0D8xKHVLsXw/ESBAx86JlT9N9Ep7S3Tf70kS4NPNZ9T+
sGuTetfM6JBcEWBv8ImHiX7ZTq88UrU0uwXItwq12x7U5hSR4mpvyRei6m5C
5+8jvU41DXwPbWchLuZVJhaVDTq9kZisr9smRLHHvQRvyw9OX+tUVuQVCiEz
ez6wXzgJqTvX91B3hXAMLp1XHclDafJSk7h3QthwPlFvNTWwt3CRwyGxCB+O
Pz7wXdZ0ZNMXSht8RbBYSGUueTILGb6hPFmCCNFRZpt51tpYk267OvOSCDPK
538saOdjNHxFFqUQIbVIZ6ApTgetDS0GDzTFWMUOfdznMgd7mD2ZS+3EcEm9
OrVLqYeZphvHdoaIER8fyd6W6CPY2t8rLFeMoTXan/sHGaCkN+OYVa0YXk/z
byRNNgJfz912Xr8Yr9XmetncMYb7Cr2qq0YEf32s2XbHlkIYJbQ0XUUQUa9c
ftWYhs+4S84XsQSR4WnO0WoMtnop7Y+WE2QnGt/8RM4gcm3Ts6pHBLJN5o1B
VQI4eNeFpqmzMD1v7puWLsSbidHXlotYjLSUTk8KFOEpc5kbDmTR03TIcznE
SHFLjI7KYDG4IFpLPpuAMcgm9ddYFDeP1ybqEahnv5+Wfp2FW0RytLUBgZIX
0uNdzeLo5fz2XIbgYr9deXsNC1eH6nP+lgSOtS9Mh6QsDruruSk8Cb76fskC
YT2LxZsPxitTCXKrFNb7HrDoUJ9tXZhBsMvMXddNzuJg2ZG/V+UQBB6XDM18
yELxssKzIp9AmJZypeARi6TtMsPwMoLKDYtsLz1h8ThO59pgnWq3etLSDgWL
eKoo7NRtAknsa4PTz1T7a0WUXyNB3j9fjm3pYBGnviS5Sk6wodXit+HnLMTp
Ab6xnQSvKlodtbpYNJt3aszvJmimXU0eqThWFlrTplT58n6ZKOxmcZ//o2DZ
AEFsfCIn7mER8yuvpW9Q5RvqK3mlYsYv42DxCIFzqH9CpZJF47Cu49oxAtGz
2o0xL1n8/18EmuvMnR17WfwLDjcoGQ==
     "]]}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 8.},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{
   "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> 
    AbsolutePointSize[6], "ScalingFunctions" -> None},
  PlotRange->{{-9, 12}, {8.06081090057915, 12.45945936969112}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.7009386929868393`*^9, 3.700938782147313*^9}, 
   3.7009388387847643`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"Show", "[", 
  RowBox[{"P1", ",", "P2"}], "]"}]], "Input",
 CellChangeTimes->{{3.700938696021535*^9, 3.700938700189583*^9}}],

Cell[BoxData[
 GraphicsBox[{{{}, {{}, 
     {RGBColor[0.368417, 0.506779, 0.709798], PointSize[
      0.012833333333333334`], AbsoluteThickness[1.6], 
      PointBox[{{-8., 8.}, {-2., 10.}, {6., 11.}}]}, {}}, {}}, {{}, {}, 
    {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[
     1.], LineBox[CompressedData["
1:eJwllHs4lXcAx1VHsS7mRJLb+77n8nsnkkhIfEOLlLCIbaKe7CJlEUlLpFxy
XzEUSRcprXJWKtZ5k9GFjnUqmZykMSexFKHEzp798X2+z+efz/N8//nSm8J9
QiarqamtUeW/ji4oHpiYoLmiHfZHGBMKDe0BNyVjNHfgxthWuYrnMro5347S
3FYNKfZTFK6fTVsoe0Nzy466dHfRFMarY8KLu2hOcWu1ZbmQQrLCp8/hHs0Z
zw5qsDClkC+YptyVR3MllxLfOiymcO18eMdbM5o72d/osd6Dwp/v6pyJG8Vt
2xLR6bKbgrTc21Ojy5gT5vrtmHKGwouTFjUCZyNulhW3tlNO4ZbmHO0buw24
GrmD9q4pNMLy9UPPyfU5X3qkqdOCRnO26KyexlxOsr3+B14wDc/bK7O6XOdw
rjbzAldk0Tj8WXLM1BId7mZZhWtANY28wftup1v43CRDH76ol0aw0sqowFyb
W7DMpu2WHoMU+w3mbXpaXOaDY30JLgwaokL6EzKmc5Whht3DEQx+5pVV3XfS
4GovtLTiGAPFN6Oe+ud4XN7wKclyGQP/lTELCZnEFUlGH28cZ1D/bpP1FekH
aY5h+tA6UwHuqg+f0D8xKHVLsXw/ESBAx86JlT9N9Ep7S3Tf70kS4NPNZ9T+
sGuTetfM6JBcEWBv8ImHiX7ZTq88UrU0uwXItwq12x7U5hSR4mpvyRei6m5C
5+8jvU41DXwPbWchLuZVJhaVDTq9kZisr9smRLHHvQRvyw9OX+tUVuQVCiEz
ez6wXzgJqTvX91B3hXAMLp1XHclDafJSk7h3QthwPlFvNTWwt3CRwyGxCB+O
Pz7wXdZ0ZNMXSht8RbBYSGUueTILGb6hPFmCCNFRZpt51tpYk267OvOSCDPK
538saOdjNHxFFqUQIbVIZ6ApTgetDS0GDzTFWMUOfdznMgd7mD2ZS+3EcEm9
OrVLqYeZphvHdoaIER8fyd6W6CPY2t8rLFeMoTXan/sHGaCkN+OYVa0YXk/z
byRNNgJfz912Xr8Yr9XmetncMYb7Cr2qq0YEf32s2XbHlkIYJbQ0XUUQUa9c
ftWYhs+4S84XsQSR4WnO0WoMtnop7Y+WE2QnGt/8RM4gcm3Ts6pHBLJN5o1B
VQI4eNeFpqmzMD1v7puWLsSbidHXlotYjLSUTk8KFOEpc5kbDmTR03TIcznE
SHFLjI7KYDG4IFpLPpuAMcgm9ddYFDeP1ybqEahnv5+Wfp2FW0RytLUBgZIX
0uNdzeLo5fz2XIbgYr9deXsNC1eH6nP+lgSOtS9Mh6QsDruruSk8Cb76fskC
YT2LxZsPxitTCXKrFNb7HrDoUJ9tXZhBsMvMXddNzuJg2ZG/V+UQBB6XDM18
yELxssKzIp9AmJZypeARi6TtMsPwMoLKDYtsLz1h8ThO59pgnWq3etLSDgWL
eKoo7NRtAknsa4PTz1T7a0WUXyNB3j9fjm3pYBGnviS5Sk6wodXit+HnLMTp
Ab6xnQSvKlodtbpYNJt3aszvJmimXU0eqThWFlrTplT58n6ZKOxmcZ//o2DZ
AEFsfCIn7mER8yuvpW9Q5RvqK3mlYsYv42DxCIFzqH9CpZJF47Cu49oxAtGz
2o0xL1n8/18EmuvMnR17WfwLDjcoGQ==
      "]]}}},
  AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948],
  Axes->{True, True},
  AxesLabel->{None, None},
  AxesOrigin->{0, 7.8500000000000005`},
  DisplayFunction->Identity,
  Frame->{{False, False}, {False, False}},
  FrameLabel->{{None, None}, {None, None}},
  FrameTicks->{{Automatic, Automatic}, {Automatic, Automatic}},
  GridLines->{None, None},
  GridLinesStyle->Directive[
    GrayLevel[0.5, 0.4]],
  Method->{},
  PlotRange->{{-8., 6.}, {8., 11.}},
  PlotRangeClipping->True,
  PlotRangePadding->{{
     Scaled[0.02], 
     Scaled[0.02]}, {
     Scaled[0.05], 
     Scaled[0.05]}},
  Ticks->{Automatic, Automatic}]], "Output",
 CellChangeTimes->{{3.7009387004624987`*^9, 3.7009387859986973`*^9}, 
   3.70093884028319*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"9.945945945945946`", "\[VeryThinSpace]", "+", 
   RowBox[{"0.20945945945945946`", " ", "x"}]}], "/.", 
  RowBox[{"x", "\[Rule]", 
   RowBox[{"{", 
    RowBox[{
     RowBox[{"-", "8"}], ",", 
     RowBox[{"-", "2"}], ",", "6"}], "}"}]}]}]], "Input",
 CellChangeTimes->{{3.700938922929449*^9, 3.700938944951849*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "8.27027027027027`", ",", "9.527027027027026`", ",", 
   "11.202702702702702`"}], "}"}]], "Output",
 CellChangeTimes->{{3.7009389350961227`*^9, 3.70093894563382*^9}}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{"{", 
   RowBox[{
   "8.27027027027027`", ",", "9.527027027027026`", ",", 
    "11.202702702702702`"}], "}"}], "-", 
  RowBox[{"{", 
   RowBox[{"8", ",", "10", ",", "11"}], "}"}]}]], "Input",
 CellChangeTimes->{{3.700938970383787*^9, 3.700938986605794*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{"0.2702702702702702`", ",", 
   RowBox[{"-", "0.4729729729729737`"}], ",", "0.20270270270270174`"}], 
  "}"}]], "Output",
 CellChangeTimes->{3.700938987490899*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{
  RowBox[{
   RowBox[{
    FractionBox["1", 
     RowBox[{"2", "\[Pi]"}]], 
    RowBox[{"Exp", "[", 
     RowBox[{"-", 
      SuperscriptBox[
       RowBox[{"(", "\[Delta]", ")"}], "2"]}], "]"}], 
    FractionBox["1", 
     RowBox[{"2", "\[Pi]"}]]}], "/.", 
   RowBox[{"\[Delta]", "\[Rule]", 
    RowBox[{"{", 
     RowBox[{"0.2702702702702702`", ",", 
      RowBox[{"-", "0.4729729729729737`"}], ",", "0.20270270270270174`"}], 
     "}"}]}]}], "//", "N"}]], "Input",
 CellChangeTimes->{{3.700939014484818*^9, 3.700939063106917*^9}, {
  3.70094105945088*^9, 3.700941077062516*^9}, {3.700952447498179*^9, 
  3.700952459948997*^9}}],

Cell[BoxData[
 RowBox[{"{", 
  RowBox[{
  "0.02354598051605746`", ",", "0.02025289432564488`", ",", 
   "0.0243106070806057`"}], "}"}]], "Output",
 CellChangeTimes->{{3.700939015330624*^9, 3.700939043432007*^9}, 
   3.700941077991856*^9, 3.7009524661724*^9}]
}, Open  ]],

Cell[BoxData[
 RowBox[{
  RowBox[{"L0", "=", 
   RowBox[{"{", 
    RowBox[{
    "0.02354598051605746`", ",", "0.02025289432564488`", ",", 
     "0.0243106070806057`"}], "}"}]}], ";"}]], "Input",
 CellChangeTimes->{{3.7009408291011*^9, 3.700940861149872*^9}, 
   3.70094108631638*^9}],

Cell[CellGroupData[{

Cell[BoxData[
 RowBox[{"evidence0", "=", 
  RowBox[{"Product", "[", 
   RowBox[{
    RowBox[{"L0", "[", 
     RowBox[{"[", "i", "]"}], "]"}], ",", 
    RowBox[{"{", 
     RowBox[{"i", ",", "1", ",", "3"}], "}"}]}], "]"}]}]], "Input",
 CellChangeTimes->{{3.700940737181302*^9, 3.7009407878375273`*^9}, {
  3.700940863266679*^9, 3.70094089621912*^9}}],

Cell[BoxData["0.000011593102644668878`"], "Output",
 CellChangeTimes->{{3.700940769412107*^9, 3.7009407962889557`*^9}, 
   3.700940863923339*^9, 3.700940896531225*^9, 3.700941089632154*^9, 
   3.700952477526828*^9}]
}, Open  ]]
}, Open  ]],

Cell[CellGroupData[{

Cell["Ratio of evidences", "Subchapter",
 CellChangeTimes->{{3.7009312250855007`*^9, 3.700931237388666*^9}, {
  3.7009394389791937`*^9, 3.7009394482944527`*^9}, {3.7009394791270742`*^9, 
  3.700939501365468*^9}, {3.700939535545888*^9, 3.700939555958934*^9}, {
  3.70094091183328*^9, 3.700940914784554*^9}}],

Cell[CellGroupData[{

Cell[BoxData[
 FractionBox["evidence0", "evidence1"]], "Input",
 CellChangeTimes->{{3.7009409704457207`*^9, 3.7009410093873167`*^9}}],

Cell[BoxData["4.534154481613225`*^27"], "Output",
 CellChangeTimes->{3.700941009891736*^9, 3.700941092143371*^9, 
  3.7009524796223907`*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 FractionBox["evidence0", "evidence2"]], "Input",
 CellChangeTimes->{{3.700941012499475*^9, 3.700941020106491*^9}}],

Cell[BoxData["743.5995347318301`"], "Output",
 CellChangeTimes->{3.700941020384047*^9, 3.700941093536951*^9, 
  3.700952481880179*^9}]
}, Open  ]],

Cell[CellGroupData[{

Cell[BoxData[
 FractionBox["evidence2", "evidence1"]], "Input",
 CellChangeTimes->{{3.700941037881441*^9, 3.700941041679965*^9}}],

Cell[BoxData["6.097575737790652`*^24"], "Output",
 CellChangeTimes->{{3.7009410222411337`*^9, 3.700941042130378*^9}, 
   3.70094109470107*^9, 3.700952483219792*^9}]
}, Open  ]]
}, Open  ]]
},
WindowSize->{808, 698},
WindowMargins->{{58, Automatic}, {Automatic, 44}},
PrintingCopies->1,
PrintingPageRange->{1, Automatic},
FrontEndVersion->"10.1 for Mac OS X x86 (32-bit, 64-bit Kernel) (March 23, \
2015)",
StyleDefinitions->"Default.nb"
]
(* End of Notebook Content *)

(* Internal cache information *)
(*CellTagsOutline
CellTagsIndex->{}
*)
(*CellTagsIndex
CellTagsIndex->{}
*)
(*NotebookFileOutline
Notebook[{
Cell[CellGroupData[{
Cell[580, 22, 439, 13, 64, "Subchapter"],
Cell[1022, 37, 183, 3, 28, "Input"],
Cell[1208, 42, 378, 13, 28, "Input"],
Cell[1589, 57, 285, 7, 28, "Input"],
Cell[1877, 66, 987, 29, 56, "Input"],
Cell[CellGroupData[{
Cell[2889, 99, 416, 13, 68, InheritFromParent],
Cell[3308, 114, 223, 5, 69, "Output"]
}, Open  ]],
Cell[3546, 122, 345, 8, 28, InheritFromParent],
Cell[CellGroupData[{
Cell[3916, 134, 448, 10, 28, "Input"],
Cell[4367, 146, 164, 2, 32, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[4580, 154, 439, 13, 64, "Subchapter"],
Cell[5022, 169, 281, 5, 28, "Input"],
Cell[5306, 176, 118, 4, 28, "Input"],
Cell[CellGroupData[{
Cell[5449, 184, 910, 29, 56, "Input"],
Cell[6362, 215, 683, 21, 66, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[7082, 241, 1162, 36, 56, "Input"],
Cell[8247, 279, 861, 29, 92, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[9145, 313, 1118, 39, 94, "Input"],
Cell[10266, 354, 218, 5, 28, "Output"]
}, Open  ]],
Cell[10499, 362, 342, 8, 28, "Input"],
Cell[CellGroupData[{
Cell[10866, 374, 401, 10, 28, "Input"],
Cell[11270, 386, 220, 3, 32, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[11539, 395, 288, 3, 62, "Subchapter"],
Cell[CellGroupData[{
Cell[11852, 402, 422, 13, 28, "Input"],
Cell[12277, 417, 316, 11, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[12630, 433, 380, 9, 28, "Input"],
Cell[13013, 444, 988, 25, 254, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[14038, 474, 599, 16, 46, "Input"],
Cell[14640, 492, 1011, 30, 265, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[15688, 527, 198, 4, 28, "Input"],
Cell[15889, 533, 783, 23, 42, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16709, 561, 72, 1, 28, "Input"],
Cell[16784, 564, 171, 3, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[16992, 572, 388, 10, 28, "Input"],
Cell[17383, 584, 2704, 56, 245, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[20124, 645, 148, 3, 28, "Input"],
Cell[20275, 650, 2768, 57, 253, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23080, 712, 349, 9, 28, "Input"],
Cell[23432, 723, 208, 5, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[23677, 733, 291, 8, 28, "Input"],
Cell[23971, 743, 202, 5, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[24210, 753, 654, 19, 69, "Input"],
Cell[24867, 774, 258, 6, 28, "Output"]
}, Open  ]],
Cell[25140, 783, 283, 8, 28, "Input"],
Cell[CellGroupData[{
Cell[25448, 795, 349, 9, 28, "Input"],
Cell[25800, 806, 215, 3, 28, "Output"]
}, Open  ]]
}, Open  ]],
Cell[CellGroupData[{
Cell[26064, 815, 306, 4, 62, "Subchapter"],
Cell[CellGroupData[{
Cell[26395, 823, 133, 2, 49, "Input"],
Cell[26531, 827, 140, 2, 32, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[26708, 834, 129, 2, 49, "Input"],
Cell[26840, 838, 134, 2, 28, "Output"]
}, Open  ]],
Cell[CellGroupData[{
Cell[27011, 845, 129, 2, 49, "Input"],
Cell[27143, 849, 164, 2, 66, "Output"]
}, Open  ]]
}, Open  ]]
}
]
*)

(* End of internal cache information *)