(* Content-type: application/vnd.wolfram.mathematica *) (*** Wolfram Notebook File ***) (* http://www.wolfram.com/nb *) (* CreatedBy='Mathematica 9.0' *) (*CacheID: 234*) (* Internal cache information: NotebookFileLineBreakTest NotebookFileLineBreakTest NotebookDataPosition[ 157, 7] NotebookDataLength[ 98455, 2471] NotebookOptionsPosition[ 93053, 2279] NotebookOutlinePosition[ 93595, 2300] CellTagsIndexPosition[ 93552, 2297] WindowFrame->Normal*) (* Beginning of Notebook Content *) Notebook[{ Cell[CellGroupData[{ Cell[TextData[StyleBox["O Limite M\[IAcute]nimo da Vari\[AHat]ncia", \ "Subtitle"]], "Subsubtitle", CellChangeTimes->{{3.591944385734125*^9, 3.5919443950627675`*^9}}], Cell[CellGroupData[{ Cell["C\[AAcute]lculo te\[OAcute]rico para uma f.d.p. de Cauchy", "Section", CellChangeTimes->{{3.5919444186734023`*^9, 3.591944436299327*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"f", "=", RowBox[{"PDF", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"\[Theta]", ",", "1"}], "]"}], ",", "x"}], "]"}]}]], "Input", CellChangeTimes->{{3.5919444423777685`*^9, 3.5919444719575033`*^9}}], Cell[BoxData[ FractionBox["1", RowBox[{"\[Pi]", " ", RowBox[{"(", RowBox[{"1", "+", SuperscriptBox[ RowBox[{"(", RowBox[{"x", "-", "\[Theta]"}], ")"}], "2"]}], ")"}]}]]], "Output", CellChangeTimes->{3.6721383882198124`*^9, 3.6730786330336027`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Integrate", "[", RowBox[{"f", ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "Infinity"}], ",", "Infinity"}], "}"}], ",", RowBox[{"Assumptions", "\[Rule]", RowBox[{"\[Theta]", "\[Element]", "Reals"}]}]}], "]"}], " ", RowBox[{"(*", " ", RowBox[{ "usando", " ", "a", " ", "propriedade", " ", "da", " ", "normaliza\[CCedilla]\[ATilde]o", " ", "para", " ", "verificar", " ", "que", " ", "achamos", " ", "como", " ", "integrar", " ", "essa", " ", "fun\[CCedilla]\[ATilde]o"}], " ", "*)"}]}]], "Input", CellChangeTimes->{{3.5919452667412095`*^9, 3.5919453188506556`*^9}, { 3.5919475156238165`*^9, 3.5919475547349033`*^9}}], Cell[BoxData["1"], "Output", CellChangeTimes->{3.6721383885200777`*^9, 3.6730786337354836`*^9}] }, Open ]], Cell["\<\ LMV \[EAcute] o inverso e oposto do valor esperado da segunda derivada - come\ \[CCedilla]amos derivando\ \>", "Text", CellChangeTimes->{{3.5919445095063562`*^9, 3.591944535804658*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d2f", "=", RowBox[{"D", "[", RowBox[{ RowBox[{"Log", "[", "f", "]"}], ",", RowBox[{"{", RowBox[{"\[Theta]", ",", "2"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5919445435706773`*^9, 3.5919446005893803`*^9}, { 3.591945390056634*^9, 3.5919453954568143`*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"-", FractionBox["2", RowBox[{"1", "+", SuperscriptBox[ RowBox[{"(", RowBox[{"x", "-", "\[Theta]"}], ")"}], "2"]}]]}], "+", FractionBox[ RowBox[{"4", " ", SuperscriptBox[ RowBox[{"(", RowBox[{"x", "-", "\[Theta]"}], ")"}], "2"]}], SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", SuperscriptBox[ RowBox[{"(", RowBox[{"x", "-", "\[Theta]"}], ")"}], "2"]}], ")"}], "2"]]}]], "Output", CellChangeTimes->{3.6721383885200777`*^9, 3.6730786337354836`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"d2f", "=", RowBox[{"FullSimplify", "[", "d2f", "]"}]}]], "Input", CellChangeTimes->{{3.5919446384038963`*^9, 3.5919446591706066`*^9}, { 3.5919475762395015`*^9, 3.591947581519679*^9}}], Cell[BoxData[ FractionBox[ RowBox[{"2", " ", RowBox[{"(", RowBox[{ RowBox[{"-", "1"}], "+", SuperscriptBox[ RowBox[{"(", RowBox[{"x", "-", "\[Theta]"}], ")"}], "2"]}], ")"}]}], SuperscriptBox[ RowBox[{"(", RowBox[{"1", "+", SuperscriptBox[ RowBox[{"(", RowBox[{"x", "-", "\[Theta]"}], ")"}], "2"]}], ")"}], "2"]]], "Output", CellChangeTimes->{3.6721383885392113`*^9, 3.6730786337888746`*^9}] }, Open ]], Cell[CellGroupData[{ Cell["\<\ A fun\[CCedilla]\[ATilde]o verossimilhan\[CCedilla]a \[EAcute] o produto de n \ fun\[CCedilla]\[OTilde]es de Cauchy e a derivada do LogL d\[AAcute] uma soma \ de n termos. Assim, terminamos com a soma de n termos identicos...\ \>", "Subsubsection", CellChangeTimes->{{3.5919476212187815`*^9, 3.5919477105052614`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"lmv", "=", RowBox[{ RowBox[{"-", "1"}], "/", RowBox[{"(", RowBox[{"n", " ", RowBox[{"Integrate", "[", RowBox[{ RowBox[{"f", " ", "d2f"}], ",", RowBox[{"{", RowBox[{"x", ",", RowBox[{"-", "Infinity"}], ",", "Infinity"}], "}"}], ",", RowBox[{"Assumptions", "\[Rule]", RowBox[{"\[Theta]", "\[Element]", "Reals"}]}]}], "]"}]}], ")"}]}]}]], "Input", CellChangeTimes->{{3.591944672483818*^9, 3.5919447883806686`*^9}, { 3.5919448646191144`*^9, 3.591944901496089*^9}, {3.5919453542836404`*^9, 3.5919453645659266`*^9}}], Cell[BoxData[ FractionBox["2", "n"]], "Output", CellChangeTimes->{3.67213838864156*^9, 3.67307863413657*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Compara\[CCedilla]\[ATilde]o com dois estimadores diferentes...", \ "Subsection", CellChangeTimes->{{3.591949296203309*^9, 3.5919493079876366`*^9}}], Cell[CellGroupData[{ Cell["Determina\[CCedilla]\[ATilde]o da vari\[AHat]ncia da mediana como \ estimativa de \[Theta]", "Subsubsection", CellChangeTimes->{{3.5919477653967857`*^9, 3.5919478355547094`*^9}, 3.5919478826900196`*^9, {3.591947924333207*^9, 3.5919479442197566`*^9}}], Cell["simula\[CCedilla]\[ATilde]o...", "Text", CellChangeTimes->{{3.591947954018005*^9, 3.591947957090119*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"u", "=", RowBox[{"RandomVariate", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"0", ",", "1"}], "]"}], ",", "10"}], "]"}]}]], "Input", CellChangeTimes->{{3.5919477355539293`*^9, 3.591947751434371*^9}, { 3.591947842128892*^9, 3.591947842746909*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{"1.2755684407166412`", ",", RowBox[{"-", "2.6955008252727843`"}], ",", RowBox[{"-", "34.178651498768104`"}], ",", "1.21774628679684`", ",", RowBox[{"-", "1.4497025629608569`"}], ",", RowBox[{"-", "0.34661459149979806`"}], ",", "0.2595230312661887`", ",", RowBox[{"-", "0.8451914411512067`"}], ",", "0.7083983528470047`", ",", "0.34820148963116565`"}], "}"}]], "Output", CellChangeTimes->{3.672138388657576*^9, 3.6730786342680936`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Median", "[", "u", "]"}]], "Input", CellChangeTimes->{{3.591947858317342*^9, 3.5919478631295033`*^9}}], Cell[BoxData[ RowBox[{"-", "0.04354578011680468`"}]], "Output", CellChangeTimes->{3.672138388657576*^9, 3.6730786342837205`*^9}] }, Open ]], Cell["\<\ repetir a simula\[CCedilla]\[ATilde]o muitas vezes (1000) e determinar a \ variancia da estimativa, para comparar com o LMV\ \>", "Text", CellChangeTimes->{{3.591947903832608*^9, 3.591947914600934*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"q", "=", RowBox[{"Table", "[", RowBox[{ RowBox[{"Median", "[", RowBox[{"RandomVariate", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"0", ",", "1"}], "]"}], ",", "10"}], "]"}], "]"}], ",", RowBox[{"{", "1000", "}"}]}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.5919448452118006`*^9, 3.591944854868598*^9}, { 3.5919454533844986`*^9, 3.591945537208825*^9}, {3.5919478703556767`*^9, 3.5919478745557947`*^9}, {3.591947978874693*^9, 3.5919480056894655`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Variance", "[", "q", "]"}]], "Input", CellChangeTimes->{{3.591947966150342*^9, 3.5919479699704456`*^9}}], Cell[BoxData["0.330393863814953`"], "Output", CellChangeTimes->{3.6721383886888275`*^9, 3.67307863432104*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Repetir o c\[AAcute]lculo para v\[AAcute]rios valores de n, pequenos", \ "Subsubsection", CellChangeTimes->{{3.591948037266345*^9, 3.591948050966726*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"varMediana", "=", RowBox[{"Table", "[", " ", RowBox[{ RowBox[{"{", RowBox[{"n", ",", RowBox[{"Variance", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"Median", "[", RowBox[{"RandomVariate", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"0", ",", "1"}], "]"}], ",", "n"}], "]"}], "]"}], ",", RowBox[{"{", "10000", "}"}]}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"n", ",", "3", ",", "30"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.591945605336722*^9, 3.5919456637763453`*^9}, 3.591948060570964*^9}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"3", ",", "15.248067842196383`"}], "}"}], ",", RowBox[{"{", RowBox[{"4", ",", "5.65177144768403`"}], "}"}], ",", RowBox[{"{", RowBox[{"5", ",", "1.0524053829353786`"}], "}"}], ",", RowBox[{"{", RowBox[{"6", ",", "0.7572138557081913`"}], "}"}], ",", RowBox[{"{", RowBox[{"7", ",", "0.6233398114601734`"}], "}"}], ",", RowBox[{"{", RowBox[{"8", ",", "0.5008992578634011`"}], "}"}], ",", RowBox[{"{", RowBox[{"9", ",", "0.4184555946714506`"}], "}"}], ",", RowBox[{"{", RowBox[{"10", ",", "0.34014384874557785`"}], "}"}], ",", RowBox[{"{", RowBox[{"11", ",", "0.3051193752691125`"}], "}"}], ",", RowBox[{"{", RowBox[{"12", ",", "0.2576402015706726`"}], "}"}], ",", RowBox[{"{", RowBox[{"13", ",", "0.23550406943397914`"}], "}"}], ",", RowBox[{"{", RowBox[{"14", ",", "0.21736907328135396`"}], "}"}], ",", RowBox[{"{", RowBox[{"15", ",", "0.20847790829420637`"}], "}"}], ",", RowBox[{"{", RowBox[{"16", ",", "0.1843796618604669`"}], "}"}], ",", RowBox[{"{", RowBox[{"17", ",", "0.1707525546301834`"}], "}"}], ",", RowBox[{"{", RowBox[{"18", ",", "0.1537587517093262`"}], "}"}], ",", RowBox[{"{", RowBox[{"19", ",", "0.1503904588087344`"}], "}"}], ",", RowBox[{"{", RowBox[{"20", ",", "0.1403996926181372`"}], "}"}], ",", RowBox[{"{", RowBox[{"21", ",", "0.13661065515622708`"}], "}"}], ",", RowBox[{"{", RowBox[{"22", ",", "0.12190908886063628`"}], "}"}], ",", RowBox[{"{", RowBox[{"23", ",", "0.12164682183807665`"}], "}"}], ",", RowBox[{"{", RowBox[{"24", ",", "0.11353501395391281`"}], "}"}], ",", RowBox[{"{", RowBox[{"25", ",", "0.11178777726311867`"}], "}"}], ",", RowBox[{"{", RowBox[{"26", ",", "0.1019344024348612`"}], "}"}], ",", RowBox[{"{", RowBox[{"27", ",", "0.10243304195579689`"}], "}"}], ",", RowBox[{"{", RowBox[{"28", ",", "0.09438946800520782`"}], "}"}], ",", RowBox[{"{", RowBox[{"29", ",", "0.09499870525068571`"}], "}"}], ",", RowBox[{"{", RowBox[{"30", ",", "0.08981960019280849`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.67213839685714*^9, 3.673078642117547*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"{", RowBox[{ RowBox[{"ListPlot", "[", RowBox[{"varMediana", ",", RowBox[{"PlotRange", "\[Rule]", "All"}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"BaseStyle", "\[Rule]", RowBox[{"{", RowBox[{"FontSize", "\[Rule]", "14"}], "}"}]}]}], "]"}], ",", RowBox[{"Plot", "[", RowBox[{ RowBox[{"2", "/", "n"}], ",", RowBox[{"{", RowBox[{"n", ",", "3", ",", "30"}], "}"}]}], "]"}]}], "}"}], "]"}]], "Input", CellChangeTimes->{{3.5919457113176656`*^9, 3.591945745994629*^9}, { 3.5919481018441105`*^9, 3.59194812620079*^9}, {3.591949375429513*^9, 3.5919494109404984`*^9}, {3.67307979810271*^9, 3.6730798424302683`*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {{}, {RGBColor[0.368417, 0.506779, 0.709798], PointSize[ 0.012833333333333334`], AbsoluteThickness[1.6], StyleBox[PointBox[CompressedData[" 1:eJxTTMoPSmViYGCQAWIQDQEcDsU3uvYz1es5QPgCDl+EC39mzhaD8kUcuCvu si2/9sEewpdwUOPw7pWwfgHlyzgotAqEpn16DOUrOGxbmvwwlv0BlK/kMNG7 5ufP47egfBWH2fyTt7w6chXKV3O4xtGVJdx+GcrXcChK+BGsW3sBytdy+D9H Nei/yjkoX8fh373iT/svnYby9RyYPrMXp68/BeUbOByru7flwMzjUL6hwyft N3YWt49C+UYOgrEcT2PXHYbyjR02JknP/mcP45s4bH8h8mfuj4NQvqmDrb64 TWkNjG/m8OTGXtl80/1QvrnDp68ljPYqML6Fw1zOR8ULBfZC+ZYOIRwF6+Xn 7YHyrRwM20/XxEvsgvKtHabctXLltYTxbRxEA7XrXmjugPJtHWxOm525Ggjj 2zmY6Mddzfq3zR4AESJ7kw== "]], FontSize->14]}, {}}, {}}, {{{}, {}, {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[ 1.], LineBox[CompressedData[" 1:eJwVzns40wscgHGPjDSVoa0RJ9dE5n738/0m1FQmnS4TqumCSJ1SilJuKSUl ndJJ6DLmUr+xaElodFEUI5F1jxKdUm4jp/PH+3z+ffV5Uf6bFRUUFIJ/97+7 XFIsVk2oozRw+p7lhjLI1Uv+IDpAw5wCalyBVzcoVyqphmRooMenvZS0M13g ObpxzZNzGhhAuT9TPbULDrvcveqYp4G9OrEbcmO7YLxqn4eaUANLo83PveV1 wffagVhRqwbC9xVVfdZd8KqxY0CFoYn8L2KqUksnVL4qlhZf1ETribLjfcxO CFNZlTsk0MKR3Rbx7fUdkLIlN4VRroWZxbren8UdcKW+L8Lpjhaa2+YIVG50 wOvEBOf9zVroNCFkrs/ugLUKN1omBrVQmyq0iY/qAO/RaRQKMQvNFl8W7NPp AIO+mjDN5ll48qywXHn/c+hqWmBrNUhH+dKMIntOO5RTQoxFcjqyWywF/3q3 Q7rbeYaLEgMb6eP3Sfd2QAFlfNEsBi6i/vLyZrXDlWSZZI0jAy30qBkV09th u1vG6vhYBnamL78ketIGCoLBfU1TZmNEq17R4Io2MEkW10RoMXHnP8xE23Ap vH4h8ynQZWLLX02BQp4UzrOmtL0zYWLQis+VDuukoNax9FOAMxNVFhu95iyT wjezbho7mIk8Y7NF91hSqHo6yTMuYOLefIpb4Y9WWDnHiyJz1cZoj68lPkmt cKismc0J0cG0EEv1OrIFdvyyN9zG10VtLZmDmuczcHHyHz4q+QNTShwF33Wf gurP/SF8FX3sMlvVe0a9CXj1/CVpTAMcnh25OZ/9GJSNJtqoJoZoOrlxznjV IxC60zeuZxuhPPPUCUXiITT6HQuz4xrj8Y3RrVfl96FPry7qZbwJjhnmzsio bYDVsQZ6+nnzMG8NC8cF9cBq7mb3lJri8gu8W2bVEpj+KnDXodb5GHD8gfr8 1nvA/eigfeG9Gf65PaeyTfUelBRKq1PVFuCzmheVI1gHy+TZVE1LC6xJ0Pvm nlYLaz2r2vTcWeib3yNSHKiBbey8ixMhltiQ1ClL4tTAUP6ObyxXK5zTbXRz RdpdWNdrFFpSa4WBMf1cg/RqkFBvB2WttEa237Cf9PgdeC/4RavussZJi2qL I0VVoGJq5OUXZYMt1nId5Ybb4Pcl6LGVgi0miHKOLPwuBs4CLw3/I7aowedv euwoBn+fGc9fMu0w8suRpMKYW0B9nxRQkG+H/b2ajObmSvgsmWlTZ2qP9NhS 85M2lfAkMS7VS2yPTrQizuGMCthSffaiqocDZkd5M+oVKmDK89PWjm0O6MoN kByMuAmrt7klXQ12RJOHcTG5vSJo9h4z431wRNdDIxGzI0WwVCHbL3WvEy6z 7zGc018OB+177L5OOmGg91HVr/vLQbSEfu5upjO6f+wPzFQqB9plLkdm5IJv 6OJP0VllYLjPp54odsGCALupNJsycLbQFSuDK1aI/7blPxKCctDFSPUHrnhi XlLJz0ghjLRaR4f6uuGVxendB+hCGDAtkM9/7YbnL8XFWN0mIXu0EReEEnhg siiKFUaCe/Plk2nhBOZlHpM3bCXh7ZVY2ecIAjlZwbXrt5Bg7mceV7iTwFfP S8dOhZBQde1YhUksgde1wgxGg0joXslm6acTuGp8tkuLPwlzSxt06OUEcj1+ zL/uRoIkMSc8+iaBTZwbL3xdSQjl7rklrSSwZ6SwdsCZBFLJZE3mHQIFS+a6 WDqS4LEu+bT6fQKNR3MbhdYkbJrqqUrtJHDm7oWbJSYkTJXprA1/SWCW5qnT W41JKC4bvPZQ9vuvZ8nwNCMSfgZf9kx9R6BSWNUWf30SUkSKhyj9BJZzPF6+ 0yHBLK2zadNXAmljvHup2iQ0bRDqSr4R6CtOeGPBJIGuxrudMETgkGtaXQyd BPEb52lvRwiM7+Em680iIbiCxl0oJ9DgwOBBiSYJiic+8XMnCIRxn+vhGiTw ebVDk5ME+oWupNNoJPwHUMuuMg== "]]}}, {}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{ FormBox["\"n\"", TraditionalForm], FormBox["\"variance\"", TraditionalForm]}, AxesOrigin->{2.578125, 0}, BaseStyle->{FontSize -> 14}, 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]], ImagePadding->All, Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& )}}, PlotRange->{{3., 30.}, {0, 15.248067842196383`}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.02], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{ 3.672138396910536*^9, 3.67307864231809*^9, {3.67307982477041*^9, 3.67307984376007*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"{", RowBox[{ RowBox[{"ListPlot", "[", RowBox[{"varMediana", ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0.5"}], "}"}]}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"BaseStyle", "\[Rule]", RowBox[{"{", RowBox[{"FontSize", "\[Rule]", "14"}], "}"}]}]}], "]"}], ",", RowBox[{"Plot", "[", RowBox[{ RowBox[{"2", "/", "n"}], ",", RowBox[{"{", RowBox[{"n", ",", "3", ",", "30"}], "}"}]}], "]"}]}], "}"}], "]"}]], "Input", CellChangeTimes->{{3.5919494440894194`*^9, 3.591949449883581*^9}, 3.673079856895153*^9}], Cell[BoxData[ GraphicsBox[{{{}, {{}, {RGBColor[0.368417, 0.506779, 0.709798], PointSize[ 0.012833333333333334`], AbsoluteThickness[1.6], StyleBox[PointBox[CompressedData[" 1:eJxTTMoPSmViYGCQAWIQDQEcDsU3uvYz1es5QPgCDl+EC39mzhaD8kUcuCvu si2/9sEewpdwUOPw7pWwfgHlyzgotAqEpn16DOUrOGxbmvwwlv0BlK/kMNG7 5ufP47egfBWH2fyTt7w6chXKV3O4xtGVJdx+GcrXcChK+BGsW3sBytdy+D9H Nei/yjkoX8fh373iT/svnYby9RyYPrMXp68/BeUbOByru7flwMzjUL6hwyft N3YWt49C+UYOgrEcT2PXHYbyjR02JknP/mcP45s4bH8h8mfuj4NQvqmDrb64 TWkNjG/m8OTGXtl80/1QvrnDp68ljPYqML6Fw1zOR8ULBfZC+ZYOIRwF6+Xn 7YHyrRwM20/XxEvsgvKtHabctXLltYTxbRxEA7XrXmjugPJtHWxOm525Ggjj 2zmY6Mddzfq3zR4AESJ7kw== "]], FontSize->14]}, {}}, {}}, {{{}, {}, {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[ 1.], LineBox[CompressedData[" 1:eJwVzns40wscgHGPjDSVoa0RJ9dE5n738/0m1FQmnS4TqumCSJ1SilJuKSUl ndJJ6DLmUr+xaElodFEUI5F1jxKdUm4jp/PH+3z+ffV5Uf6bFRUUFIJ/97+7 XFIsVk2oozRw+p7lhjLI1Uv+IDpAw5wCalyBVzcoVyqphmRooMenvZS0M13g ObpxzZNzGhhAuT9TPbULDrvcveqYp4G9OrEbcmO7YLxqn4eaUANLo83PveV1 wffagVhRqwbC9xVVfdZd8KqxY0CFoYn8L2KqUksnVL4qlhZf1ETribLjfcxO CFNZlTsk0MKR3Rbx7fUdkLIlN4VRroWZxbren8UdcKW+L8Lpjhaa2+YIVG50 wOvEBOf9zVroNCFkrs/ugLUKN1omBrVQmyq0iY/qAO/RaRQKMQvNFl8W7NPp AIO+mjDN5ll48qywXHn/c+hqWmBrNUhH+dKMIntOO5RTQoxFcjqyWywF/3q3 Q7rbeYaLEgMb6eP3Sfd2QAFlfNEsBi6i/vLyZrXDlWSZZI0jAy30qBkV09th u1vG6vhYBnamL78ketIGCoLBfU1TZmNEq17R4Io2MEkW10RoMXHnP8xE23Ap vH4h8ynQZWLLX02BQp4UzrOmtL0zYWLQis+VDuukoNax9FOAMxNVFhu95iyT wjezbho7mIk8Y7NF91hSqHo6yTMuYOLefIpb4Y9WWDnHiyJz1cZoj68lPkmt cKismc0J0cG0EEv1OrIFdvyyN9zG10VtLZmDmuczcHHyHz4q+QNTShwF33Wf gurP/SF8FX3sMlvVe0a9CXj1/CVpTAMcnh25OZ/9GJSNJtqoJoZoOrlxznjV IxC60zeuZxuhPPPUCUXiITT6HQuz4xrj8Y3RrVfl96FPry7qZbwJjhnmzsio bYDVsQZ6+nnzMG8NC8cF9cBq7mb3lJri8gu8W2bVEpj+KnDXodb5GHD8gfr8 1nvA/eigfeG9Gf65PaeyTfUelBRKq1PVFuCzmheVI1gHy+TZVE1LC6xJ0Pvm nlYLaz2r2vTcWeib3yNSHKiBbey8ixMhltiQ1ClL4tTAUP6ObyxXK5zTbXRz RdpdWNdrFFpSa4WBMf1cg/RqkFBvB2WttEa237Cf9PgdeC/4RavussZJi2qL I0VVoGJq5OUXZYMt1nId5Ybb4Pcl6LGVgi0miHKOLPwuBs4CLw3/I7aowedv euwoBn+fGc9fMu0w8suRpMKYW0B9nxRQkG+H/b2ajObmSvgsmWlTZ2qP9NhS 85M2lfAkMS7VS2yPTrQizuGMCthSffaiqocDZkd5M+oVKmDK89PWjm0O6MoN kByMuAmrt7klXQ12RJOHcTG5vSJo9h4z431wRNdDIxGzI0WwVCHbL3WvEy6z 7zGc018OB+177L5OOmGg91HVr/vLQbSEfu5upjO6f+wPzFQqB9plLkdm5IJv 6OJP0VllYLjPp54odsGCALupNJsycLbQFSuDK1aI/7blPxKCctDFSPUHrnhi XlLJz0ghjLRaR4f6uuGVxendB+hCGDAtkM9/7YbnL8XFWN0mIXu0EReEEnhg siiKFUaCe/Plk2nhBOZlHpM3bCXh7ZVY2ecIAjlZwbXrt5Bg7mceV7iTwFfP S8dOhZBQde1YhUksgde1wgxGg0joXslm6acTuGp8tkuLPwlzSxt06OUEcj1+ zL/uRoIkMSc8+iaBTZwbL3xdSQjl7rklrSSwZ6SwdsCZBFLJZE3mHQIFS+a6 WDqS4LEu+bT6fQKNR3MbhdYkbJrqqUrtJHDm7oWbJSYkTJXprA1/SWCW5qnT W41JKC4bvPZQ9vuvZ8nwNCMSfgZf9kx9R6BSWNUWf30SUkSKhyj9BJZzPF6+ 0yHBLK2zadNXAmljvHup2iQ0bRDqSr4R6CtOeGPBJIGuxrudMETgkGtaXQyd BPEb52lvRwiM7+Em680iIbiCxl0oJ9DgwOBBiSYJiic+8XMnCIRxn+vhGiTw ebVDk5ME+oWupNNoJPwHUMuuMg== "]]}}, {}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{ FormBox["\"n\"", TraditionalForm], FormBox["\"variance\"", TraditionalForm]}, AxesOrigin->{0, 0}, BaseStyle->{FontSize -> 14}, 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]], ImagePadding->All, Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& )}}, PlotRange->{{0, 30.}, {0, 0.5}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, {0, 0}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.6721383969417887`*^9, 3.673078642340226*^9, 3.6730798591017714`*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Para n grande, h\[AAcute] outra forma de estimar mais precisa que a mediana\ \>", "Subsubsection", CellChangeTimes->{{3.591948180028284*^9, 3.59194819692078*^9}}], Cell["\<\ Procedimento: 1) O conjunto de dados \[EAcute] ordenado - como estamos \ simulando, encadeamos as duas opera\[CCedilla]\[OTilde]es\ \>", "Text", CellChangeTimes->{{3.591948223069481*^9, 3.5919482268976154`*^9}, { 3.5919482590644865`*^9, 3.591948308413871*^9}, {3.5919493210139985`*^9, 3.591949331232295*^9}}], Cell[BoxData[ RowBox[{ RowBox[{"u", "=", RowBox[{"Sort", "[", RowBox[{"RandomVariate", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"0", ",", "1"}], "]"}], ",", "200"}], "]"}], "]"}]}], ";"}]], "Input", CellChangeTimes->{{3.591948241998007*^9, 3.59194824968622*^9}, { 3.5919483313644915`*^9, 3.591948331732501*^9}}], Cell["\<\ e 2) calcula-se a m\[EAcute]dia dos 24% dos dados em torno da mediana; Floor \ e Ceiling arredondam para Integer\ \>", "Text", CellChangeTimes->{{3.591948357437225*^9, 3.5919483753977156`*^9}, { 3.5919484381174707`*^9, 3.59194845478595*^9}, {3.5919493357744102`*^9, 3.591949337374454*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Mean", "[", RowBox[{"u", "[", RowBox[{"[", RowBox[{ RowBox[{"Floor", "[", RowBox[{"0.38", " ", "200"}], "]"}], ";;", RowBox[{"Ceiling", "[", RowBox[{"0.62", " ", "200"}], "]"}]}], "]"}], "]"}], " ", "]"}]], "Input", CellChangeTimes->{{3.591948399376381*^9, 3.591948428763198*^9}}], Cell[BoxData["0.11471871520597865`"], "Output", CellChangeTimes->{3.6721383969574413`*^9, 3.673078642355854*^9}] }, Open ]], Cell["\<\ D\[AAcute] para fazer em uma instru\[CCedilla]\[ATilde]o s\[OAcute], embora \ fique meio ileg\[IAcute]vel\ \>", "Text", CellChangeTimes->{{3.5919484900169125`*^9, 3.5919485037233133`*^9}, { 3.591948535770177*^9, 3.5919485424323874`*^9}, {3.591948761304441*^9, 3.591948762136469*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Mean", "[", RowBox[{ RowBox[{"Sort", "[", RowBox[{"RandomVariate", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"0", ",", "1"}], "]"}], ",", "200"}], "]"}], "]"}], "[", RowBox[{"[", RowBox[{ RowBox[{"Floor", "[", RowBox[{"0.38", " ", "200"}], "]"}], ";;", RowBox[{"Ceiling", "[", RowBox[{"0.62", " ", "200"}], "]"}]}], "]"}], "]"}], " ", "]"}]], "Input", CellChangeTimes->{3.591948520335743*^9}], Cell[BoxData[ RowBox[{"-", "0.1683609829427297`"}]], "Output", CellChangeTimes->{3.6721383969740458`*^9, 3.673078642355854*^9}] }, Open ]], Cell["\<\ repete a simula\[CCedilla]ao 100 vezes e calcula a variancia da estimativaa\ \>", "Text", CellChangeTimes->{{3.591948676382082*^9, 3.591948705988917*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Variance", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"Mean", "[", RowBox[{ RowBox[{"Sort", "[", RowBox[{"RandomVariate", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"0", ",", "1"}], "]"}], ",", "200"}], "]"}], "]"}], "[", RowBox[{"[", RowBox[{ RowBox[{"Floor", "[", RowBox[{"0.38", " ", "200"}], "]"}], ";;", RowBox[{"Ceiling", "[", RowBox[{"0.62", " ", "200"}], "]"}]}], "]"}], "]"}], " ", "]"}], ",", RowBox[{"{", "1000", "}"}]}], "]"}], "]"}]], "Input", CellChangeTimes->{{3.591946114678878*^9, 3.591946135443455*^9}, 3.5919486684698606`*^9}], Cell[BoxData["0.010845647137166068`"], "Output", CellChangeTimes->{3.672138397010805*^9, 3.673078642402731*^9}] }, Open ]], Cell["fazemos agora para diversos valores do n\[UAcute]mero de dados n", \ "Text", CellChangeTimes->{{3.591948720451333*^9, 3.5919487489181232`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"varEstim24", "=", RowBox[{"Table", "[", " ", RowBox[{ RowBox[{"{", RowBox[{"n", ",", RowBox[{"Variance", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"Mean", "[", RowBox[{ RowBox[{"Sort", "[", RowBox[{"RandomVariate", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"0", ",", "1"}], "]"}], ",", "n"}], "]"}], "]"}], "[", RowBox[{"[", RowBox[{ RowBox[{"Floor", "[", RowBox[{"0.38", " ", "n"}], "]"}], ";;", RowBox[{"Ceiling", "[", RowBox[{"0.62", " ", "n"}], "]"}]}], "]"}], "]"}], " ", "]"}], ",", RowBox[{"{", "1000", "}"}]}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"n", ",", "100", ",", "2000", ",", "100"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5919461559020243`*^9, 3.5919461891369476`*^9}, { 3.5919462208058276`*^9, 3.591946222235871*^9}, {3.5919487786769266`*^9, 3.5919488116978416`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"100", ",", "0.023953695189549047`"}], "}"}], ",", RowBox[{"{", RowBox[{"200", ",", "0.012298874146572271`"}], "}"}], ",", RowBox[{"{", RowBox[{"300", ",", "0.007649532297569435`"}], "}"}], ",", RowBox[{"{", RowBox[{"400", ",", "0.00617688953913628`"}], "}"}], ",", RowBox[{"{", RowBox[{"500", ",", "0.004460782814801153`"}], "}"}], ",", RowBox[{"{", RowBox[{"600", ",", "0.0039827927781873`"}], "}"}], ",", RowBox[{"{", RowBox[{"700", ",", "0.0033977732752177663`"}], "}"}], ",", RowBox[{"{", RowBox[{"800", ",", "0.002730827527523522`"}], "}"}], ",", RowBox[{"{", RowBox[{"900", ",", "0.00256529929635679`"}], "}"}], ",", RowBox[{"{", RowBox[{"1000", ",", "0.0022293897579021863`"}], "}"}], ",", RowBox[{"{", RowBox[{"1100", ",", "0.002091637705983606`"}], "}"}], ",", RowBox[{"{", RowBox[{"1200", ",", "0.0017530553447685625`"}], "}"}], ",", RowBox[{"{", RowBox[{"1300", ",", "0.0017237414088648404`"}], "}"}], ",", RowBox[{"{", RowBox[{"1400", ",", "0.0016446687686248432`"}], "}"}], ",", RowBox[{"{", RowBox[{"1500", ",", "0.0016004419615106371`"}], "}"}], ",", RowBox[{"{", RowBox[{"1600", ",", "0.001350021906176669`"}], "}"}], ",", RowBox[{"{", RowBox[{"1700", ",", "0.0012646207986980555`"}], "}"}], ",", RowBox[{"{", RowBox[{"1800", ",", "0.0013085219917949944`"}], "}"}], ",", RowBox[{"{", RowBox[{"1900", ",", "0.0012375543301905308`"}], "}"}], ",", RowBox[{"{", RowBox[{"2000", ",", "0.0011113518484843596`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.672138399293005*^9, 3.6730786448885064`*^9}] }, Open ]], Cell["fazemos o mesmo para a mediana, para comparar", "Text", CellChangeTimes->{{3.591948919356866*^9, 3.591948929755151*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"varMedianaNGrande", "=", RowBox[{"Table", "[", " ", RowBox[{ RowBox[{"{", RowBox[{"n", ",", RowBox[{"Variance", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"Median", "[", RowBox[{"RandomVariate", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"0", ",", "1"}], "]"}], ",", "n"}], "]"}], "]"}], ",", RowBox[{"{", "1000", "}"}]}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"n", ",", "100", ",", "2000", ",", "100"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.591948855457058*^9, 3.591948900566312*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"100", ",", "0.024786255184017322`"}], "}"}], ",", RowBox[{"{", RowBox[{"200", ",", "0.012272814678910896`"}], "}"}], ",", RowBox[{"{", RowBox[{"300", ",", "0.007484952953499866`"}], "}"}], ",", RowBox[{"{", RowBox[{"400", ",", "0.00638790743730034`"}], "}"}], ",", RowBox[{"{", RowBox[{"500", ",", "0.004868584923430223`"}], "}"}], ",", RowBox[{"{", RowBox[{"600", ",", "0.0041199526316446754`"}], "}"}], ",", RowBox[{"{", RowBox[{"700", ",", "0.0037075474418471924`"}], "}"}], ",", RowBox[{"{", RowBox[{"800", ",", "0.0029582634838269278`"}], "}"}], ",", RowBox[{"{", RowBox[{"900", ",", "0.0029653797698438824`"}], "}"}], ",", RowBox[{"{", RowBox[{"1000", ",", "0.002449662747192123`"}], "}"}], ",", RowBox[{"{", RowBox[{"1100", ",", "0.0022448812577414855`"}], "}"}], ",", RowBox[{"{", RowBox[{"1200", ",", "0.00195858089185349`"}], "}"}], ",", RowBox[{"{", RowBox[{"1300", ",", "0.0017846238236782447`"}], "}"}], ",", RowBox[{"{", RowBox[{"1400", ",", "0.001629408046591152`"}], "}"}], ",", RowBox[{"{", RowBox[{"1500", ",", "0.0015888484753810318`"}], "}"}], ",", RowBox[{"{", RowBox[{"1600", ",", "0.0015611300774280412`"}], "}"}], ",", RowBox[{"{", RowBox[{"1700", ",", "0.0013722181902746063`"}], "}"}], ",", RowBox[{"{", RowBox[{"1800", ",", "0.0013221159998322563`"}], "}"}], ",", RowBox[{"{", RowBox[{"1900", ",", "0.0012754781929020015`"}], "}"}], ",", RowBox[{"{", RowBox[{"2000", ",", "0.001119400235956071`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.6721384014671974`*^9, 3.6730786475820837`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"varMedianaNGrande", ",", "varEstim24"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"BaseStyle", "\[Rule]", RowBox[{"{", RowBox[{"FontSize", "\[Rule]", "14"}], "}"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.591946208499489*^9, 3.591946213541652*^9}, { 3.59194638741446*^9, 3.5919463895105176`*^9}, {3.591946433165729*^9, 3.5919464362018147`*^9}, {3.5919465049957275`*^9, 3.5919465576531906`*^9}, {3.591948984600648*^9, 3.5919491459251447`*^9}, 3.6730798686127777`*^9}], Cell[BoxData[ GraphicsBox[{{}, {{}, {RGBColor[0.368417, 0.506779, 0.709798], PointSize[0.012833333333333334`], AbsoluteThickness[1.6], StyleBox[PointBox[CompressedData[" 1:eJxTTMoPSmViYGAQAWIQDQGRDnPS5udOTJxpD+FnOiixZvjWKnVC+AeKHDau jj7UtaIOKl/pwFHffsNEqwrCd6h3+MGadEHkYzFUfZMDnzBbrOiDAgj/QavD 809RCV7xeVD9nQ6pR/3nvrXOgPAVehyuPjRTueUF5Tv0O6yMDl0TIJgC4RtM dFCcv7vodnIS1PxJDlt4BGUKuRMg/IApDpnHQ20dbWOh9k11mDUzfaL1pigI v2C6wwmNWH5ZDiifYaaD8XS1FbcmR0L4E2Y5vI9jFmioDoO6Z45DdcpRu+iV oRD+hrkO05tYPPY9DYG6b77DlYwnN6XCg+wBwflbDw== "]], FontSize->14]}, {RGBColor[0.880722, 0.611041, 0.142051], PointSize[0.012833333333333334`], AbsoluteThickness[1.6], StyleBox[PointBox[CompressedData[" 1:eJxTTMoPSmViYGAQAWIQDQGRDo8Kr7gGts+wh/AzHQ6uffpRzaATwj9Q5OAx 5/AG2dB6qHylwz2vPxoffCohfId6B/3td10qXIug6pscrKxWnXYKKIDwH7Q6 BB8QWbz0ajZUf6cD34W6u1/j0iB8hR6HCptFky4xp0LN63c4xm/aHOecBOEb THTYFB1VWK2UCDV/koOF39vrN3bEQPgBUxw0bk+xe2oL5T+Y6vApas8S209R EH7BdIeleydU7LWA8hlmOjD/ffo7TS4Mwp8wy0E7zNLOZEcI1D1zHFQrvE+F FYRC+BvmOrgpM99Y6waVd5jvsCnu79cI0yB7AAr7Ydc= "]], FontSize->14]}, {}}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{ FormBox["\"n\"", TraditionalForm], FormBox["\"variance\"", TraditionalForm]}, AxesOrigin->{70.3125, 0}, BaseStyle->{FontSize -> 14}, 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]], ImagePadding->All, Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& )}}, PlotRange->{{100., 2000.}, {0, 0.024786255184017322`}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.02], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.672138401504965*^9, 3.673078647654153*^9, 3.6730798702129707`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"comparacao", "=", RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"varMedianaNGrande", ",", "varEstim24"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0.004"}], "}"}]}], ",", RowBox[{"PlotLegends", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"BaseStyle", "\[Rule]", RowBox[{"{", RowBox[{"FontSize", "\[Rule]", "14"}], "}"}]}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5919491691237764`*^9, 3.5919491724958706`*^9}, 3.6730798790304384`*^9}], Cell[BoxData[ TemplateBox[{GraphicsBox[{{}, {{{}, { Hue[0.67, 0.6, 0.6], Directive[ PointSize[0.012833333333333334`], RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], FontSize -> 14], PointBox[CompressedData[" 1:eJxTTMoPSmViYGAQAWIQDQGRDnPS5udOTJxpD+FnOiixZvjWKnVC+AeKHDau jj7UtaIOKl/pwFHffsNEqwrCd6h3+MGadEHkYzFUfZMDnzBbrOiDAgj/QavD 809RCV7xeVD9nQ6pR/3nvrXOgPAVehyuPjRTueUF5Tv0O6yMDl0TIJgC4RtM dFCcv7vodnIS1PxJDlt4BGUKuRMg/IApDpnHQ20dbWOh9k11mDUzfaL1pigI v2C6wwmNWH5ZDiifYaaD8XS1FbcmR0L4E2Y5vI9jFmioDoO6Z45DdcpRu+iV oRD+hrkO05tYPPY9DYG6b77DlYwnN6XCg+wBwflbDw== "]]}, { Hue[0.9060679774997897, 0.6, 0.6], Directive[ PointSize[0.012833333333333334`], RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], FontSize -> 14], PointBox[CompressedData[" 1:eJxTTMoPSmViYGAQAWIQDQGRDo8Kr7gGts+wh/AzHQ6uffpRzaATwj9Q5OAx 5/AG2dB6qHylwz2vPxoffCohfId6B/3td10qXIug6pscrKxWnXYKKIDwH7Q6 BB8QWbz0ajZUf6cD34W6u1/j0iB8hR6HCptFky4xp0LN63c4xm/aHOecBOEb THTYFB1VWK2UCDV/koOF39vrN3bEQPgBUxw0bk+xe2oL5T+Y6vApas8S209R EH7BdIeleydU7LWA8hlmOjD/ffo7TS4Mwp8wy0E7zNLOZEcI1D1zHFQrvE+F FYRC+BvmOrgpM99Y6waVd5jvsCnu79cI0yB7AAr7Ydc= "]]}, {}}}, {}}, {DisplayFunction -> Identity, PlotRangePadding -> {{ Scaled[0.02], Scaled[0.02]}, {0, 0}}, AxesOrigin -> {0, 0}, PlotRange -> {{0, 2000.}, {0, 0.004}}, PlotRangeClipping -> True, ImagePadding -> All, DisplayFunction -> Identity, AspectRatio -> NCache[GoldenRatio^(-1), 0.6180339887498948], Axes -> {True, True}, AxesLabel -> { FormBox["\"n\"", TraditionalForm], FormBox["\"variance\"", TraditionalForm]}, AxesOrigin -> {0, 0}, BaseStyle -> {FontSize -> 14}, 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 -> {"CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& )}}, PlotRange -> {{0, 2000.}, {0, 0.004}}, PlotRangeClipping -> True, PlotRangePadding -> {{ Scaled[0.02], Scaled[0.02]}, {0, 0}}, Ticks -> {Automatic, Automatic}}],FormBox[ FormBox[ TemplateBox[{"\"mediana\"", "\"media central\""}, "PointLegend", DisplayFunction -> (FormBox[ StyleBox[ StyleBox[ PaneBox[ TagBox[ GridBox[{{ TagBox[ GridBox[{{ GraphicsBox[{{}, { Directive[ EdgeForm[ Directive[ Opacity[0.3], GrayLevel[0]]], PointSize[0.462], RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], FontSize -> 14], { Directive[ EdgeForm[ Directive[ Opacity[0.3], GrayLevel[0]]], PointSize[0.462], RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], FontSize -> 14], PointBox[ NCache[{ Scaled[{ Rational[1, 2], Rational[1, 2]}]}, { Scaled[{0.5, 0.5}]}]]}}}, AspectRatio -> Full, ImageSize -> {10, 10}, PlotRangePadding -> None, ImagePadding -> Automatic, BaselinePosition -> (Scaled[0.1] -> Baseline)], #}, { GraphicsBox[{{}, { Directive[ EdgeForm[ Directive[ Opacity[0.3], GrayLevel[0]]], PointSize[0.462], RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], FontSize -> 14], { Directive[ EdgeForm[ Directive[ Opacity[0.3], GrayLevel[0]]], PointSize[0.462], RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], FontSize -> 14], PointBox[ NCache[{ Scaled[{ Rational[1, 2], Rational[1, 2]}]}, { Scaled[{0.5, 0.5}]}]]}}}, AspectRatio -> Full, ImageSize -> {10, 10}, PlotRangePadding -> None, ImagePadding -> Automatic, BaselinePosition -> (Scaled[0.1] -> Baseline)], #2}}, GridBoxAlignment -> { "Columns" -> {Center, Left}, "Rows" -> {{Baseline}}}, AutoDelete -> False, GridBoxDividers -> { "Columns" -> {{False}}, "Rows" -> {{False}}}, GridBoxItemSize -> {"Columns" -> {{All}}, "Rows" -> {{All}}}, GridBoxSpacings -> { "Columns" -> {{0.5}}, "Rows" -> {{0.8}}}], "Grid"]}}, GridBoxAlignment -> {"Columns" -> {{Left}}, "Rows" -> {{Top}}}, AutoDelete -> False, GridBoxItemSize -> { "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, GridBoxSpacings -> {"Columns" -> {{1}}, "Rows" -> {{0}}}], "Grid"], Alignment -> Left, AppearanceElements -> None, ImageMargins -> {{5, 5}, {5, 5}}, ImageSizeAction -> "ResizeToFit"], LineIndent -> 0, StripOnInput -> False], { FontFamily -> "Arial"}, Background -> Automatic, StripOnInput -> False], TraditionalForm]& ), InterpretationFunction :> (RowBox[{"PointLegend", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Directive", "[", RowBox[{ RowBox[{"PointSize", "[", "0.012833333333333334`", "]"}], ",", InterpretationBox[ ButtonBox[ TooltipBox[ GraphicsBox[{{ GrayLevel[0], RectangleBox[{0, 0}]}, { GrayLevel[0], RectangleBox[{1, -1}]}, { RGBColor[0.368417, 0.506779, 0.709798], RectangleBox[{0, -1}, {2, 1}]}}, AspectRatio -> 1, Frame -> True, FrameStyle -> RGBColor[ 0.24561133333333335`, 0.3378526666666667, 0.4731986666666667], FrameTicks -> None, PlotRangePadding -> None, ImageSize -> Dynamic[{ Automatic, 1.35 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}]], "RGBColor[0.368417, 0.506779, 0.709798]"], Appearance -> None, BaseStyle -> {}, BaselinePosition -> Baseline, DefaultBaseStyle -> {}, ButtonFunction :> With[{Typeset`box$ = EvaluationBox[]}, If[ Not[ AbsoluteCurrentValue["Deployed"]], SelectionMove[Typeset`box$, All, Expression]; FrontEnd`Private`$ColorSelectorInitialAlpha = 1; FrontEnd`Private`$ColorSelectorInitialColor = RGBColor[0.368417, 0.506779, 0.709798]; FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; MathLink`CallFrontEnd[ FrontEnd`AttachCell[Typeset`box$, FrontEndResource["RGBColorValueSelector"], { 0, {Left, Bottom}}, {Left, Top}, "ClosingActions" -> { "SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> Automatic, Method -> "Preemptive"], RGBColor[0.368417, 0.506779, 0.709798], Editable -> False, Selectable -> False], ",", RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}], ",", RowBox[{"FontSize", "\[Rule]", "14"}]}], "]"}], ",", RowBox[{"Directive", "[", RowBox[{ RowBox[{"PointSize", "[", "0.012833333333333334`", "]"}], ",", InterpretationBox[ ButtonBox[ TooltipBox[ GraphicsBox[{{ GrayLevel[0], RectangleBox[{0, 0}]}, { GrayLevel[0], RectangleBox[{1, -1}]}, { RGBColor[0.880722, 0.611041, 0.142051], RectangleBox[{0, -1}, {2, 1}]}}, AspectRatio -> 1, Frame -> True, FrameStyle -> RGBColor[ 0.587148, 0.40736066666666665`, 0.09470066666666668], FrameTicks -> None, PlotRangePadding -> None, ImageSize -> Dynamic[{ Automatic, 1.35 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}]], "RGBColor[0.880722, 0.611041, 0.142051]"], Appearance -> None, BaseStyle -> {}, BaselinePosition -> Baseline, DefaultBaseStyle -> {}, ButtonFunction :> With[{Typeset`box$ = EvaluationBox[]}, If[ Not[ AbsoluteCurrentValue["Deployed"]], SelectionMove[Typeset`box$, All, Expression]; FrontEnd`Private`$ColorSelectorInitialAlpha = 1; FrontEnd`Private`$ColorSelectorInitialColor = RGBColor[0.880722, 0.611041, 0.142051]; FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; MathLink`CallFrontEnd[ FrontEnd`AttachCell[Typeset`box$, FrontEndResource["RGBColorValueSelector"], { 0, {Left, Bottom}}, {Left, Top}, "ClosingActions" -> { "SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> Automatic, Method -> "Preemptive"], RGBColor[0.880722, 0.611041, 0.142051], Editable -> False, Selectable -> False], ",", RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}], ",", RowBox[{"FontSize", "\[Rule]", "14"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{#, ",", #2}], "}"}], ",", RowBox[{"LegendMarkers", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"False", ",", "Automatic"}], "}"}], ",", RowBox[{"{", RowBox[{"False", ",", "Automatic"}], "}"}]}], "}"}]}], ",", RowBox[{"Joined", "\[Rule]", RowBox[{"{", RowBox[{"False", ",", "False"}], "}"}]}], ",", RowBox[{"LabelStyle", "\[Rule]", RowBox[{"{", "}"}]}], ",", RowBox[{"LegendLayout", "\[Rule]", "\"Column\""}]}], "]"}]& ), Editable -> True], TraditionalForm], TraditionalForm]}, "Legended", DisplayFunction->(GridBox[{{ TagBox[ ItemBox[ PaneBox[ TagBox[#, "SkipImageSizeLevel"], Alignment -> {Center, Baseline}, BaselinePosition -> Baseline], DefaultBaseStyle -> "Labeled"], "SkipImageSizeLevel"], ItemBox[#2, DefaultBaseStyle -> "LabeledLabel"]}}, GridBoxAlignment -> {"Columns" -> {{Center}}, "Rows" -> {{Center}}}, AutoDelete -> False, GridBoxItemSize -> Automatic, BaselinePosition -> {1, 1}]& ), Editable->True, InterpretationFunction->(RowBox[{"Legended", "[", RowBox[{#, ",", RowBox[{"Placed", "[", RowBox[{#2, ",", "After"}], "]"}]}], "]"}]& )]], "Output", CellChangeTimes->{3.6721384015518413`*^9, 3.673078647774827*^9, 3.673079879847658*^9}] }, Open ]], Cell["\<\ Para que os pontos simulados variem mais suavemente, basta repetir os \ calculos com tiragens de 10000 simula\[CCedilla]\[OTilde]es...\ \>", "Text", CellChangeTimes->{{3.591949227877437*^9, 3.591949257908275*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"{", RowBox[{"comparacao", ",", RowBox[{"Plot", "[", RowBox[{ RowBox[{"2", "/", "n"}], ",", RowBox[{"{", RowBox[{"n", ",", "200", ",", "2000"}], "}"}]}], "]"}]}], "}"}], "]"}]], "Input", CellChangeTimes->{{3.59194917671799*^9, 3.5919492073628397`*^9}}], Cell[BoxData[ TemplateBox[{GraphicsBox[{{{}, {{{}, { Hue[0.67, 0.6, 0.6], Directive[ PointSize[0.012833333333333334`], RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], FontSize -> 14], PointBox[CompressedData[" 1:eJxTTMoPSmViYGAQAWIQDQGRDnPS5udOTJxpD+FnOiixZvjWKnVC+AeKHDau jj7UtaIOKl/pwFHffsNEqwrCd6h3+MGadEHkYzFUfZMDnzBbrOiDAgj/QavD 809RCV7xeVD9nQ6pR/3nvrXOgPAVehyuPjRTueUF5Tv0O6yMDl0TIJgC4RtM dFCcv7vodnIS1PxJDlt4BGUKuRMg/IApDpnHQ20dbWOh9k11mDUzfaL1pigI v2C6wwmNWH5ZDiifYaaD8XS1FbcmR0L4E2Y5vI9jFmioDoO6Z45DdcpRu+iV oRD+hrkO05tYPPY9DYG6b77DlYwnN6XCg+wBwflbDw== "]]}, { Hue[0.9060679774997897, 0.6, 0.6], Directive[ PointSize[0.012833333333333334`], RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], FontSize -> 14], PointBox[CompressedData[" 1:eJxTTMoPSmViYGAQAWIQDQGRDo8Kr7gGts+wh/AzHQ6uffpRzaATwj9Q5OAx 5/AG2dB6qHylwz2vPxoffCohfId6B/3td10qXIug6pscrKxWnXYKKIDwH7Q6 BB8QWbz0ajZUf6cD34W6u1/j0iB8hR6HCptFky4xp0LN63c4xm/aHOecBOEb THTYFB1VWK2UCDV/koOF39vrN3bEQPgBUxw0bk+xe2oL5T+Y6vApas8S209R EH7BdIeleydU7LWA8hlmOjD/ffo7TS4Mwp8wy0E7zNLOZEcI1D1zHFQrvE+F FYRC+BvmOrgpM99Y6waVd5jvsCnu79cI0yB7AAr7Ydc= "]]}, {}}}, {}}, {{{{}, {}, { Directive[ Opacity[1.], RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6]], LineBox[CompressedData[" 1:eJwVxXk8VAkcAPDaLklsSSlilaKGdCgi/X4l1tHxkBzlCGHMJOPN85ppiMrI UY5xN+ZFJR3kzFYqDV27XUTa7sSglatUyrG7f3w/Xz2/fc57fhk3bpz3f/5/ UyNH7Xg5D3WmYF1pjwjSzofo25mH45TRFm7XXBF4fQxYErSDxIeekWq6QiE4 q062ZPmQaFZholgYKgS7lWc39waRuNa9MXmVnxBWCrv20QISrwby9EgHISgp h16Ok5JYq5V0nqclhMoltHVRK4mVjzZO870pABV2vFdXKB/VPiuVnlYVwNWO 4lROHIWm8lL1pfdoEGh0PPuUQmGf+HNF8XUazK1/0w7LpVB7BZFsUUHDZVla Ib+YwvXbau/sk9FQ7kpfi3xKIdUy3GlL0XC+DtuTdSOQ2HRUnqdPQy7z1Lyy OgJ9TPOnnoyPgANuQ2+HO2j8oavfKg+hwO6AtOxVH42HLBTHG/0omMXA4Zoh GnNNF2l3elJQoohdLJq2H1m5M2/oO1LQSs3c+3PZfpz2VzLnA4sCh3TWjyFq Pz558vfHsB4+aDZ4zfo2QYA6Bu+sCmg+VDnI7fp1hag6w3iNQQ4Jm/ePbEle IkRrG37nQCoJrWfMXYxXCXHZcHVbXQIJquMvebFthXjv9dApSkRCcPWJ8Pdc IWq4nZLO8SVBW5+SNlwRojAuIZEwJOHIiEFfmcsBJFg/CsNqwmF76bGs8DgR XmyTKFwGeHBFUrzp6ocovG0w2r1eEgaOkYvE5YnR6O9W9LGhPxS62YO7tl6P QfvY37uqX3HBXnPirRMvDqHlULjadwMOxE+y2r554DBqOgXTJhfYMMyiPev0 YrF0NONxqkUwsGb1ffpnnRh3iqya1SYHQdLqkC8v9sXh4wkHHBOe74GGDNWl xnlHcexRm47WwwBorO7vXn45Hlt3R+bJu/xhkvqzBT8VCcg/Y8M1UPKHLyZD 7RETkjAgarTwva0fzBbP175idgxflWyImi7eDVpa3BjvwONo2HMsI7vFF17x hsiq6GRst1/OK0Rf6BxXz2V7pODX5iPaY4QP1J5bqiA6U7BqBoe31c4blMMX f2PzU9E6LbIcXL3g07OCDYEqabj456WIdd67oER6w9kyNw2z2J7sgdidoDTm mfd5lQTpk5No5pwnNEUnHXz9pwTHnxb4u7zxgCh/L8OTbukYs6hJ4sLygFut Jcq2fekoPvOmt5l0h1Y5U1EflYG5b5VSaxrcgJIv+sJXzsSxuRta71i4Qb63 UahSQSY6TOZRpswOuLCzsYRvkoVmc91nm2vsgG2D3w3r72Vh+tkizXdiV2jo vcHudsvGEMFNPTslV6ifOm+jbVs2OpaQKmWp28GpdtDlsCgHj8jPuz6dvR1q WsKMcpRzkZ9B7DE75wK1idz1aadzkWvtU5Vl5gJrfEY6aixP4Gi3RL37gTMk SIzda++dQF9FwUj0XmegA96MP7tbijYJKvMtlJ1hY+ay8MA+KapfcLxtVeYE BzvUhW6iPLTyiMq38XIC12bBNdYsGRrKh1nMJCeYqn+w7nm+DKfLNRRtlwlI L3DwdTotQwuBvDGxigAdvdkj9wtlWKKj/m1lJQGmOsVrrl2Q4dPgEEVMGQG+ mi/P5VXJULxT4L3gIgF/qJil+N2X4XAHu2hvPgGBgz27uvtkOP77qK1xEgH9 5NWhgM8yvNLAUX2ZQIBoIDbz9aAM032eTYuPJyCtV/vJox8yzOoKK2oXE3Cj y8G6bCKDpS+cg/JjCNB4U7gkQpNB11spfQtoAk7uCr/TM4/BmV8aZzRRBCx9 aeUfNJ9BYsQmOJZPAD5vkrovYLC3K7Gyk0cAt/GXXy2NGLTkrL5fxiXgG/Ho YsUyBk0oaVoAh4BDj3PsjVYwaMR5nDInhIDsB8sPz1/D4Pv+QFZUEAELHYd1 M80ZbKbL764IJKDk/t0aVUsG7ZuystsDCFhrJ/GIs2LwVqKyLMefgPo73l/H gMHyxaMtW/wI+Bflt8L7 "]]}}}, {}}}, {DisplayFunction -> Identity, PlotRangePadding -> {{ Scaled[0.02], Scaled[0.02]}, {0, 0}}, AxesOrigin -> {0, 0}, PlotRange -> {{0, 2000.}, {0, 0.004}}, PlotRangeClipping -> True, ImagePadding -> All, DisplayFunction -> Identity, AspectRatio -> NCache[GoldenRatio^(-1), 0.6180339887498948], Axes -> {True, True}, AxesLabel -> { FormBox["\"n\"", TraditionalForm], FormBox["\"variance\"", TraditionalForm]}, AxesOrigin -> {0, 0}, BaseStyle -> {FontSize -> 14}, 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 -> {"CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& )}}, PlotRange -> {{0, 2000.}, {0, 0.004}}, PlotRangeClipping -> True, PlotRangePadding -> {{ Scaled[0.02], Scaled[0.02]}, {0, 0}}, Ticks -> {Automatic, Automatic}}],FormBox[ FormBox[ TemplateBox[{"\"mediana\"", "\"media central\""}, "PointLegend", DisplayFunction -> (FormBox[ StyleBox[ StyleBox[ PaneBox[ TagBox[ GridBox[{{ TagBox[ GridBox[{{ GraphicsBox[{{}, { Directive[ EdgeForm[ Directive[ Opacity[0.3], GrayLevel[0]]], PointSize[0.462], RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], FontSize -> 14], { Directive[ EdgeForm[ Directive[ Opacity[0.3], GrayLevel[0]]], PointSize[0.462], RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], FontSize -> 14], PointBox[ NCache[{ Scaled[{ Rational[1, 2], Rational[1, 2]}]}, { Scaled[{0.5, 0.5}]}]]}}}, AspectRatio -> Full, ImageSize -> {10, 10}, PlotRangePadding -> None, ImagePadding -> Automatic, BaselinePosition -> (Scaled[0.1] -> Baseline)], #}, { GraphicsBox[{{}, { Directive[ EdgeForm[ Directive[ Opacity[0.3], GrayLevel[0]]], PointSize[0.462], RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], FontSize -> 14], { Directive[ EdgeForm[ Directive[ Opacity[0.3], GrayLevel[0]]], PointSize[0.462], RGBColor[0.880722, 0.611041, 0.142051], AbsoluteThickness[1.6], FontSize -> 14], PointBox[ NCache[{ Scaled[{ Rational[1, 2], Rational[1, 2]}]}, { Scaled[{0.5, 0.5}]}]]}}}, AspectRatio -> Full, ImageSize -> {10, 10}, PlotRangePadding -> None, ImagePadding -> Automatic, BaselinePosition -> (Scaled[0.1] -> Baseline)], #2}}, GridBoxAlignment -> { "Columns" -> {Center, Left}, "Rows" -> {{Baseline}}}, AutoDelete -> False, GridBoxDividers -> { "Columns" -> {{False}}, "Rows" -> {{False}}}, GridBoxItemSize -> {"Columns" -> {{All}}, "Rows" -> {{All}}}, GridBoxSpacings -> { "Columns" -> {{0.5}}, "Rows" -> {{0.8}}}], "Grid"]}}, GridBoxAlignment -> {"Columns" -> {{Left}}, "Rows" -> {{Top}}}, AutoDelete -> False, GridBoxItemSize -> { "Columns" -> {{Automatic}}, "Rows" -> {{Automatic}}}, GridBoxSpacings -> {"Columns" -> {{1}}, "Rows" -> {{0}}}], "Grid"], Alignment -> Left, AppearanceElements -> None, ImageMargins -> {{5, 5}, {5, 5}}, ImageSizeAction -> "ResizeToFit"], LineIndent -> 0, StripOnInput -> False], { FontFamily -> "Arial"}, Background -> Automatic, StripOnInput -> False], TraditionalForm]& ), InterpretationFunction :> (RowBox[{"PointLegend", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Directive", "[", RowBox[{ RowBox[{"PointSize", "[", "0.012833333333333334`", "]"}], ",", InterpretationBox[ ButtonBox[ TooltipBox[ GraphicsBox[{{ GrayLevel[0], RectangleBox[{0, 0}]}, { GrayLevel[0], RectangleBox[{1, -1}]}, { RGBColor[0.368417, 0.506779, 0.709798], RectangleBox[{0, -1}, {2, 1}]}}, AspectRatio -> 1, Frame -> True, FrameStyle -> RGBColor[ 0.24561133333333335`, 0.3378526666666667, 0.4731986666666667], FrameTicks -> None, PlotRangePadding -> None, ImageSize -> Dynamic[{ Automatic, 1.35 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}]], "RGBColor[0.368417, 0.506779, 0.709798]"], Appearance -> None, BaseStyle -> {}, BaselinePosition -> Baseline, DefaultBaseStyle -> {}, ButtonFunction :> With[{Typeset`box$ = EvaluationBox[]}, If[ Not[ AbsoluteCurrentValue["Deployed"]], SelectionMove[Typeset`box$, All, Expression]; FrontEnd`Private`$ColorSelectorInitialAlpha = 1; FrontEnd`Private`$ColorSelectorInitialColor = RGBColor[0.368417, 0.506779, 0.709798]; FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; MathLink`CallFrontEnd[ FrontEnd`AttachCell[Typeset`box$, FrontEndResource["RGBColorValueSelector"], { 0, {Left, Bottom}}, {Left, Top}, "ClosingActions" -> { "SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> Automatic, Method -> "Preemptive"], RGBColor[0.368417, 0.506779, 0.709798], Editable -> False, Selectable -> False], ",", RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}], ",", RowBox[{"FontSize", "\[Rule]", "14"}]}], "]"}], ",", RowBox[{"Directive", "[", RowBox[{ RowBox[{"PointSize", "[", "0.012833333333333334`", "]"}], ",", InterpretationBox[ ButtonBox[ TooltipBox[ GraphicsBox[{{ GrayLevel[0], RectangleBox[{0, 0}]}, { GrayLevel[0], RectangleBox[{1, -1}]}, { RGBColor[0.880722, 0.611041, 0.142051], RectangleBox[{0, -1}, {2, 1}]}}, AspectRatio -> 1, Frame -> True, FrameStyle -> RGBColor[ 0.587148, 0.40736066666666665`, 0.09470066666666668], FrameTicks -> None, PlotRangePadding -> None, ImageSize -> Dynamic[{ Automatic, 1.35 CurrentValue["FontCapHeight"]/ AbsoluteCurrentValue[Magnification]}]], "RGBColor[0.880722, 0.611041, 0.142051]"], Appearance -> None, BaseStyle -> {}, BaselinePosition -> Baseline, DefaultBaseStyle -> {}, ButtonFunction :> With[{Typeset`box$ = EvaluationBox[]}, If[ Not[ AbsoluteCurrentValue["Deployed"]], SelectionMove[Typeset`box$, All, Expression]; FrontEnd`Private`$ColorSelectorInitialAlpha = 1; FrontEnd`Private`$ColorSelectorInitialColor = RGBColor[0.880722, 0.611041, 0.142051]; FrontEnd`Private`$ColorSelectorUseMakeBoxes = True; MathLink`CallFrontEnd[ FrontEnd`AttachCell[Typeset`box$, FrontEndResource["RGBColorValueSelector"], { 0, {Left, Bottom}}, {Left, Top}, "ClosingActions" -> { "SelectionDeparture", "ParentChanged", "EvaluatorQuit"}]]]], BaseStyle -> Inherited, Evaluator -> Automatic, Method -> "Preemptive"], RGBColor[0.880722, 0.611041, 0.142051], Editable -> False, Selectable -> False], ",", RowBox[{"AbsoluteThickness", "[", "1.6`", "]"}], ",", RowBox[{"FontSize", "\[Rule]", "14"}]}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{#, ",", #2}], "}"}], ",", RowBox[{"LegendMarkers", "\[Rule]", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"False", ",", "Automatic"}], "}"}], ",", RowBox[{"{", RowBox[{"False", ",", "Automatic"}], "}"}]}], "}"}]}], ",", RowBox[{"Joined", "\[Rule]", RowBox[{"{", RowBox[{"False", ",", "False"}], "}"}]}], ",", RowBox[{"LabelStyle", "\[Rule]", RowBox[{"{", "}"}]}], ",", RowBox[{"LegendLayout", "\[Rule]", "\"Column\""}]}], "]"}]& ), Editable -> True], TraditionalForm], TraditionalForm]}, "Legended", DisplayFunction->(GridBox[{{ TagBox[ ItemBox[ PaneBox[ TagBox[#, "SkipImageSizeLevel"], Alignment -> {Center, Baseline}, BaselinePosition -> Baseline], DefaultBaseStyle -> "Labeled"], "SkipImageSizeLevel"], ItemBox[#2, DefaultBaseStyle -> "LabeledLabel"]}}, GridBoxAlignment -> {"Columns" -> {{Center}}, "Rows" -> {{Center}}}, AutoDelete -> False, GridBoxItemSize -> Automatic, BaselinePosition -> {1, 1}]& ), Editable->True, InterpretationFunction->(RowBox[{"Legended", "[", RowBox[{#, ",", RowBox[{"Placed", "[", RowBox[{#2, ",", "After"}], "]"}]}], "]"}]& )]], "Output", CellChangeTimes->{3.672138401605234*^9, 3.6730786478125896`*^9, 3.673079889604308*^9}] }, Open ]] }, Open ]] }, Open ]] }, Open ]], Cell["Maximum Likelihood estimate of Cauchy localization parameter", "Section", CellChangeTimes->{{3.673078490013589*^9, 3.6730785062114573`*^9}}], Cell[CellGroupData[{ Cell["Cauchy distribution location parameter with unity width", "Section", CellChangeTimes->{{3.6707447525522213`*^9, 3.670744775623436*^9}}], Cell[CellGroupData[{ Cell["preparing / testing", "Subsection", CellChangeTimes->{{3.670745015622259*^9, 3.6707450367204742`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"PDF", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"x0", ",", "\[Gamma]"}], "]"}], ",", "x"}], "]"}]], "Input", CellChangeTimes->{{3.670582737081394*^9, 3.6705827487964554`*^9}, { 3.670582792318639*^9, 3.670582808275345*^9}, 3.670582848798712*^9, { 3.6730785961168203`*^9, 3.6730786197378516`*^9}}], Cell[BoxData[ FractionBox["1", RowBox[{"\[Pi]", " ", RowBox[{"(", RowBox[{"1", "+", FractionBox[ SuperscriptBox[ RowBox[{"(", RowBox[{"x", "-", "x0"}], ")"}], "2"], SuperscriptBox["\[Gamma]", "2"]]}], ")"}], " ", "\[Gamma]"}]]], "Output", CellChangeTimes->{ 3.6705828090069375`*^9, 3.670582849645963*^9, 3.670586207378999*^9, 3.670605104084545*^9, 3.670744780690857*^9, 3.673078028492376*^9, { 3.6730785779202995`*^9, 3.673078599770576*^9}, 3.6730786478282156`*^9}] }, Open ]], Cell[BoxData[ RowBox[{ RowBox[{ RowBox[{"logLCauchy", "[", "data_", "]"}], "[", "a_", "]"}], ":=", RowBox[{"-", RowBox[{"Total", "[", RowBox[{"Log", "[", RowBox[{"1", "+", RowBox[{ RowBox[{"(", RowBox[{"data", "-", "a"}], ")"}], "^", "2"}]}], "]"}], " ", "]"}]}]}]], "Input"], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"w", "=", RowBox[{"RandomVariate", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"10", ",", "1"}], "]"}], ",", "10"}], "]"}]}]], "Input", CellChangeTimes->{{3.670583047493455*^9, 3.6705830780980043`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ "10.255374883302288`", ",", "9.69270170358252`", ",", "10.518751633969817`", ",", "9.218126956325765`", ",", "5.694841316293542`", ",", "9.023215965786695`", ",", "11.351704175869152`", ",", "6.250467582926023`", ",", "4.625233593574584`", ",", "7.82236825034389`"}], "}"}]], "Output", CellChangeTimes->{3.6705830787361717`*^9, 3.6705835522604785`*^9, 3.670586207394625*^9, 3.6706051041001716`*^9, 3.670744784981822*^9, 3.6707449230130606`*^9, 3.6730780285080013`*^9, 3.6730785780205684`*^9, 3.6730786478594675`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"logLc", "=", RowBox[{"logLCauchy", "[", "w", "]"}]}]], "Input", CellChangeTimes->{{3.670583827524956*^9, 3.6705838761190896`*^9}}], Cell[BoxData[ RowBox[{"logLCauchy", "[", RowBox[{"{", RowBox[{ "10.255374883302288`", ",", "9.69270170358252`", ",", "10.518751633969817`", ",", "9.218126956325765`", ",", "5.694841316293542`", ",", "9.023215965786695`", ",", "11.351704175869152`", ",", "6.250467582926023`", ",", "4.625233593574584`", ",", "7.82236825034389`"}], "}"}], "]"}]], "Output", CellChangeTimes->{{3.67058385340241*^9, 3.6705838766812882`*^9}, 3.670586207432392*^9, 3.6706051041314244`*^9, 3.670744855399873*^9, 3.67074492501453*^9, 3.6730780285080013`*^9, 3.673078578036195*^9, 3.6730786478750935`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Plot", "[", RowBox[{ RowBox[{"logLc", "[", "a", "]"}], ",", RowBox[{"{", RowBox[{"a", ",", "5", ",", "15"}], "}"}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"BaseStyle", "\[Rule]", RowBox[{"{", RowBox[{"FontSize", "\[Rule]", "14"}], "}"}]}]}], "]"}]], "Input", CellChangeTimes->{{3.670583855288988*^9, 3.670583867012478*^9}, { 3.6730798981680794`*^9, 3.6730799112838926`*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {}, {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[ 1.], LineBox[CompressedData[" 1:eJwV1nk4VWsXAHBucTUowyFEnH2Oss8gswax7kmRkiJDMs+RJErTV8kYynGT KFFUF6WMJcldXUoqQpIilbGIUgiJ7+2P/ezn9+y913qHtd5nMz2CrL3/EBER +USu3/cIhSJJEREGbH3isWaToiEmOPfNZc9iwPW6ZNMJOUNMzlKWMBdnQFmV yadrMoZ4hRfzh3A+A7aF/pgze74hVoPTD1VFBhxjymY9nDZA8Z3iH0x0GND/ cf+Ud7cBxt51LDnmxQDjgFvrlQoMMMFx1o6Zx+R9v9RXWzYYYNK/ac96s+Tg 9uZjNntC9NFqPDHJebc8pP3RpVN5Vw9jQnXOt+gvgo2aTxQGGXp4T1apxnCB AritDfj74AFdzIl2N9NoVYAL470D8FMHk5SHBkSKFaFv88dZktE6ODK6tLnt pBLUdfvXt6vooFT2vwvXWy+GjD3Sz9UrtVEj5GJQKlcZhnkRh8u3aOOEfBt1 aEwZNB1fNyh2aqHpdLWiy30V6Lpx8aD5bi2MlxjR45xdAl621p/MWpdjfO+J 7CgXVeh3OVN6uEcTq0K4F7yU1eBw+uurcV18/HU2IDomRg0OLGdUPennISNr x+K9A2rgYpR3p0eMh7kmu/9rsGdC8PCdAFzExUcbCnd4b2eC6QWhRRCDi902 IyMTjkxI3lhXrirNRTW/Q8tYLkx41vXKM2ouF1OF0Qn7vJjAqd7q5/WLg9Ef MuyUgplw02Sifl8XBz0invd7xjGhaNjz/rxCDh4/JRM5Hs8E2y5JqR/5HMw4 Z6ty6hQTKi0GLnXncfD19Tar20ImlLXXnXl0hYNbmvtKJM4xQf7nQadrqRw0 Vhc5lp/NhD3jW/a8Oc5BpRot2R8VTBD9ZRpsZMPBozsqjuRXMmHBRNh1ky0c 7Pxi1uuJTEg3sXdeb8nBXAXXuw1VTGh/n2DqYsbBFf6nXPOeMGH0nFhf2WoO 2kn2X3duZUKovutFVXUO/m2TbVr9nQku4JDQOUnjWJ/mzUOjTFATyCl4jNPo eKR8kfYPJhjWqX7vGaWRebWx/+IkEx5Nvb41+ZXGgrEZYZgoBS/m3xpy6KOx Ps3pLWchBV3K+sWDzTTOeS+3P4lDwXO7ao/MYhpHnml0m/AoqI2bbq8opLHj 7mrrIT4F2+o1lr27RWPJGXfNjdoU3K7LN9K5QaOLWX7v7BUUDKVzu2Sukuc3 Te0PrKNg88DktQfnaMy8YP9wqRkFi1MPSGqm0BgX66/70pyC9ECTiMxkGl09 hAt1NlFgXH4xLymJxrny7TUD1hTUN2wwK44nz/8XssLVjYLZncHHh47TaOEf fU3Sg4LmsUvu+cdo1LNPY1R4UtB7PD8i+Cj5Xvvfrwq+FCyx7Emec4SMp3tu blMgBeKqCToHw8h4GlUWhQdRULTlP3eP/WQ8lVpRy4MpcJKN6dmyj+RLtXNP CKWg/H7a9zUhJN6mLMX1hylQjFL2dgiicXRFaczIEQp4Fk6hu3bT+F798WjW UQrkp1M+RgXSWDoz2DgTToFyk+DP2gAa3YpXxt2NpSDnbq6l0I+M/9Kmcd84 CkwvbS597kuj/ilXH/kECuhlRokM4rm+UYKQRAqEmSWZN71JPKXGSV4KBfte f3Zge9JoXhMXU3mOgp8b9SzSPGh8E2LK2JJGgeBjj6gs8fSzO/yQdAqM5sk9 ZrjTKDwYXD47g4LD70xb091opJZyzVIyKThweW0RTbz+WIbb3SwK3MLKH211 pbGV6zBocYWCm/pX/xh0odG/VfpQ+1UK3rqZzjtFPBX5VHz3PxQUyxp/0CZO 1I46M5NDQUGHVlSbM4nfYayWlEdB4dviHyeJS+LGb1A3KHhSZ7dyDXFrV8Cj dbcoSFNdxCtyIvGF6javCii417e8I4R4yujdO78iCizPu7qvIj79KXXXZDEF Ij4yt8WI1VKsJ+JLKZC6lN/+cgeNRYL50Sp3SPw3iS/ziE2/PJS5VUaB+pTU 5UjilgvHMqGcgoBtJ0w8iXear+Q13aNgrESzcD3x1Mi3Ms/7FHyJdhzTJD59 +ca60UoKlnauXahMrLbZpykayfyzqSnJ3/kmVV0V/iPzi9K4P/t3vn9eD+RW UXBIL9xa5Hc+mzMHVj+kQNfXuuK3/UQsxeoeUdA2em9SjHjyhvjfLo8p8P/6 SlKKOGE7LvlaS+KVlo8uIV4ifuh6+FMKTOaGFOkQFxTprpCtI/2xfK7pRmKB 62D1lXoKNk1E5PkRN8/7Z6tBAwWOfR09ccQ+ZW4dNY3EKxljhcTjXkoB219Q wPhToe0tcbx084/+Zgoe271OXkDWc0nlqcgjLRR06wpYpsQF/mbSC1opiHho GH2UWLBINCPzNekf7/P37xM3V5VztNsoOFKyrlaU7KfPntA7/7WT/SiVyrUg HlfWNN3WQYHV2zqnVGLl/Zedwz6QereyNhGQerlJ7eiX6CL1lZoTmEEMzxlh 57sp0NNfGDRD7KURK7zfR0HT2V+fnpN6HGv+S8XqE1mP/d89TEi9xob/zH3f T8ENyZb8YuL8N7urZg1RUHZCtiyP1DvEaGw5+4Ws56bK/XzSD026ne1Lhyko aRgVKyUeS7Ad2zBCgeHuyOQm0k/GJka0cJKCVtVyHzMfGuX3G0evnaJgeCxM foZ46AZ0jf2iILZHKqCC9GeG0rqLzqIsMLf2FDHfSePMqKU0V4IFXytya7N2 kfrmbdndMYcFxw4dPh1DzoMCT+unSfNYwLs98C2YnBdujXZR4wtYsHjvs3P2 e2jEG64TD+VYILVnzMgnlMbULne7g4tYMLFgZX8YOY/2KHkV8xRZ0Oehqp1I ziu1WL/AM8rke/X82PoDNIZ7Bn9wY7GgUnW1ufB/ZL+UImp/apF87PVWqrE0 Km2NWnpLh0X+Fxz04k7S+C0mJsJDjwURwYVXJ+NozBqNX1NryIKcJxX051M0 zmpMLkwxZkH2SFmGNDmfq2OupWlvYoFDzt0w78s0plfmjHVbsmBhDbdMI5vG 0NE8m1QrFvj83b1/+AqNbM9bkjPWLCgzKgk4k0NjpHFZ+LPtLPD1+FdJsYDG daO1O319WRAafqKRgyT+wN6eX34seNCxK4f5H+mnzsXuyf4s2DZzeY1KNfHz QIcHgSwQZjtdoh6T+eVKmy0OZcE93736Pg1k/5y3s58fZ8EGM0X5hg80Gjz8 2GFwngUi/SZ1HAkO3r6X5Fh3gQWfP9dvLpvLQf2iVS2eF1lg+SXsfxaSHNTL TKgTXmLBqgdH3h+R4aDOQa2K/mssiHJJLVJU4eBy/oG0jGIWPPI2Xamry0GN s39uE68j8XjC7i4PDmp2GzYn17PAtik6a8ybxNf128ZqYMGYq8EryZ0cFDQ+ 3gYvWDAT6XtlQxAHnSTjbA+9ZkFRhtqzscMcFEbNtx/qYQFn0M98WQoHx/dJ O7ZMs2Au1b7oQx0Hp6v/euMlwobVFy7l/F4xMcZex++ibLjpJleo8ZKDMkVN jlJibJCoSqVPtHOQN3Rmh8U8Nsw8nK7wH+Cgm4+8c+UiNqxZVp0oQ/6XHtst drumxQanp9JSlRZcZGglF9jpsGFl4S+mhRUX3edIiv6px4b8dSOH22y4OHlv +rKvIRs4erclFjhzkcfs7NIwZsPTeKW4oiAuCgdyfK9vZIPfePaGibNctD9u EFTgw4aroi0SR3u52Ju79Wh5OhvOplT02J/j4RtB8KBlBhvu2C+50pDOw/o2 odOHTDZsTDj/0TKLh7clG1ZJZLOhJebLjG0+D6ODN/+wzWVDgvKARXoVD9VX btwzXErW45uztcgwD70em3poPGdD5mXVImUrPn7oNVyfIqoO1Uu+UmsVNNGx FNR6d6vDPu85r8RFtVB4VVXZ8I06mE33BuECbcw6dvOBlfFS0L/HcFJbooM6 A64H2/OWwq0Ueo6trS5mwIM3WlLLYDT9tJxYpB7W3Bl3jg5ZBmGdulzrYn3c jdbfrd4tgxobN+uZKQP8VuSk82GtBkD5kPha7gocMTgZ0VugARHhUyXxfiux 6Yr30WJ5GqQyl9jUR67CQD2FC+GHaHgqntUSkrgao20Ha6t6abj5y9wBso3w iNverLqNHBCP1Px6MnUNfp1SPLXxLgcc95+3HzprjMx759KesLnwcWS2tViE CWp6PrSSPsmFgJpaCVMu4N10OYVno1x4+U+iimwaIBPl4jbv4EGnQ3Ony0/A 9WKfVZoqeRBWEnpSxuEv/NTWXVXL4UNk0MeJ1rK/0C/9lOzJRD7Ee4Sublko wEei7C8WSXxo+TR+30ZagCzf8ifzz/CBJ2Yi1igjwA7tvnBhCh/CXLZ3PpET oE0NfD13kQ+BeY9a7y0WoPG378+uXedDmZr/o7+XCVDW3DG6uoYPl5mBlJyJ APfkf3WPruXD6j9V3OJBgPUyMWvMn5J81U+8RAUCjH1bPPK0ng8fLb2rPpsK cHqvpOeLl3wwTz/8Ci0E+PEimnR28aG+2cjfxU6A62fZL77Sw4eYve9VG+0F mO03OObdxweLzOSLa7cL0EVX6eanfj4kbd07vMxJgC8ehygPD5Px6rRNDboL UEtz7njRdz7wz7aNuHgK8PSZSy9CR/mwlhlS+txLgOYudXHj48QH84QFvgK8 VuXpUz7Jh6VnA66r7RTgLHryryNTfLA7JJ4o9Beg22mhivE0Hz5Vxq8W2SXA yu9LJ2Zm+PBueF5hUKAA/w8pTbcr "]]}}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{ FormBox["\"a\"", TraditionalForm], FormBox["\"LogL\"", TraditionalForm]}, AxesOrigin->{6., -37.}, BaseStyle->{FontSize -> 14}, 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]], ImagePadding->All, Method->{ "DefaultBoundaryStyle" -> Automatic, "DefaultMeshStyle" -> AbsolutePointSize[6], "ScalingFunctions" -> None}, PlotRange->{{5, 15}, {-36.85497105677836, -12.677736833831876`}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, { Scaled[0.05], Scaled[0.05]}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{{3.6705838719155283`*^9, 3.6705838789048815`*^9}, 3.670586207463648*^9, 3.670605104169193*^9, 3.670744858100457*^9, 3.6707449269356995`*^9, 3.673078028530136*^9, 3.6730785780518217`*^9, 3.67307864789723*^9, 3.673079912583273*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["estimating location parameter from the maximum", "Subsection", CellChangeTimes->{{3.670744873303566*^9, 3.6707448933618555`*^9}}], Cell[BoxData[ RowBox[{" ", RowBox[{ RowBox[{"ahat", "[", "data_", "]"}], ":=", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"u", ",", "a"}], "}"}], ",", RowBox[{ RowBox[{"u", "=", RowBox[{"data", "-", "a"}]}], ";", "\[IndentingNewLine]", RowBox[{"a", "/.", RowBox[{"FindRoot", "[", RowBox[{ RowBox[{ RowBox[{"Total", "[", RowBox[{"data", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"u", "^", "2"}]}], ")"}]}], "]"}], "\[Equal]", RowBox[{"a", " ", RowBox[{"Total", "[", RowBox[{"1", "/", RowBox[{"(", RowBox[{"1", "+", RowBox[{"u", "^", "2"}]}], ")"}]}], "]"}]}]}], ",", RowBox[{"{", RowBox[{"a", ",", RowBox[{"Median", "[", "data", "]"}]}], "}"}]}], "]"}]}]}]}], "\[IndentingNewLine]", "]"}]}]}]], "Input", CellChangeTimes->{3.6705856215683546`*^9, 3.670744832436139*^9}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ahat", "[", "w", "]"}]], "Input", CellChangeTimes->{{3.6707449065983953`*^9, 3.6707449104858313`*^9}}], Cell[BoxData["9.383425406724028`"], "Output", CellChangeTimes->{{3.6707449111019163`*^9, 3.6707449295405035`*^9}, 3.673078028561393*^9, 3.673078578073958*^9, 3.673078647912858*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["demonstration", "Subsection", CellChangeTimes->{{3.6707449912106905`*^9, 3.6707449932443447`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"Manipulate", "[", RowBox[{ RowBox[{"randomize", ";", " ", RowBox[{"Module", "[", RowBox[{ RowBox[{"{", RowBox[{"u", ",", "logL", ",", "\[Sigma]0", ",", "min", ",", "max"}], "}"}], ",", RowBox[{ RowBox[{"u", "=", RowBox[{"RandomVariate", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"a0", ",", "1"}], "]"}], ",", "n"}], "]"}]}], " ", ";", "\[IndentingNewLine]", RowBox[{"logL", "=", RowBox[{"logLCauchy", "[", "u", "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"\[Sigma]0", "=", RowBox[{"Sqrt", "[", RowBox[{"2", "/", "n"}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"min", "=", RowBox[{"Max", "[", RowBox[{"0", ",", RowBox[{"a0", "-", RowBox[{"5", "\[Sigma]0"}]}]}], "]"}]}], ";", "\[IndentingNewLine]", RowBox[{"max", " ", "=", " ", RowBox[{"a0", "+", RowBox[{"5", "\[Sigma]0"}]}]}], ";", "\[IndentingNewLine]", RowBox[{"Plot", "[", RowBox[{ RowBox[{"{", " ", RowBox[{ RowBox[{"logL", "[", "a", "]"}], ",", RowBox[{ RowBox[{ RowBox[{"-", RowBox[{ RowBox[{"(", RowBox[{"a", "-", "a0"}], ")"}], "^", "2"}]}], "/", RowBox[{"(", RowBox[{"2", RowBox[{"\[Sigma]0", "^", "2"}]}], ")"}]}], "+", RowBox[{"logL", "[", "a0", "]"}]}]}], "}"}], ",", RowBox[{"{", RowBox[{"a", ",", "min", ",", "max"}], "}"}], ",", RowBox[{"BaseStyle", "\[Rule]", RowBox[{"{", RowBox[{"FontSize", "\[Rule]", "14"}], "}"}]}], ",", RowBox[{"PlotStyle", "\[Rule]", RowBox[{"{", RowBox[{"Blue", ",", "Red"}], "}"}]}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"a", ",", "\"\\""}], "}"}]}], ",", RowBox[{"PlotLabels", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"Epilog", "\[Rule]", RowBox[{"{", RowBox[{"Dashed", ",", "Red", ",", RowBox[{"Line", "[", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"a0", ",", RowBox[{"logL", "[", "max", "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"a0", ",", RowBox[{"logL", "[", "a0", "]"}]}], "}"}]}], "}"}], "]"}]}], "}"}]}]}], "]"}]}]}], "\[IndentingNewLine]", "]"}]}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"a0", ",", "3"}], "}"}], ",", "1", ",", "10"}], "}"}], ",", RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"n", ",", "10"}], "}"}], ",", "1", ",", "100", ",", "1"}], "}"}], ",", RowBox[{"Row", "[", RowBox[{"{", RowBox[{"Button", "[", RowBox[{"\"\\"", ",", RowBox[{"{", RowBox[{"randomize", "=", RowBox[{"RandomReal", "[", "]"}]}], "}"}], ",", RowBox[{"Appearance", "\[Rule]", "\"\\""}], ",", RowBox[{"ImageSize", "\[Rule]", RowBox[{"{", RowBox[{"240", ",", "28"}], "}"}]}]}], "]"}], "}"}], "]"}], ",", RowBox[{"ControlPlacement", "\[Rule]", "Bottom"}]}], "]"}], "\[IndentingNewLine]"}]], "Input", CellChangeTimes->{{3.670745055457469*^9, 3.67074507743085*^9}, { 3.6707452271581335`*^9, 3.670745315957334*^9}, {3.6707454812411165`*^9, 3.6707454846412396`*^9}}], Cell[BoxData[ TagBox[ StyleBox[ DynamicModuleBox[{$CellContext`a0$$ = 3, $CellContext`n$$ = 100, Typeset`show$$ = True, Typeset`bookmarkList$$ = {}, Typeset`bookmarkMode$$ = "Menu", Typeset`animator$$, Typeset`animvar$$ = 1, Typeset`name$$ = "\"untitled\"", Typeset`specs$$ = {{{ Hold[$CellContext`a0$$], 3}, 1, 10}, {{ Hold[$CellContext`n$$], 10}, 1, 100, 1}, { Hold[ Row[{ Button[ "Sort again", {$CellContext`randomize = RandomReal[]}, Appearance -> "Palette", ImageSize -> {240, 28}]}]], Manipulate`Dump`ThisIsNotAControl}}, Typeset`size$$ = { 396., {119., 125.}}, Typeset`update$$ = 0, Typeset`initDone$$, Typeset`skipInitDone$$ = True, $CellContext`a0$3045$$ = 0, $CellContext`n$3050$$ = 0}, DynamicBox[Manipulate`ManipulateBoxes[ 1, StandardForm, "Variables" :> {$CellContext`a0$$ = 3, $CellContext`n$$ = 10}, "ControllerVariables" :> { Hold[$CellContext`a0$$, $CellContext`a0$3045$$, 0], Hold[$CellContext`n$$, $CellContext`n$3050$$, 0]}, "OtherVariables" :> { Typeset`show$$, Typeset`bookmarkList$$, Typeset`bookmarkMode$$, Typeset`animator$$, Typeset`animvar$$, Typeset`name$$, Typeset`specs$$, Typeset`size$$, Typeset`update$$, Typeset`initDone$$, Typeset`skipInitDone$$}, "Body" :> ($CellContext`randomize; Module[{$CellContext`u$, $CellContext`logL$, $CellContext`\[Sigma]0$, \ $CellContext`min$, $CellContext`max$}, $CellContext`u$ = RandomVariate[ CauchyDistribution[$CellContext`a0$$, 1], $CellContext`n$$]; $CellContext`logL$ = \ $CellContext`logLCauchy[$CellContext`u$]; $CellContext`\[Sigma]0$ = Sqrt[2/$CellContext`n$$]; $CellContext`min$ = Max[0, $CellContext`a0$$ - 5 $CellContext`\[Sigma]0$]; $CellContext`max$ = $CellContext`a0$$ + 5 $CellContext`\[Sigma]0$; Plot[{ $CellContext`logL$[$CellContext`a], (-($CellContext`a - \ $CellContext`a0$$)^2)/( 2 $CellContext`\[Sigma]0$^2) + \ $CellContext`logL$[$CellContext`a0$$]}, {$CellContext`a, $CellContext`min$, \ $CellContext`max$}, BaseStyle -> {FontSize -> 14}, PlotStyle -> {Blue, Red}, AxesLabel -> {$CellContext`a, "logL"}, PlotLabels -> {"exp", "limit"}, Epilog -> {Dashed, Red, Line[{{$CellContext`a0$$, $CellContext`logL$[$CellContext`max$]}, {$CellContext`a0$$, $CellContext`logL$[$CellContext`a0$$]}}]}]]), "Specifications" :> {{{$CellContext`a0$$, 3}, 1, 10}, {{$CellContext`n$$, 10}, 1, 100, 1}, Row[{ Button[ "Sort again", {$CellContext`randomize = RandomReal[]}, Appearance -> "Palette", ImageSize -> {240, 28}]}]}, "Options" :> {ControlPlacement -> Bottom}, "DefaultOptions" :> {}], ImageSizeCache->{451., {201., 208.}}, SingleEvaluation->True], Deinitialization:>None, DynamicModuleValues:>{}, SynchronousInitialization->True, UndoTrackedVariables:>{Typeset`show$$, Typeset`bookmarkMode$$}, UnsavedVariables:>{Typeset`initDone$$}, UntrackedVariables:>{Typeset`size$$}], "Manipulate", Deployed->True, StripOnInput->False], Manipulate`InterpretManipulate[1]]], "Output", CellChangeTimes->{3.6707453307491817`*^9, 3.6707454929681907`*^9, 3.6730780285770187`*^9, 3.6730785781052136`*^9, 3.673078648197596*^9}] }, Open ]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["\<\ Comparison between the efficiencies of the Median, trimmed median and Maximum \ Likelihood\ \>", "Section", CellChangeTimes->{{3.6730786967012205`*^9, 3.673078735555438*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"varMaxL", "=", RowBox[{"Table", "[", " ", RowBox[{ RowBox[{"{", RowBox[{"n", ",", RowBox[{"Variance", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"ahat", "[", RowBox[{"RandomVariate", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"0", ",", "1"}], "]"}], ",", "n"}], "]"}], "]"}], ",", RowBox[{"{", "100", "}"}]}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{"n", ",", "15", ",", "30"}], "}"}]}], "]"}]}]], "Input", CellChangeTimes->{{3.673078814856469*^9, 3.673078815593113*^9}, { 3.6730790102850323`*^9, 3.6730790104221554`*^9}, {3.6730790486850405`*^9, 3.673079051546935*^9}, {3.6730790824337673`*^9, 3.6730790894389143`*^9}, { 3.673079140243559*^9, 3.6730792453212748`*^9}}], Cell[BoxData[ RowBox[{"{", RowBox[{ RowBox[{"{", RowBox[{"15", ",", "0.15338881324534537`"}], "}"}], ",", RowBox[{"{", RowBox[{"16", ",", "1030.611809846219`"}], "}"}], ",", RowBox[{"{", RowBox[{"17", ",", "0.12185736975924141`"}], "}"}], ",", RowBox[{"{", RowBox[{"18", ",", "0.1206938655985441`"}], "}"}], ",", RowBox[{"{", RowBox[{"19", ",", "0.15473751471030592`"}], "}"}], ",", RowBox[{"{", RowBox[{"20", ",", "0.10373646400786186`"}], "}"}], ",", RowBox[{"{", RowBox[{"21", ",", "0.10873470743578326`"}], "}"}], ",", RowBox[{"{", RowBox[{"22", ",", "0.11376189695797588`"}], "}"}], ",", RowBox[{"{", RowBox[{"23", ",", "0.12564789746818358`"}], "}"}], ",", RowBox[{"{", RowBox[{"24", ",", "0.09580265089719003`"}], "}"}], ",", RowBox[{"{", RowBox[{"25", ",", "0.12658637132213013`"}], "}"}], ",", RowBox[{"{", RowBox[{"26", ",", "0.09599007668202884`"}], "}"}], ",", RowBox[{"{", RowBox[{"27", ",", "0.0723087198351809`"}], "}"}], ",", RowBox[{"{", RowBox[{"28", ",", "0.05902319907361563`"}], "}"}], ",", RowBox[{"{", RowBox[{"29", ",", "0.07812103440111509`"}], "}"}], ",", RowBox[{"{", RowBox[{"30", ",", "0.08853275662395085`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{ 3.673079007298946*^9, 3.6730790425137043`*^9, 3.6730790912*^9, { 3.673079149362887*^9, 3.6730792531156845`*^9}}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"{", RowBox[{ RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"varMediana", ",", "varMaxL"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0.5"}], "}"}]}], ",", RowBox[{"AxesLabel", "\[Rule]", RowBox[{"{", RowBox[{"\"\\"", ",", "\"\\""}], "}"}]}], ",", RowBox[{"BaseStyle", "\[Rule]", RowBox[{"{", RowBox[{"FontSize", "\[Rule]", "14"}], "}"}]}]}], "]"}], ",", RowBox[{"Plot", "[", RowBox[{ RowBox[{"2", "/", "n"}], ",", RowBox[{"{", RowBox[{"n", ",", "3", ",", "30"}], "}"}]}], "]"}]}], "}"}], "]"}]], "Input", CellChangeTimes->{{3.673079121011117*^9, 3.673079130561592*^9}, 3.6730799390285554`*^9}], Cell[BoxData[ GraphicsBox[{{{}, {{}, {RGBColor[0.368417, 0.506779, 0.709798], PointSize[ 0.012833333333333334`], AbsoluteThickness[1.6], StyleBox[PointBox[CompressedData[" 1:eJxTTMoPSmViYGCQAWIQDQEcDsU3uvYz1es5QPgCDl+EC39mzhaD8kUcuCvu si2/9sEewpdwUOPw7pWwfgHlyzgotAqEpn16DOUrOGxbmvwwlv0BlK/kMNG7 5ufP47egfBWH2fyTt7w6chXKV3O4xtGVJdx+GcrXcChK+BGsW3sBytdy+D9H Nei/yjkoX8fh373iT/svnYby9RyYPrMXp68/BeUbOByru7flwMzjUL6hwyft N3YWt49C+UYOgrEcT2PXHYbyjR02JknP/mcP45s4bH8h8mfuj4NQvqmDrb64 TWkNjG/m8OTGXtl80/1QvrnDp68ljPYqML6Fw1zOR8ULBfZC+ZYOIRwF6+Xn 7YHyrRwM20/XxEvsgvKtHabctXLltYTxbRxEA7XrXmjugPJtHWxOm525Ggjj 2zmY6Mddzfq3zR4AESJ7kw== "]], FontSize->14]}, {RGBColor[0.880722, 0.611041, 0.142051], PointSize[ 0.012833333333333334`], AbsoluteThickness[1.6], StyleBox[ PointBox[{{15., 0.15338881324534537`}, {16., 1030.611809846219}, {17., 0.12185736975924141`}, {18., 0.1206938655985441}, {19., 0.15473751471030592`}, {20., 0.10373646400786186`}, {21., 0.10873470743578326`}, {22., 0.11376189695797588`}, {23., 0.12564789746818358`}, {24., 0.09580265089719003}, {25., 0.12658637132213013`}, {26., 0.09599007668202884}, {27., 0.0723087198351809}, {28., 0.05902319907361563}, {29., 0.07812103440111509}, {30., 0.08853275662395085}}], FontSize->14]}, {}}, {}}, {{{}, {}, {RGBColor[0.368417, 0.506779, 0.709798], AbsoluteThickness[1.6], Opacity[ 1.], LineBox[CompressedData[" 1:eJwVzns40wscgHGPjDSVoa0RJ9dE5n738/0m1FQmnS4TqumCSJ1SilJuKSUl ndJJ6DLmUr+xaElodFEUI5F1jxKdUm4jp/PH+3z+ffV5Uf6bFRUUFIJ/97+7 XFIsVk2oozRw+p7lhjLI1Uv+IDpAw5wCalyBVzcoVyqphmRooMenvZS0M13g ObpxzZNzGhhAuT9TPbULDrvcveqYp4G9OrEbcmO7YLxqn4eaUANLo83PveV1 wffagVhRqwbC9xVVfdZd8KqxY0CFoYn8L2KqUksnVL4qlhZf1ETribLjfcxO CFNZlTsk0MKR3Rbx7fUdkLIlN4VRroWZxbren8UdcKW+L8Lpjhaa2+YIVG50 wOvEBOf9zVroNCFkrs/ugLUKN1omBrVQmyq0iY/qAO/RaRQKMQvNFl8W7NPp AIO+mjDN5ll48qywXHn/c+hqWmBrNUhH+dKMIntOO5RTQoxFcjqyWywF/3q3 Q7rbeYaLEgMb6eP3Sfd2QAFlfNEsBi6i/vLyZrXDlWSZZI0jAy30qBkV09th u1vG6vhYBnamL78ketIGCoLBfU1TZmNEq17R4Io2MEkW10RoMXHnP8xE23Ap vH4h8ynQZWLLX02BQp4UzrOmtL0zYWLQis+VDuukoNax9FOAMxNVFhu95iyT wjezbho7mIk8Y7NF91hSqHo6yTMuYOLefIpb4Y9WWDnHiyJz1cZoj68lPkmt cKismc0J0cG0EEv1OrIFdvyyN9zG10VtLZmDmuczcHHyHz4q+QNTShwF33Wf gurP/SF8FX3sMlvVe0a9CXj1/CVpTAMcnh25OZ/9GJSNJtqoJoZoOrlxznjV IxC60zeuZxuhPPPUCUXiITT6HQuz4xrj8Y3RrVfl96FPry7qZbwJjhnmzsio bYDVsQZ6+nnzMG8NC8cF9cBq7mb3lJri8gu8W2bVEpj+KnDXodb5GHD8gfr8 1nvA/eigfeG9Gf65PaeyTfUelBRKq1PVFuCzmheVI1gHy+TZVE1LC6xJ0Pvm nlYLaz2r2vTcWeib3yNSHKiBbey8ixMhltiQ1ClL4tTAUP6ObyxXK5zTbXRz RdpdWNdrFFpSa4WBMf1cg/RqkFBvB2WttEa237Cf9PgdeC/4RavussZJi2qL I0VVoGJq5OUXZYMt1nId5Ybb4Pcl6LGVgi0miHKOLPwuBs4CLw3/I7aowedv euwoBn+fGc9fMu0w8suRpMKYW0B9nxRQkG+H/b2ajObmSvgsmWlTZ2qP9NhS 85M2lfAkMS7VS2yPTrQizuGMCthSffaiqocDZkd5M+oVKmDK89PWjm0O6MoN kByMuAmrt7klXQ12RJOHcTG5vSJo9h4z431wRNdDIxGzI0WwVCHbL3WvEy6z 7zGc018OB+177L5OOmGg91HVr/vLQbSEfu5upjO6f+wPzFQqB9plLkdm5IJv 6OJP0VllYLjPp54odsGCALupNJsycLbQFSuDK1aI/7blPxKCctDFSPUHrnhi XlLJz0ghjLRaR4f6uuGVxendB+hCGDAtkM9/7YbnL8XFWN0mIXu0EReEEnhg siiKFUaCe/Plk2nhBOZlHpM3bCXh7ZVY2ecIAjlZwbXrt5Bg7mceV7iTwFfP S8dOhZBQde1YhUksgde1wgxGg0joXslm6acTuGp8tkuLPwlzSxt06OUEcj1+ zL/uRoIkMSc8+iaBTZwbL3xdSQjl7rklrSSwZ6SwdsCZBFLJZE3mHQIFS+a6 WDqS4LEu+bT6fQKNR3MbhdYkbJrqqUrtJHDm7oWbJSYkTJXprA1/SWCW5qnT W41JKC4bvPZQ9vuvZ8nwNCMSfgZf9kx9R6BSWNUWf30SUkSKhyj9BJZzPF6+ 0yHBLK2zadNXAmljvHup2iQ0bRDqSr4R6CtOeGPBJIGuxrudMETgkGtaXQyd BPEb52lvRwiM7+Em680iIbiCxl0oJ9DgwOBBiSYJiic+8XMnCIRxn+vhGiTw ebVDk5ME+oWupNNoJPwHUMuuMg== "]]}}, {}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->{True, True}, AxesLabel->{ FormBox["\"n\"", TraditionalForm], FormBox["\"variance\"", TraditionalForm]}, AxesOrigin->{0, 0}, BaseStyle->{FontSize -> 14}, 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]], ImagePadding->All, Method->{"CoordinatesToolOptions" -> {"DisplayFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& ), "CopiedValueFunction" -> ({ (Part[{{Identity, Identity}, {Identity, Identity}}, 1, 2][#]& )[ Part[#, 1]], (Part[{{Identity, Identity}, {Identity, Identity}}, 2, 2][#]& )[ Part[#, 2]]}& )}}, PlotRange->{{0, 30.}, {0, 0.5}}, PlotRangeClipping->True, PlotRangePadding->{{ Scaled[0.02], Scaled[0.02]}, {0, 0}}, Ticks->{Automatic, Automatic}]], "Output", CellChangeTimes->{3.673079134030777*^9, 3.673079232383156*^9, 3.6730793803848567`*^9, 3.6730799396452074`*^9}] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell["Example with n=100", "Section", CellChangeTimes->{{3.67307943406485*^9, 3.6730794385991406`*^9}}], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"varaHat100ML", "=", RowBox[{"Variance", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"ahat", "[", RowBox[{"RandomVariate", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"0", ",", "1"}], "]"}], ",", "100"}], "]"}], "]"}], ",", RowBox[{"{", "1000", "}"}]}], "]"}], "]"}]}]], "Input", CellChangeTimes->{{3.673079446815154*^9, 3.6730794671986313`*^9}, { 3.6730796037119293`*^9, 3.67307960416007*^9}}], Cell[BoxData["0.019717072842624153`"], "Output", CellChangeTimes->{3.6730794868490067`*^9, 3.6730796412787027`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{ RowBox[{"2", "/", "100"}], "//", "N"}]], "Input", CellChangeTimes->{{3.67307949157034*^9, 3.673079496328627*^9}}], Cell[BoxData["0.02`"], "Output", CellChangeTimes->{3.673079496968281*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"varaHatTrimMedian", "=", RowBox[{"Variance", "[", RowBox[{"Table", "[", RowBox[{ RowBox[{"Mean", "[", RowBox[{ RowBox[{"Sort", "[", RowBox[{"RandomVariate", "[", RowBox[{ RowBox[{"CauchyDistribution", "[", RowBox[{"0", ",", "1"}], "]"}], ",", "100"}], "]"}], "]"}], "[", RowBox[{"[", RowBox[{ RowBox[{"Floor", "[", RowBox[{"0.38", " ", "100"}], "]"}], ";;", RowBox[{"Ceiling", "[", RowBox[{"0.62", " ", "100"}], "]"}]}], "]"}], "]"}], " ", "]"}], ",", RowBox[{"{", "1000", "}"}]}], "]"}], "]"}]}]], "Input", CellChangeTimes->{{3.673079547034832*^9, 3.673079547402975*^9}, { 3.6730796304136977`*^9, 3.673079671629961*^9}, {3.673079714901377*^9, 3.673079731809615*^9}}], Cell[BoxData["0.0246341915173824`"], "Output", CellChangeTimes->{{3.6730797083297358`*^9, 3.673079734895982*^9}}] }, Open ]] }, Open ]] }, Open ]] }, WindowSize->{1312, 874}, WindowMargins->{{56, Automatic}, {11, Automatic}}, PrintingCopies->1, PrintingPageRange->{32000, 32000}, PrintingOptions->{"Magnification"->1., "PaperOrientation"->"Portrait", "PaperSize"->{595.1999999999999, 841.68}}, Magnification:>1.1 Inherited, FrontEndVersion->"10.4 for Microsoft Windows (64-bit) (April 11, 2016)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[579, 22, 167, 2, 46, "Subsubtitle"], Cell[CellGroupData[{ Cell[771, 28, 144, 1, 78, "Section"], Cell[CellGroupData[{ Cell[940, 33, 253, 6, 32, "Input"], Cell[1196, 41, 281, 8, 53, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[1514, 54, 714, 16, 53, "Input"], Cell[2231, 72, 96, 1, 33, "Output"] }, Open ]], Cell[2342, 76, 196, 4, 32, "Text"], Cell[CellGroupData[{ Cell[2563, 84, 314, 8, 32, "Input"], Cell[2880, 94, 576, 20, 56, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3493, 119, 213, 4, 32, "Input"], Cell[3709, 125, 460, 15, 56, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4206, 145, 328, 5, 75, "Subsubsection"], Cell[CellGroupData[{ Cell[4559, 154, 624, 17, 32, "Input"], Cell[5186, 173, 110, 2, 48, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[5345, 181, 155, 2, 53, "Subsection"], Cell[CellGroupData[{ Cell[5525, 187, 261, 3, 43, "Subsubsection"], Cell[5789, 192, 112, 1, 32, "Text"], Cell[CellGroupData[{ Cell[5926, 197, 304, 7, 32, "Input"], Cell[6233, 206, 498, 9, 33, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[6768, 220, 127, 2, 32, "Input"], Cell[6898, 224, 130, 2, 33, "Output"] }, Open ]], Cell[7043, 229, 213, 4, 32, "Text"], Cell[7259, 235, 564, 13, 32, "Input"], Cell[CellGroupData[{ Cell[7848, 252, 129, 2, 32, "Input"], Cell[7980, 256, 110, 1, 33, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[8139, 263, 161, 2, 43, "Subsubsection"], Cell[CellGroupData[{ Cell[8325, 269, 674, 18, 32, "Input"], Cell[9002, 289, 2304, 59, 78, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11343, 353, 843, 21, 32, "Input"], Cell[12189, 376, 4237, 86, 274, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[16463, 467, 778, 22, 32, "Input"], Cell[17244, 491, 4158, 83, 282, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[21451, 580, 173, 3, 43, "Subsubsection"], Cell[21627, 585, 324, 6, 32, "Text"], Cell[21954, 593, 365, 10, 32, "Input"], Cell[22322, 605, 305, 6, 32, "Text"], Cell[CellGroupData[{ Cell[22652, 615, 345, 10, 32, "Input"], Cell[23000, 627, 113, 1, 33, "Output"] }, Open ]], Cell[23128, 631, 299, 6, 32, "Text"], Cell[CellGroupData[{ Cell[23452, 641, 500, 15, 32, "Input"], Cell[23955, 658, 129, 2, 33, "Output"] }, Open ]], Cell[24099, 663, 165, 3, 32, "Text"], Cell[CellGroupData[{ Cell[24289, 670, 706, 19, 32, "Input"], Cell[24998, 691, 112, 1, 33, "Output"] }, Open ]], Cell[25125, 695, 150, 2, 32, "Text"], Cell[CellGroupData[{ Cell[25300, 701, 1100, 29, 53, "Input"], Cell[26403, 732, 1756, 43, 78, "Output"] }, Open ]], Cell[28174, 778, 127, 1, 32, "Text"], Cell[CellGroupData[{ Cell[28326, 783, 672, 18, 32, "Input"], Cell[29001, 803, 1761, 43, 78, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[30799, 851, 732, 16, 32, "Input"], Cell[31534, 869, 2458, 57, 267, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[34029, 931, 766, 19, 53, "Input"], Cell[34798, 952, 12919, 264, 277, "Output"] }, Open ]], Cell[47732, 1219, 224, 4, 32, "Text"], Cell[CellGroupData[{ Cell[47981, 1227, 339, 10, 32, "Input"], Cell[48323, 1239, 15172, 304, 277, "Output"] }, Open ]] }, Open ]] }, Open ]] }, Open ]], Cell[63546, 1549, 147, 1, 78, "Section"], Cell[CellGroupData[{ Cell[63718, 1554, 142, 1, 62, "Section"], Cell[CellGroupData[{ Cell[63885, 1559, 109, 1, 53, "Subsection"], Cell[CellGroupData[{ Cell[64019, 1564, 352, 7, 32, "Input"], Cell[64374, 1573, 523, 13, 66, "Output"] }, Open ]], Cell[64912, 1589, 327, 11, 32, "Input"], Cell[CellGroupData[{ Cell[65264, 1604, 256, 6, 32, "Input"], Cell[65523, 1612, 584, 11, 33, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[66144, 1628, 158, 3, 32, "Input"], Cell[66305, 1633, 624, 12, 33, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[66966, 1650, 523, 13, 32, "Input"], Cell[67492, 1665, 6115, 113, 284, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[73656, 1784, 136, 1, 53, "Subsection"], Cell[73795, 1787, 1044, 30, 73, "Input"], Cell[CellGroupData[{ Cell[74864, 1821, 127, 2, 32, "Input"], Cell[74994, 1825, 185, 2, 33, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[75228, 1833, 105, 1, 53, "Subsection"], Cell[CellGroupData[{ Cell[75358, 1838, 3811, 97, 213, "Input"], Cell[79172, 1937, 3493, 70, 427, "Output"] }, Open ]] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[82726, 2014, 185, 4, 78, "Section"], Cell[CellGroupData[{ Cell[82936, 2022, 845, 20, 32, "Input"], Cell[83784, 2044, 1451, 37, 56, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[85272, 2086, 839, 24, 32, "Input"], Cell[86114, 2112, 4866, 95, 282, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[91029, 2213, 104, 1, 78, "Section"], Cell[CellGroupData[{ Cell[91158, 2218, 489, 12, 32, "Input"], Cell[91650, 2232, 116, 1, 33, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[91803, 2238, 140, 3, 32, "Input"], Cell[91946, 2243, 74, 1, 33, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[92057, 2249, 839, 22, 32, "Input"], Cell[92899, 2273, 114, 1, 33, "Output"] }, Open ]] }, Open ]] }, Open ]] } ] *)