31:[] 2f:{"title":"Cauchy distribution","cells":[{"id":38,"value":"# Cauchy distribution","pinned":false,"hidden":false,"mode":"md"},{"id":523,"value":"

${tex`X \\sim \\mathrm{Cauchy}(\\mu,\\tau^2)`}, for ${tex`-\\infty \\leq X \\leq \\infty`}","pinned":false,"hidden":false,"mode":"md"},{"id":406,"value":"The density of a ${tex`X \\sim \\mathrm{Cauchy}(\\mu,\\tau^2)`} random variable is\n\n${tex`f(x) = \\frac{1}{\\pi\\tau\\Big(1+\\big( \\frac{x-\\mu}{\\tau} \\big)^2 \\Big)}`}","pinned":false,"hidden":false,"mode":"md"},{"id":408,"value":"moments = tex`\n\\begin{aligned}\n&E( X) = \\text{undefined} \\\\[1em]\n&Var( X) = \\text{undefined}\n\\end{aligned}\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":118,"value":"viewof myinputs = Inputs.form([\n Inputs.range([-5, 5], {value: 0, step: 0.1, label: tex`\\text{location }\\mu`}),\n Inputs.range([0.1, 5], {value: 1, step: 0.1, label: tex`\\text{scale }\\tau`}),\n Inputs.range([-10, 10], {value: -1.96, step: 0.01, label: \"quantile\"})\n ])","pinned":false,"hidden":false,"mode":"ojs"},{"id":450,"value":"viewof shownormal = Inputs.toggle({value: false, label: \"show normal\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":462,"value":"tailprobs = \n{\n var tailprobs;\n if (shownormal){\n tailprobs = tex`\\text{Normal distribution: } P(X\\leq ${myinputs[2]})=${normcdf.toPrecision(4)} \\\\ \\text{Cauchy distribution: } P(X\\leq ${myinputs[2]})=${studentcdf.toPrecision(4)}`\n }\n else{\n tailprobs = tex`\\text{Cauchy distribution: } P(X\\leq ${myinputs[2]})=${studentcdf.toPrecision(4)}`\n }\n return tailprobs;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":305,"value":"plt = Plot.plot({\n color: {\n legend: shownormal\n },\n x: {\n label: \"x\",\n axis: true\n },\n y: {\n axis: false,\n domain: [0,1.05*jstat.normal.pdf(myinputs[0], myinputs[0], myinputs[1])]\n },\n marks: [\n Plot.ruleY([0]),\n Plot.line(data, {filter: d => d.dist == \"student-t\", x: \"x\", y: \"pdf\", stroke : \"dist\", strokeWidth: 2}),\n Plot.line(data, {filter: d => d.dist == \"normal\" && shownormal, x: \"x\", y: \"pdf\", stroke : \"dist\", strokeWidth: 2}),\n Plot.areaY(data, {filter: d => d.x <= myinputs[2] && d.dist == \"normal\" && shownormal, x: \"x\", y: \"pdf\", fill: \"steelblue\", opacity: 0.2}),\n Plot.areaY(data, {filter: d => d.x <= myinputs[2] && d.dist == \"student-t\", x: \"x\", y: \"pdf\", fill: \"orange\", opacity: 0.2}),\n \n ]\n })","pinned":false,"hidden":false,"mode":"ojs"},{"id":515,"value":"## Relations to other distributions\n\n- If ${tex`X\\sim N(0,1)`} and ${tex`Y\\sim N(0,1)`} are independent, then ${tex`X/Y\\sim\\mathrm{Cauchy}(0,1)`}\n- If ${tex`X \\sim \\mathrm{Cauchy}(0,1)`} then ${tex`\\mu + \\tau X \\sim \\mathrm{Cauchy}(\\mu,\\tau^2)`}","pinned":false,"hidden":false,"mode":"md"},{"id":109,"value":"
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n
\n\n---\n## 🙈 Appendix 🙈\n\nHere is only code and details that you probably don't want to read.","pinned":false,"hidden":false,"mode":"md"},{"id":221,"value":"#### Imports","pinned":false,"hidden":false,"mode":"md"},{"id":26,"value":"jstat = require('jstat')","pinned":false,"hidden":false,"mode":"ojs"},{"id":223,"value":"#### Data","pinned":false,"hidden":false,"mode":"md"},{"id":17,"value":"data = {\n const x = d3.range(-10, 10, 0.01);\n const normpdf = x.map(x => ({x: x, pdf: jstat.normal.pdf(x, myinputs[0], myinputs[1]), dist: \"normal\"}));\n const studentpdf = x.map(x => ({x: x, pdf: jstat.studentt.pdf((x-myinputs[0])/myinputs[1], 1)/myinputs[1], dist: \"student-t\"}));\n const data = normpdf.concat(studentpdf)\n return data\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":19,"value":"normcdf = jstat.normal.cdf(myinputs[2], myinputs[0], myinputs[1]);","pinned":false,"hidden":false,"mode":"ojs"},{"id":23,"value":"studentcdf = jstat.studentt.cdf((myinputs[2]-myinputs[0])/myinputs[1], 1);","pinned":false,"hidden":false,"mode":"ojs"},{"id":229,"value":"#### Style","pinned":false,"hidden":false,"mode":"md"},{"id":202,"value":"","pinned":false,"hidden":false,"mode":"html"},{"id":58,"value":"viewof md_colour = Inputs.color({label: \"Highlight color cells\", value: \"#F5F5DC\"})","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W31","stdlib":"1","titling":"auto"}