(* 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[ 48460, 1299] NotebookOptionsPosition[ 45075, 1177] NotebookOutlinePosition[ 45454, 1193] CellTagsIndexPosition[ 45411, 1190] 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.591944475848307*^9, 3.5919494851885653`*^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.5919453200406866`*^9, 3.591949485810581*^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.591944568806431*^9, 3.5919446016675315`*^9}, 3.5919454183354206`*^9, 3.5919494858225822`*^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.5919446526546326`*^9, 3.5919446601237826`*^9}, 3.591945423037551*^9, 3.5919494858705826`*^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.5919447113140473`*^9, 3.591944732143257*^9}, 3.591944763504307*^9, 3.5919447949591126`*^9, {3.5919448770416245`*^9, 3.5919449035274167`*^9}, 3.591945366421976*^9, 3.5919454279396873`*^9, 3.591949486238591*^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[{ RowBox[{"-", "0.391032142089397`"}], ",", RowBox[{"-", "1.1330178576933734`"}], ",", "0.08615097094802154`", ",", RowBox[{"-", "1.0216794354504193`"}], ",", "0.5731266134399043`", ",", RowBox[{"-", "4.800251447293904`"}], ",", RowBox[{"-", "2.5660235189122464`"}], ",", RowBox[{"-", "1.3008613105143811`"}], ",", "2.6222240095291176`", ",", "1.3297340922896033`"}], "}"}]], "Output", CellChangeTimes->{3.591947843926944*^9, 3.5919494862845936`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Median", "[", "u", "]"}]], "Input", CellChangeTimes->{{3.591947858317342*^9, 3.5919478631295033`*^9}}], Cell[BoxData[ RowBox[{"-", "0.7063557887699081`"}]], "Output", CellChangeTimes->{3.591947863959501*^9, 3.591949486302594*^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.35988448403800766`"], "Output", CellChangeTimes->{{3.5919479705024605`*^9, 3.5919480098095536`*^9}, 3.5919494863565955`*^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", ",", "4.297793531521028`"}], "}"}], ",", RowBox[{"{", RowBox[{"4", ",", "2.9707776894512548`"}], "}"}], ",", RowBox[{"{", RowBox[{"5", ",", "1.1056877723442524`"}], "}"}], ",", RowBox[{"{", RowBox[{"6", ",", "0.7508760458193187`"}], "}"}], ",", RowBox[{"{", RowBox[{"7", ",", "0.610455409832632`"}], "}"}], ",", RowBox[{"{", RowBox[{"8", ",", "0.4829696041213253`"}], "}"}], ",", RowBox[{"{", RowBox[{"9", ",", "0.380848391167129`"}], "}"}], ",", RowBox[{"{", RowBox[{"10", ",", "0.3384488292403201`"}], "}"}], ",", RowBox[{"{", RowBox[{"11", ",", "0.30222984449897755`"}], "}"}], ",", RowBox[{"{", RowBox[{"12", ",", "0.24987717207909843`"}], "}"}], ",", RowBox[{"{", RowBox[{"13", ",", "0.24208906914804718`"}], "}"}], ",", RowBox[{"{", RowBox[{"14", ",", "0.21211977523883627`"}], "}"}], ",", RowBox[{"{", RowBox[{"15", ",", "0.20661417458628592`"}], "}"}], ",", RowBox[{"{", RowBox[{"16", ",", "0.1850042066214656`"}], "}"}], ",", RowBox[{"{", RowBox[{"17", ",", "0.1766017086834271`"}], "}"}], ",", RowBox[{"{", RowBox[{"18", ",", "0.1587747096282328`"}], "}"}], ",", RowBox[{"{", RowBox[{"19", ",", "0.15292895998916317`"}], "}"}], ",", RowBox[{"{", RowBox[{"20", ",", "0.1399867032080739`"}], "}"}], ",", RowBox[{"{", RowBox[{"21", ",", "0.13618850612445108`"}], "}"}], ",", RowBox[{"{", RowBox[{"22", ",", "0.12310090287194443`"}], "}"}], ",", RowBox[{"{", RowBox[{"23", ",", "0.12446939128031967`"}], "}"}], ",", RowBox[{"{", RowBox[{"24", ",", "0.1112067961673509`"}], "}"}], ",", RowBox[{"{", RowBox[{"25", ",", "0.11214296359023762`"}], "}"}], ",", RowBox[{"{", RowBox[{"26", ",", "0.10965998833239621`"}], "}"}], ",", RowBox[{"{", RowBox[{"27", ",", "0.10171919101510997`"}], "}"}], ",", RowBox[{"{", RowBox[{"28", ",", "0.0984565824500083`"}], "}"}], ",", RowBox[{"{", RowBox[{"29", ",", "0.09444945762677892`"}], "}"}], ",", RowBox[{"{", RowBox[{"30", ",", "0.0885920617051254`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{{3.5919456431677713`*^9, 3.5919456657143974`*^9}, 3.5919480735673256`*^9, 3.591949494586824*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"{", RowBox[{ RowBox[{"ListPlot", "[", RowBox[{"varMediana", ",", RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}], ",", 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}}], Cell[BoxData[ GraphicsBox[{{{}, {RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData[" 1:eJxTTMoPSmViYGCQAWIQDQEcDqnxs09+MBB0gPAFHJI574urH2GH8kUcfj6O Wfl0w0d7CF/C4djMFmNd9hdQvozDoe1zTtzseAzlKzgYFa4K//niHpSv5DCH yfjdxfgbUL6Kw09LW3fVlVehfDWHnY/rb+wOvwzlazhkJPAs/fn7PJSv5VBz P8r32J9zUL6OQ7PWviV71U5D+XoO8tVxgaFFp6B8Awezx+1HzNcdh/INHW66 zG58NOsYlG/kwNCw5MyuoCNQvrGD2+OHErqTD0P5Jg78y4WniL4+COWbOhQv /vR3fh6Mb+bgmfsxv6t5P5Rv7tD08tkdy7swvoXD4e0LC3gq9kD5lg4iuyR2 pm2F8a0cOBwm+q8VgvGtHW6KnLzmwrULyrdxuHdvU2SxyU4o39Zhw/GFn6/p 7oDy7Rx08req/1i7zR4AyL2CpA== "]]}, {}}, {{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwV0Gk01AsYBnCHibrINkZDcq0xlV2W8fe+Ksot0aRFZQllTwslriVaZKkp lJu6ZRSSpT8GV4qxxrEbEjVHaHJz1EQryXU/POf58Pv0PJo+Yawj4mJiYp5L +b9zVYU2exbkcY94p6eztgBYNVv/5cYo4L7H9REFDq+hIzhR1petiCK39l/J GSPgeFx4sDNLEUdnaCvkk0agPtyp0DJHEWOuP3e/Fz0C5TFyW2XKFLFXKfza mM8I3LqaHc/tV0R/s2vcKZMRCCgvn5VSUUIhd5JC6RsGyvz4y6I7Ssg523hh ij4MzOQteV8LqcjUlY8cbB6Cvw1upKpUUHE5Y9Hufc0QiLe9O2n1lIqPE2Jz pR4PwXOpFIjqpuL2Rbai160h2Huhd2hhloqVhSnr4sKG4GS8p/QyQhmN58Jz z6oNQeHpyONK3cpYmnKmWDLqBaj6FtkYz9KQ2O5438JlEI6FTa7nztPQhf2W I3IcBF60joYNRQW/ZZL1pN0gBGTcltiirIKNSuV2joaDUNmc2rHPcsnd6pKq ZAeBpR/qGRetguxUsSxu5wAkT29I6JJYhc8me3Jndw3A3JnSthAqHSf1u6PN gvhgqTG8u0CdjoNW5/eU+fAhvIUiGNej41XWPXLjQT58pB6YOWBNx9DG4SGX HXyYIClqTp50ZDx6w2w05EPXlHuIbgEduxLrzB9+7geOF0VOwFTFJGFJ3h/n +8Fpq7ubi68a1in/Wt5A9oFS3w6j4Hx1PK7NN5LZ0guCrBFKcpMGml11z51R 7wHZtOXB+VKauCtjcThDvgseTYlcU+haeEQiaj/HqQOk9H8JpPW0UXC4W+5n bTtk8wRBXk46aFR35Jw40QZGP2TCzd118WZxXuuD+VaYT1sZ/SpOD+9Umkqw eS1wQcyCoZmzFhM8V1v8LGwGOYWHe9+V6KOfOaeE8awJemNF8fH9Blh87ZCk QX8jpCdf0M+eYKDPpZiSgRWNkDY73Zkksx55TpKl37EBHGyi1ZSMNmD9tvR3 dik8WFizV7jGzhAHMm4XiX+oh2lvSumCrxHmV/cOnHeph7Df4+cMmca4VeXy /V0pdeAYxz5RzDNGrYW6bVpXnkF1arJ/5m4T1J1jIj/1KawLXLf62YgJVsQG qF96VAs2n7kurmGmyCFXyki2PIG1H1gvjMXMMLiGGmk/UwPKVTfUWZfMMDsj Yk+HZQ3cNWwae0U3x9eryIiHkf+AQUTW0QKOOYruX1ne3V0N1xlxdg36Fnjz R4T6VdNquHHFJdOhxgKtoxrsz7Grlv5tKVyxaSO2Bz1Y0SxWBba5a8FyYCP+ HJSoig2phNe7ja4/8LREDY3A4HuTXPA45WTt89YS9We1vVeFcuGsvbNX0hkr nEh1U1k9XQGZgm+bPy5aYdSAycKHqApwZqTm16VbI9uR75pOqYDxCVMvgY4N OgiHX0dklkPUsfgBosgGVzde/ilvWg7FeVXtksDEMWMd/fz2MhgtzoyVf87E QNjE+RJaBjPaFy8G7LTFteHRfTG0Mri7zVuWMWqLt0yzQ42fkOBs+4a1PoDA HMY3P8NAEoS1x3NSggg0dRD/1OJPQjwhJnofQiC1TrzK6ygJZXYaaQ9PEBjT +kh0zZeEVfYerXrRBG4S0ug/PEiYcHjJ1LxCoPdHXdM+FglRLr06tAoCqxW4 OqW2JFB7vE9FVBJ4cjKbv5NJQrGriMevJtDKub72gzUJo7vkvNKfEnh/KsXC yJIEBzfnv+RbCVx/s7K1zIQEefc2WelhAn0SC3ya9EgofLn/UNArAueoWmx/ XRI2H5gsbBMQuLJL4stvOiScPii1LWmcwB+idj+WJgmvPBwSlk0TmPTnt+Fx tSUX8Hv8PhK4rtS+IUmVBDkvP42mT0v+Pnl0A52ETd6JtQlfl/YzhbxIGgkj o4rSY98JDLERnF+jTEL4YY67/TyBHbSc2CYlElaOmRTcWyAwq1uzNEiRhHwf 3tfFRQJdA3bTFBRI+A+If7KD "]]}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{3., 0.0885920617051254}, ImageSize->{729.9999705950452, Automatic}, Method->{}, PlotRangeClipping->True]], "Output", CellChangeTimes->{{3.591945730274194*^9, 3.591945746546646*^9}, 3.59194808380361*^9, 3.5919481271148157`*^9, 3.5919494136625757`*^9, 3.591949494726829*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"Show", "[", RowBox[{"{", RowBox[{ RowBox[{"ListPlot", "[", RowBox[{"varMediana", ",", RowBox[{"PlotRange", "\[Rule]", RowBox[{"{", RowBox[{"0", ",", "0.5"}], "}"}]}]}], "]"}], ",", RowBox[{"Plot", "[", RowBox[{ RowBox[{"2", "/", "n"}], ",", RowBox[{"{", RowBox[{"n", ",", "3", ",", "30"}], "}"}]}], "]"}]}], "}"}], "]"}]], "Input", CellChangeTimes->{{3.5919494440894194`*^9, 3.591949449883581*^9}}], Cell[BoxData[ GraphicsBox[{{{}, {RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData[" 1:eJxTTMoPSmViYGCQAWIQDQEcDqnxs09+MBB0gPAFHJI574urH2GH8kUcfj6O Wfl0w0d7CF/C4djMFmNd9hdQvozDoe1zTtzseAzlKzgYFa4K//niHpSv5DCH yfjdxfgbUL6Kw09LW3fVlVehfDWHnY/rb+wOvwzlazhkJPAs/fn7PJSv5VBz P8r32J9zUL6OQ7PWviV71U5D+XoO8tVxgaFFp6B8Awezx+1HzNcdh/INHW66 zG58NOsYlG/kwNCw5MyuoCNQvrGD2+OHErqTD0P5Jg78y4WniL4+COWbOhQv /vR3fh6Mb+bgmfsxv6t5P5Rv7tD08tkdy7swvoXD4e0LC3gq9kD5lg4iuyR2 pm2F8a0cOBwm+q8VgvGtHW6KnLzmwrULyrdxuHdvU2SxyU4o39Zhw/GFn6/p 7oDy7Rx08req/1i7zR4AyL2CpA== "]]}, {}}, {{}, {}, {Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwV0Gk01AsYBnCHibrINkZDcq0xlV2W8fe+Ksot0aRFZQllTwslriVaZKkp lJu6ZRSSpT8GV4qxxrEbEjVHaHJz1EQryXU/POf58Pv0PJo+Yawj4mJiYp5L +b9zVYU2exbkcY94p6eztgBYNVv/5cYo4L7H9REFDq+hIzhR1petiCK39l/J GSPgeFx4sDNLEUdnaCvkk0agPtyp0DJHEWOuP3e/Fz0C5TFyW2XKFLFXKfza mM8I3LqaHc/tV0R/s2vcKZMRCCgvn5VSUUIhd5JC6RsGyvz4y6I7Ssg523hh ij4MzOQteV8LqcjUlY8cbB6Cvw1upKpUUHE5Y9Hufc0QiLe9O2n1lIqPE2Jz pR4PwXOpFIjqpuL2Rbai160h2Huhd2hhloqVhSnr4sKG4GS8p/QyQhmN58Jz z6oNQeHpyONK3cpYmnKmWDLqBaj6FtkYz9KQ2O5438JlEI6FTa7nztPQhf2W I3IcBF60joYNRQW/ZZL1pN0gBGTcltiirIKNSuV2joaDUNmc2rHPcsnd6pKq ZAeBpR/qGRetguxUsSxu5wAkT29I6JJYhc8me3Jndw3A3JnSthAqHSf1u6PN gvhgqTG8u0CdjoNW5/eU+fAhvIUiGNej41XWPXLjQT58pB6YOWBNx9DG4SGX HXyYIClqTp50ZDx6w2w05EPXlHuIbgEduxLrzB9+7geOF0VOwFTFJGFJ3h/n +8Fpq7ubi68a1in/Wt5A9oFS3w6j4Hx1PK7NN5LZ0guCrBFKcpMGml11z51R 7wHZtOXB+VKauCtjcThDvgseTYlcU+haeEQiaj/HqQOk9H8JpPW0UXC4W+5n bTtk8wRBXk46aFR35Jw40QZGP2TCzd118WZxXuuD+VaYT1sZ/SpOD+9Umkqw eS1wQcyCoZmzFhM8V1v8LGwGOYWHe9+V6KOfOaeE8awJemNF8fH9Blh87ZCk QX8jpCdf0M+eYKDPpZiSgRWNkDY73Zkksx55TpKl37EBHGyi1ZSMNmD9tvR3 dik8WFizV7jGzhAHMm4XiX+oh2lvSumCrxHmV/cOnHeph7Df4+cMmca4VeXy /V0pdeAYxz5RzDNGrYW6bVpXnkF1arJ/5m4T1J1jIj/1KawLXLf62YgJVsQG qF96VAs2n7kurmGmyCFXyki2PIG1H1gvjMXMMLiGGmk/UwPKVTfUWZfMMDsj Yk+HZQ3cNWwae0U3x9eryIiHkf+AQUTW0QKOOYruX1ne3V0N1xlxdg36Fnjz R4T6VdNquHHFJdOhxgKtoxrsz7Grlv5tKVyxaSO2Bz1Y0SxWBba5a8FyYCP+ HJSoig2phNe7ja4/8LREDY3A4HuTXPA45WTt89YS9We1vVeFcuGsvbNX0hkr nEh1U1k9XQGZgm+bPy5aYdSAycKHqApwZqTm16VbI9uR75pOqYDxCVMvgY4N OgiHX0dklkPUsfgBosgGVzde/ilvWg7FeVXtksDEMWMd/fz2MhgtzoyVf87E QNjE+RJaBjPaFy8G7LTFteHRfTG0Mri7zVuWMWqLt0yzQ42fkOBs+4a1PoDA HMY3P8NAEoS1x3NSggg0dRD/1OJPQjwhJnofQiC1TrzK6ygJZXYaaQ9PEBjT +kh0zZeEVfYerXrRBG4S0ug/PEiYcHjJ1LxCoPdHXdM+FglRLr06tAoCqxW4 OqW2JFB7vE9FVBJ4cjKbv5NJQrGriMevJtDKub72gzUJo7vkvNKfEnh/KsXC yJIEBzfnv+RbCVx/s7K1zIQEefc2WelhAn0SC3ya9EgofLn/UNArAueoWmx/ XRI2H5gsbBMQuLJL4stvOiScPii1LWmcwB+idj+WJgmvPBwSlk0TmPTnt+Fx tSUX8Hv8PhK4rtS+IUmVBDkvP42mT0v+Pnl0A52ETd6JtQlfl/YzhbxIGgkj o4rSY98JDLERnF+jTEL4YY67/TyBHbSc2CYlElaOmRTcWyAwq1uzNEiRhHwf 3tfFRQJdA3bTFBRI+A+If7KD "]]}}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{0, 0}, Method->{}, PlotRange->{{0, 30.}, {0, 0.5}}, PlotRangeClipping->True, PlotRangePadding->{{0.6, 0.6}, {0., 0.}}]], "Output", CellChangeTimes->{3.591949450263593*^9, 3.59194949475883*^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[ RowBox[{"-", "0.1276922692171784`"}]], "Output", CellChangeTimes->{3.5919484325693054`*^9, 3.5919494947788305`*^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.10844626371594702`"}]], "Output", CellChangeTimes->{{3.591948523321827*^9, 3.591948531240048*^9}, 3.5919494947948313`*^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.011762077577393398`"], "Output", CellChangeTimes->{3.5919461369514956`*^9, 3.5919487099350147`*^9, 3.5919494948568316`*^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.021699160104204855`"}], "}"}], ",", RowBox[{"{", RowBox[{"200", ",", "0.012405650364490529`"}], "}"}], ",", RowBox[{"{", RowBox[{"300", ",", "0.007425972443160326`"}], "}"}], ",", RowBox[{"{", RowBox[{"400", ",", "0.006116947605791854`"}], "}"}], ",", RowBox[{"{", RowBox[{"500", ",", "0.00454146039140622`"}], "}"}], ",", RowBox[{"{", RowBox[{"600", ",", "0.0038910522839241232`"}], "}"}], ",", RowBox[{"{", RowBox[{"700", ",", "0.0032704540061168305`"}], "}"}], ",", RowBox[{"{", RowBox[{"800", ",", "0.002679679246756394`"}], "}"}], ",", RowBox[{"{", RowBox[{"900", ",", "0.002441129398048492`"}], "}"}], ",", RowBox[{"{", RowBox[{"1000", ",", "0.002279613069110258`"}], "}"}], ",", RowBox[{"{", RowBox[{"1100", ",", "0.00206993325400291`"}], "}"}], ",", RowBox[{"{", RowBox[{"1200", ",", "0.002080924067152367`"}], "}"}], ",", RowBox[{"{", RowBox[{"1300", ",", "0.0017372578067082946`"}], "}"}], ",", RowBox[{"{", RowBox[{"1400", ",", "0.0015478255788381379`"}], "}"}], ",", RowBox[{"{", RowBox[{"1500", ",", "0.0014114447744974846`"}], "}"}], ",", RowBox[{"{", RowBox[{"1600", ",", "0.0013315536696652584`"}], "}"}], ",", RowBox[{"{", RowBox[{"1700", ",", "0.001321943990767372`"}], "}"}], ",", RowBox[{"{", RowBox[{"1800", ",", "0.0013057869457953558`"}], "}"}], ",", RowBox[{"{", RowBox[{"1900", ",", "0.0012091800018010559`"}], "}"}], ",", RowBox[{"{", RowBox[{"2000", ",", "0.001147552489049451`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{3.591946191851022*^9, 3.5919462243439255`*^9, 3.591946301208062*^9, 3.591946417403291*^9, 3.5919488237101774`*^9, 3.5919494974329033`*^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.024993928441981183`"}], "}"}], ",", RowBox[{"{", RowBox[{"200", ",", "0.011942487129775654`"}], "}"}], ",", RowBox[{"{", RowBox[{"300", ",", "0.007799257692211399`"}], "}"}], ",", RowBox[{"{", RowBox[{"400", ",", "0.005991055725440552`"}], "}"}], ",", RowBox[{"{", RowBox[{"500", ",", "0.005278901701437471`"}], "}"}], ",", RowBox[{"{", RowBox[{"600", ",", "0.004074782838774228`"}], "}"}], ",", RowBox[{"{", RowBox[{"700", ",", "0.0035803476956317868`"}], "}"}], ",", RowBox[{"{", RowBox[{"800", ",", "0.003022823968717223`"}], "}"}], ",", RowBox[{"{", RowBox[{"900", ",", "0.0026768869320222175`"}], "}"}], ",", RowBox[{"{", RowBox[{"1000", ",", "0.002500757624006736`"}], "}"}], ",", RowBox[{"{", RowBox[{"1100", ",", "0.0021244285579737647`"}], "}"}], ",", RowBox[{"{", RowBox[{"1200", ",", "0.00217626689804539`"}], "}"}], ",", RowBox[{"{", RowBox[{"1300", ",", "0.001972640876998141`"}], "}"}], ",", RowBox[{"{", RowBox[{"1400", ",", "0.0017003832254462886`"}], "}"}], ",", RowBox[{"{", RowBox[{"1500", ",", "0.0017570509233490069`"}], "}"}], ",", RowBox[{"{", RowBox[{"1600", ",", "0.001574429678393771`"}], "}"}], ",", RowBox[{"{", RowBox[{"1700", ",", "0.0015134847610835831`"}], "}"}], ",", RowBox[{"{", RowBox[{"1800", ",", "0.0014024546393668737`"}], "}"}], ",", RowBox[{"{", RowBox[{"1900", ",", "0.001349529472864001`"}], "}"}], ",", RowBox[{"{", RowBox[{"2000", ",", "0.0011959983528175083`"}], "}"}]}], "}"}]], "Output", CellChangeTimes->{{3.5919488971962185`*^9, 3.591948907194496*^9}, 3.5919494999869757`*^9}] }, Open ]], Cell[CellGroupData[{ Cell[BoxData[ RowBox[{"ListPlot", "[", RowBox[{ RowBox[{"{", RowBox[{"varMedianaNGrande", ",", "varEstim24"}], "}"}], ",", RowBox[{"PlotRange", "\[Rule]", "All"}]}], "]"}]], "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}}], Cell[BoxData[ GraphicsBox[{{}, { {RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData[" 1:eJxTTMoPSmViYGAQAWIQDQGRDpc52VWZZsy0h/AzHbzY7Wp9Szsg/ANFDrt3 dnrLfqqHylc6fMrsOSbcVQHhO9Q7uN1fmRQyvxSqvsnBz7/OYMeGAgj/QauD U0NUR0dILlR/p8PvpFNOQYczIHyFHodshoZd196mQs3rdwjh3M1WUZMC4RtM dDj2ePdT+/REqPmTHN7mf0z5dhnKD5ji4PZxa+NbjQSofVMdTmUzT3tzOxrC L5ju4O1mIjD7ZAzU/pkOn+cs+LfodCSEP2GWw5nqwMfMZyKg7pnjsCXplnHQ rzAIf8NcBxf9OzZhMlC+w3yHhUHvihNnBNsDAFpCY3w= "]]}, {RGBColor[0.6, 0.24, 0.4428931686004542], PointBox[CompressedData[" 1:eJxTTMoPSmViYGAQAWIQDQGRDnn6/377Wkyzh/AzHaR4zpsoZXRC+AeKHLJO 9PhsyqqDylc6cHNunCfKVwnhO9Q7vP627a/ArCKo+iaHmevPnVR4kA/hP2h1 MF3wa/7iU1lQ/Z0OEjMDrm/8nArhK/Q4bHlzMDb7fzLUvH6HfcpdnzauSYLw DSY6zHmk/+XblwSo+ZMcUr25ZzHzJEL4AVMcbjzfFTutLAZq31SH689z+Hlj IiH8gukO85sTPdkVwqH2z3SQbNx898OFUAh/wiyHzCb33kUroHyFOQ5PTCbz 3UuB8jfMdYitYVZZfT4Y6r75DkdjfKT0zwbZAwDvo2PA "]]}}, {}}, AspectRatio->NCache[GoldenRatio^(-1), 0.6180339887498948], Axes->True, AxesLabel->{None, None}, AxesOrigin->{100., 0.001147552489049451}, ImageSize->{419.3333286179448, Automatic}, Method->{}, PlotRangeClipping->True]], "Output", CellChangeTimes->{ 3.591946226817993*^9, 3.59194630440215*^9, 3.59194639031254*^9, { 3.5919464208733883`*^9, 3.591946437257845*^9}, {3.5919464964234915`*^9, 3.591946559909255*^9}, 3.5919490016071224`*^9, {3.5919490543785887`*^9, 3.5919490744171467`*^9}, 3.59194911446826*^9, 3.591949146283144*^9, 3.5919495000249767`*^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[{"\"\\"", ",", "\"\\""}], "}"}]}]}], "]"}]}]], "Input", CellChangeTimes->{{3.5919491691237764`*^9, 3.5919491724958706`*^9}}], Cell[BoxData[ TemplateBox[{GraphicsBox[{{}, {{ RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData[" 1:eJxTTMoPSmViYGAQAWIQDQGRDpc52VWZZsy0h/AzHbzY7Wp9Szsg/ANFDrt3 dnrLfqqHylc6fMrsOSbcVQHhO9Q7uN1fmRQyvxSqvsnBz7/OYMeGAgj/QauD U0NUR0dILlR/p8PvpFNOQYczIHyFHodshoZd196mQs3rdwjh3M1WUZMC4RtM dDj2ePdT+/REqPmTHN7mf0z5dhnKD5ji4PZxa+NbjQSofVMdTmUzT3tzOxrC L5ju4O1mIjD7ZAzU/pkOn+cs+LfodCSEP2GWw5nqwMfMZyKg7pnjsCXplnHQ rzAIf8NcBxf9OzZhMlC+w3yHhUHvihNnBNsDAFpCY3w= "]]}, { RGBColor[0.6, 0.24, 0.4428931686004542], PointBox[CompressedData[" 1:eJxTTMoPSmViYGAQAWIQDQGRDnn6/377Wkyzh/AzHaR4zpsoZXRC+AeKHLJO 9PhsyqqDylc6cHNunCfKVwnhO9Q7vP627a/ArCKo+iaHmevPnVR4kA/hP2h1 MF3wa/7iU1lQ/Z0OEjMDrm/8nArhK/Q4bHlzMDb7fzLUvH6HfcpdnzauSYLw DSY6zHmk/+XblwSo+ZMcUr25ZzHzJEL4AVMcbjzfFTutLAZq31SH689z+Hlj IiH8gukO85sTPdkVwqH2z3SQbNx898OFUAh/wiyHzCb33kUroHyFOQ5PTCbz 3UuB8jfMdYitYVZZfT4Y6r75DkdjfKT0zwbZAwDvo2PA "]]}}, {}}, AspectRatio -> NCache[GoldenRatio^(-1), 0.6180339887498948], Axes -> True, AxesLabel -> {None, None}, AxesOrigin -> {0, 0}, Method -> {}, PlotRange -> {{0, 2000.}, {0, 0.004}}, PlotRangeClipping -> True, PlotRangePadding -> {{40., 40.}, {0., 0.}}], TemplateBox[{"\"mediana\"", "\"media central\""}, "PointLegend", DisplayFunction -> (StyleBox[ StyleBox[ PaneBox[ TagBox[ GridBox[{{ TagBox[ GridBox[{{ GraphicsBox[{{}, { Directive[ EdgeForm[{ Opacity[0.3], GrayLevel[0]}], RGBColor[0.24720000000000014`, 0.24, 0.6]], { Directive[ EdgeForm[{ Opacity[0.3], GrayLevel[0]}], RGBColor[0.24720000000000014`, 0.24, 0.6]], PointBox[ NCache[{ Scaled[{ Rational[1, 2], Rational[1, 2]}]}, { Scaled[{0.5, 0.5}]}]]}}}, AspectRatio -> Full, ImageSize -> {10, 10}, PlotRangePadding -> None, ImagePadding -> 1, BaselinePosition -> (Scaled[0.1] -> Baseline)], #}, { GraphicsBox[{{}, { Directive[ EdgeForm[{ Opacity[0.3], GrayLevel[0]}], RGBColor[0.6, 0.24, 0.4428931686004542]], { Directive[ EdgeForm[{ Opacity[0.3], GrayLevel[0]}], RGBColor[0.6, 0.24, 0.4428931686004542]], PointBox[ NCache[{ Scaled[{ Rational[1, 2], Rational[1, 2]}]}, { Scaled[{0.5, 0.5}]}]]}}}, AspectRatio -> Full, ImageSize -> {10, 10}, PlotRangePadding -> None, ImagePadding -> 1, 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 -> "Times"}, Background -> Automatic, StripOnInput -> False]& ), Editable -> True, InterpretationFunction :> (RowBox[{"PointLegend", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Directive", "[", RowBox[{"RGBColor", "[", RowBox[{"0.24720000000000014`", ",", "0.24`", ",", "0.6`"}], "]"}], "]"}], ",", RowBox[{"Directive", "[", RowBox[{"RGBColor", "[", RowBox[{"0.6`", ",", "0.24`", ",", "0.4428931686004542`"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{#, ",", #2}], "}"}], ",", RowBox[{"LegendLayout", "\[Rule]", "\"Column\""}], ",", RowBox[{"LegendMarkers", "\[Rule]", "False"}]}], "]"}]& )]}, "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.5919491554914*^9, 3.5919491729878855`*^9}, 3.5919495001029787`*^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[{{{}, {{ RGBColor[0.24720000000000014`, 0.24, 0.6], PointBox[CompressedData[" 1:eJxTTMoPSmViYGAQAWIQDQGRDpc52VWZZsy0h/AzHbzY7Wp9Szsg/ANFDrt3 dnrLfqqHylc6fMrsOSbcVQHhO9Q7uN1fmRQyvxSqvsnBz7/OYMeGAgj/QauD U0NUR0dILlR/p8PvpFNOQYczIHyFHodshoZd196mQs3rdwjh3M1WUZMC4RtM dDj2ePdT+/REqPmTHN7mf0z5dhnKD5ji4PZxa+NbjQSofVMdTmUzT3tzOxrC L5ju4O1mIjD7ZAzU/pkOn+cs+LfodCSEP2GWw5nqwMfMZyKg7pnjsCXplnHQ rzAIf8NcBxf9OzZhMlC+w3yHhUHvihNnBNsDAFpCY3w= "]]}, { RGBColor[0.6, 0.24, 0.4428931686004542], PointBox[CompressedData[" 1:eJxTTMoPSmViYGAQAWIQDQGRDnn6/377Wkyzh/AzHaR4zpsoZXRC+AeKHLJO 9PhsyqqDylc6cHNunCfKVwnhO9Q7vP627a/ArCKo+iaHmevPnVR4kA/hP2h1 MF3wa/7iU1lQ/Z0OEjMDrm/8nArhK/Q4bHlzMDb7fzLUvH6HfcpdnzauSYLw DSY6zHmk/+XblwSo+ZMcUr25ZzHzJEL4AVMcbjzfFTutLAZq31SH689z+Hlj IiH8gukO85sTPdkVwqH2z3SQbNx898OFUAh/wiyHzCb33kUroHyFOQ5PTCbz 3UuB8jfMdYitYVZZfT4Y6r75DkdjfKT0zwbZAwDvo2PA "]]}}, {}}, {{}, {}, { Hue[0.67, 0.6, 0.6], LineBox[CompressedData[" 1:eJwVxXk4lAkcAGBJRbVlVbrIsjoxOjxNT9jfr8Iqsp9zHONoDBMK45vPmGlK coyjnOU2H6Kl6EHJEquDko0NiRzlqOiSJGnaanf/eJ9XmxNk56OooKDg8Z// r3Jz1Uys4uNI09uyigkJqHMcN1nuCsEVoW2+L1dLwCp5GYPnRGLapuUqWmIx TF+tAz1PEt00mp7+HCgG2eNDtu94JF705sbt4IhhUqdCIBSR6PzXtzXkATGk V1tfl+aSaBapWsRfK4aRvmjLkhESGfxrSl6NIhCvn+W8DBTgnHDP0qIlIiit G8gIkFJYZT2weEuLENgBGoNvkykce6VZXt4ghCUabJ3gbAqZpCh29xUhhBzv LxOUUyhyG2sMkgnBBPtuHu+iUPOy1ogFJYT25p7XSVqhGBt+/nqerhCmHnTC 1ZpQ1DeqUMyPCwWTsZaxf8aEuPa2bv8tfwr+XBdSOzApRCX6hbSTQwE4rU2o /yzEJvb6FeOuFOxtPsqQLApDP0KtRteKgv3n1agvjDCcqUnkjupR4OTJVpRT YchI6RkNnhBAcM87jU9zRchVHtxZKBRA8V11u/daYvxqbmC4MYuErQnOrKTN YgyuFAxPpZBQZ5PNNtghRv2CPwZvx5PQ3q3J87MQ4+LP8lxKQsLsqK5k+IgY 1/GLz630IsFaYceFjloxTjxLjCI2kTC96zd5pf0xPKWqlB9cHwJ7L0oLQ6QS fGRa9tR+ig/drua2daMnUH7A5MUvacFAe24/U5VwEufdDBvseB8Ib1O7uDYN ETguEg/UDByBDfnF93L6TuH6E0pzZzcGgNtymbv1VCTyQ5wCDC/5QalKCfe2 djQq339xO2X3YRjpHZl9bRKDYxVj95bO5wFzs6FCf5AUe7NPYXyvD0QPFew0 yIvF1m1my9e2caGRby7fei0OM5bppd166Q3LouRbv7yIx3GOGWejsjfYBQ9N h849jfdszPOGLThQNtijV8s8g4xHnYIfYg7BnKqPKR6+idgRyD2d2eMF6GUW XX0yCd/5mB2+gF7Q9msn5eeSjI/26S/6TnhC8P3sd8R4Mh7PifSysfSA5KgZ BX9BClp8/V4Aju5g+rXPyndxKvrIaV8TDzYMjbDcjbNTsfnbJZepaDdQjXMv +bAjDUsM63l0qSs0cGbiB1vTkLyb6GD/xAW8Sjca5bPOone5abS9ngskVwyr W0yexQ/N9HA36QzVukY3mk6cwyPWqyLrO1jgfH/1HGphOu5nZj+8s5sFsUax x5QL0/F6v/phI9oJylWS6gWGGTh9fq/KrhVO4EmvYza1ZOBYbd+ioRhHWFBo EfaGlYmimm0rLJUdQb93jq3Fs0xcEKmjUJniAE/bzb0jJVmo/zhzf5e6A0i3 65hkLczGPdFCF2apPdQ7SW1Si7LRP+FmaQbTHoIuh8/WG+dg5UGG8pv7drCS VPW70ZKD8c1BUyeP2oH17OZlvx/KRd6+JWq7F9qB8FV/pO9kLl5UOV5nWmkL Wqrap1mSPDz0aiDd3N0WHGYm2vSWy3CiMV6bnmcL4fsdH/UWyJBIau9/do0A QVB3iG2RDOVKF1oSqgnwP+e4tPWCDPe5ab/dfpUAh2FHy+uXZMh80jYQUUnA JpFTXV61DLPoGEedMgIelLBknHsynN7Tkn+0gAAdZVffN5MyLI74uMfgNAEr GX2KPh9kGBGYoNIfT8BiB1d68KMMQxS1FsTFETBDu/a2y2WYyl9U9DyGgFam 24FKJRq1UpK9CyIIIHlsRugqGov91CZ0hAQcPjPQOrGGxskhydKHFAEeV9g8 niaNn7RX+UYLCLD8zs531qFxlnW5apxPgEaGu5qxPo3VXZ/uVh4h4MeGwctX GDSer49J4QYQMH/U3Up/G41fo4aTVvoTMMnwiNLcSWNhfMuWEzwCnjs8+Sl9 F41GDcZ3t/kS0Cf2aFhiTGNtrWXmcy4Bf+c/cZGa0linNJqX5U1A0x2Pme9A Y9WGbz0HOQT8CxyMv/k= "]]}}}, AspectRatio -> NCache[GoldenRatio^(-1), 0.6180339887498948], Axes -> True, AxesLabel -> {None, None}, AxesOrigin -> {0, 0}, Method -> {}, PlotRange -> {{0, 2000.}, {0, 0.004}}, PlotRangeClipping -> True, PlotRangePadding -> {{40., 40.}, {0., 0.}}], TemplateBox[{"\"mediana\"", "\"media central\""}, "PointLegend", DisplayFunction -> (StyleBox[ StyleBox[ PaneBox[ TagBox[ GridBox[{{ TagBox[ GridBox[{{ GraphicsBox[{{}, { Directive[ EdgeForm[{ Opacity[0.3], GrayLevel[0]}], RGBColor[0.24720000000000014`, 0.24, 0.6]], { Directive[ EdgeForm[{ Opacity[0.3], GrayLevel[0]}], RGBColor[0.24720000000000014`, 0.24, 0.6]], PointBox[ NCache[{ Scaled[{ Rational[1, 2], Rational[1, 2]}]}, { Scaled[{0.5, 0.5}]}]]}}}, AspectRatio -> Full, ImageSize -> {10, 10}, PlotRangePadding -> None, ImagePadding -> 1, BaselinePosition -> (Scaled[0.1] -> Baseline)], #}, { GraphicsBox[{{}, { Directive[ EdgeForm[{ Opacity[0.3], GrayLevel[0]}], RGBColor[0.6, 0.24, 0.4428931686004542]], { Directive[ EdgeForm[{ Opacity[0.3], GrayLevel[0]}], RGBColor[0.6, 0.24, 0.4428931686004542]], PointBox[ NCache[{ Scaled[{ Rational[1, 2], Rational[1, 2]}]}, { Scaled[{0.5, 0.5}]}]]}}}, AspectRatio -> Full, ImageSize -> {10, 10}, PlotRangePadding -> None, ImagePadding -> 1, 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 -> "Times"}, Background -> Automatic, StripOnInput -> False]& ), Editable -> True, InterpretationFunction :> (RowBox[{"PointLegend", "[", RowBox[{ RowBox[{"{", RowBox[{ RowBox[{"Directive", "[", RowBox[{"RGBColor", "[", RowBox[{"0.24720000000000014`", ",", "0.24`", ",", "0.6`"}], "]"}], "]"}], ",", RowBox[{"Directive", "[", RowBox[{"RGBColor", "[", RowBox[{"0.6`", ",", "0.24`", ",", "0.4428931686004542`"}], "]"}], "]"}]}], "}"}], ",", RowBox[{"{", RowBox[{#, ",", #2}], "}"}], ",", RowBox[{"LegendLayout", "\[Rule]", "\"Column\""}], ",", RowBox[{"LegendMarkers", "\[Rule]", "False"}]}], "]"}]& )]}, "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.591949207924857*^9, 3.591949500156979*^9}] }, Open ]] }, Open ]] }, Open ]] }, Open ]] }, Open ]] }, WindowSize->{1312, 625}, WindowMargins->{{0, Automatic}, {Automatic, 0}}, Magnification->1.5000001192092896`, FrontEndVersion->"9.0 for Microsoft Windows (64-bit) (January 25, 2013)", StyleDefinitions->"Default.nb" ] (* End of Notebook Content *) (* Internal cache information *) (*CellTagsOutline CellTagsIndex->{} *) (*CellTagsIndex CellTagsIndex->{} *) (*NotebookFileOutline Notebook[{ Cell[CellGroupData[{ Cell[579, 22, 167, 2, 56, "Subsubtitle"], Cell[CellGroupData[{ Cell[771, 28, 152, 3, 117, "Section"], Cell[CellGroupData[{ Cell[948, 35, 253, 6, 45, "Input"], Cell[1204, 43, 279, 8, 76, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[1520, 56, 714, 16, 104, "Input"], Cell[2237, 74, 94, 1, 44, "Output"] }, Open ]], Cell[2346, 78, 196, 4, 45, "Text"], Cell[CellGroupData[{ Cell[2567, 86, 314, 8, 45, "Input"], Cell[2884, 96, 628, 21, 82, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[3549, 122, 213, 4, 45, "Input"], Cell[3765, 128, 512, 16, 85, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[4314, 149, 328, 5, 86, "Subsubsection"], Cell[CellGroupData[{ Cell[4667, 158, 624, 17, 45, "Input"], Cell[5294, 177, 292, 5, 65, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[5635, 188, 161, 3, 65, "Subsection"], Cell[CellGroupData[{ Cell[5821, 195, 269, 5, 51, "Subsubsection"], Cell[6093, 202, 112, 1, 45, "Text"], Cell[CellGroupData[{ Cell[6230, 207, 304, 7, 45, "Input"], Cell[6537, 216, 516, 10, 73, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[7090, 231, 127, 2, 45, "Input"], Cell[7220, 235, 127, 2, 44, "Output"] }, Open ]], Cell[7362, 240, 213, 4, 45, "Text"], Cell[7578, 246, 564, 13, 45, "Input"], Cell[CellGroupData[{ Cell[8167, 263, 129, 2, 45, "Input"], Cell[8299, 267, 145, 2, 44, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[8493, 275, 167, 3, 51, "Subsubsection"], Cell[CellGroupData[{ Cell[8685, 282, 674, 18, 104, "Input"], Cell[9362, 302, 2360, 60, 157, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[11759, 367, 543, 15, 45, "Input"], Cell[12305, 384, 3102, 57, 699, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[15444, 446, 499, 15, 45, "Input"], Cell[15946, 463, 3016, 56, 378, "Output"] }, Open ]] }, Open ]], Cell[CellGroupData[{ Cell[19011, 525, 173, 3, 51, "Subsubsection"], Cell[19187, 530, 324, 6, 45, "Text"], Cell[19514, 538, 365, 10, 45, "Input"], Cell[19882, 550, 305, 6, 45, "Text"], Cell[CellGroupData[{ Cell[20212, 560, 345, 10, 45, "Input"], Cell[20560, 572, 131, 2, 44, "Output"] }, Open ]], Cell[20706, 577, 299, 6, 45, "Text"], Cell[CellGroupData[{ Cell[21030, 587, 500, 15, 75, "Input"], Cell[21533, 604, 158, 3, 44, "Output"] }, Open ]], Cell[21706, 610, 165, 3, 45, "Text"], Cell[CellGroupData[{ Cell[21896, 617, 706, 19, 104, "Input"], Cell[22605, 638, 143, 2, 44, "Output"] }, Open ]], Cell[22763, 643, 156, 3, 45, "Text"], Cell[CellGroupData[{ Cell[22944, 650, 1100, 29, 191, "Input"], Cell[24047, 681, 1854, 45, 157, "Output"] }, Open ]], Cell[25916, 729, 135, 3, 45, "Text"], Cell[CellGroupData[{ Cell[26076, 736, 672, 18, 104, "Input"], Cell[26751, 756, 1787, 44, 157, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[28575, 805, 466, 9, 45, "Input"], Cell[29044, 816, 1498, 30, 405, "Output"] }, Open ]], Cell[CellGroupData[{ Cell[30579, 851, 501, 13, 75, "Input"], Cell[31083, 866, 5613, 123, 369, "Output"] }, Open ]], Cell[36711, 992, 224, 4, 45, "Text"], Cell[CellGroupData[{ Cell[36960, 1000, 339, 10, 45, "Input"], Cell[37302, 1012, 7709, 158, 401, "Output"] }, Open ]] }, Open ]] }, Open ]] }, Open ]] }, Open ]] } ] *) (* End of internal cache information *)