2f:T566,trReport = Plot.plot({ width: 400, height: 300, y: {doman: [0,200], axis:null}, x: {domain: [0,200], axis:null}, marks: [ Plot.text([[5,190]], {text: ["Base Rate: "+ (rate.fn + rate.tp).toFixed(2) ], fontSize: 16, textAnchor: "start", dy:10}), Plot.text([[5, 170]], {text: ["Select Rate: " + (rate.tp + rate.fp).toFixed(2)], fontSize: 16, textAnchor: "start", dy:10}), Plot.text([[5, 150]], {text:["False Positive: " + (rate.fp/(rate.fp+rate.tn)).toFixed(2)], fontSize: 16, textAnchor: "start", dy:10}), Plot.text([[5, 130]], {text:["False Negative: " + (1 - rate.tp/(rate.tp+rate.fn)).toFixed(2)], fontSize: 16, textAnchor: "start", dy:10}), Plot.text([[5, 110]], {text:["PPV: " + (rate.tp/(rate.tp+rate.fp)).toFixed(2)], fontSize: 16, textAnchor: "start", dy:10}), Plot.text([[5, 90]], {text:["NPV: " + (rate.tn/(rate.tn+rate.fn)).toFixed(2)], fontSize: 16, textAnchor: "start", dy:10}), Plot.text([[5, 70]], {text:["TPR 'Hit Rate': " + (rate.tp/(rate.tp+rate.fn)).toFixed(2)], fontSize: 16, textAnchor: "start", dy:10}), Plot.text([[5, 50]], {text:["FPR 'False Alarm Rate': " + (rate.fp/(rate.fp+rate.tn)).toFixed(2)], fontSize: 16, textAnchor: "start", dy:10}), Plot.text([[5, 30]], {text:["Accuracy: " + (rate.tp+rate.tn).toFixed(2)], fontSize: 16, textAnchor: "start", dy:10}), Plot.text([[0,0]], {text: [" "]}) ] })30:[] 2d:{"title":"Taylor-Russell Diagram and Computations","cells":[{"id":0,"value":"# Taylor-Russell Diagram and Computations","pinned":false,"hidden":false,"mode":"md"},{"id":38,"value":"viewof rho = Inputs.range([-1,1], {label: \"Rho\", step: 0.01, start: 0})","pinned":false,"hidden":false,"mode":"ojs"},{"id":49,"value":"viewof cutoff = Inputs.range([-3,3], {label: \"Cutoff\", step: 0.01, start: 0})","pinned":false,"hidden":false,"mode":"ojs"},{"id":55,"value":"viewof truecut = Inputs.range([-3,3], {label: \"truecut\", step: 0.01, start: 0})","pinned":false,"hidden":false,"mode":"ojs"},{"id":214,"value":"\n \n \n \n

\n ${trPlot}\n

\n ${trReport}\n

\n \n","pinned":false,"hidden":false,"mode":"html"},{"id":36,"value":"trPlot = Plot.plot({\n x: {domain: [-3,3]},\n y: {domain: [-3,3]},\n marks: [\n Plot.dot(data, {x: \"z1\", y: \"z2\", r: 1}),\n Plot.ruleX([cutoff], {strokeWidth: 3, stroke: 'blue'}),\n Plot.ruleY([truecut], {strokeWidth: 3, stroke: 'blue'} )\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":161,"value":"$2f","pinned":false,"hidden":false,"mode":"ojs"},{"id":135,"value":"| | | Decision | | |\n|-------|----------|------------|----------|---|\n| | | Negative | Positive | |\n| **Truth** | Positive | ${rate.fn.toFixed(2)} | ${rate.tp.toFixed(2)} | |\n| | Negative | ${rate.tn.toFixed(2)} | ${rate.fp.toFixed(2)} | |","pinned":false,"hidden":false,"mode":"md"},{"id":63,"value":"## Setup","pinned":false,"hidden":false,"mode":"md"},{"id":21,"value":"mu1 = 0","pinned":false,"hidden":false,"mode":"ojs"},{"id":23,"value":"mu2 = 0","pinned":false,"hidden":false,"mode":"ojs"},{"id":8,"value":"s1=1","pinned":false,"hidden":false,"mode":"ojs"},{"id":10,"value":"s2 = 1","pinned":false,"hidden":false,"mode":"ojs"},{"id":12,"value":"lambda = s2/s1*rho","pinned":false,"hidden":false,"mode":"ojs"},{"id":17,"value":"nu = Math.sqrt(1-rho*rho)*s1*s2","pinned":false,"hidden":false,"mode":"ojs"},{"id":19,"value":"function biRand() {let z1 = mu1 + s1 * jStat.normal.sample(0,1);\n let z2 = mu2 + lambda*(z1 - mu1) + nu* jStat.normal.sample(0,1);\n return {z1,z2}\n }","pinned":false,"hidden":false,"mode":"ojs"},{"id":28,"value":"n = 10000","pinned":false,"hidden":false,"mode":"ojs"},{"id":30,"value":"data = {let dd = [];\n for(let i = 1; i <= n; i++)\n dd.push(biRand())\n return dd;\n }\n ","pinned":false,"hidden":false,"mode":"ojs"},{"id":74,"value":"rate = {let fp = 0; let tp = 0; let fn = 0; let tn = 0;\n for(let i = 0; i <= data.length-1; i++)\n {let x1 = data[i].z1;\n let x2 = data[i].z2;\n if(x1 >= cutoff && x2 >= truecut) tp++\n else if(x1 >= cutoff && x2 < truecut) fp++\n else if(x1 < cutoff && x2 >= truecut) fn++\n else tn++\n }\n return({fp: fp/n, tp: tp/n, fn: fn/n, tn: tn/n})}","pinned":false,"hidden":false,"mode":"ojs"},{"id":6,"value":"jStat = require('jStat')","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W30","stdlib":"1","titling":"auto"}