2f:T415,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: "+ (rateAll.fn + rateAll.tp).toFixed(2) ], fontSize: 16, textAnchor: "start", dy:10}), Plot.text([[5, 170]], {text: ["Select Rate: " + (rateAll.tp + rateAll.fp).toFixed(2)], fontSize: 16, textAnchor: "start", dy:10}), Plot.text([[5,140]], {text: ["Base Rate: "+ (rate1.fn + rate1.tp).toFixed(2) ], fontSize: 16, textAnchor: "start", dy:10, fill: graphColors.blue}), Plot.text([[5, 120]], {text: ["Select Rate: " + (rate1.tp + rate1.fp).toFixed(2)], fontSize: 16, textAnchor: "start", dy:10, fill: graphColors.blue}), Plot.text([[5,90]], {text: ["Base Rate: "+ (rate2.fn + rate2.tp).toFixed(2) ], fontSize: 16, textAnchor: "start", dy:10, fill: graphColors.orange}), Plot.text([[5, 70]], {text: ["Select Rate: " + (rate2.tp + rate2.fp).toFixed(2)], fontSize: 16, textAnchor: "start", dy:10, fill: graphColors.orange}), Plot.ruleY([0]) ] })30:[] 2d:{"title":"Taylor-Russell Diagram and Computations for Two Groups","cells":[{"id":0,"value":"# Taylor-Russell Diagram and Computations for Two Groups","pinned":false,"hidden":false,"mode":"md"},{"id":38,"value":"viewof rho = Inputs.range([0,1], {label: \"Rho\", 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":383,"value":"viewof mu2 = Inputs.range([0,1.5], {label: \"Relative Prevalence\", 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,4]},\n marks: [\n Plot.dot(data1, {x: \"z1\", y: \"z2\", r: 1, fill: graphColors.blue, opacity: 0.5}),\n Plot.dot(data2, {x: \"z1\", y: \"z2\", r: 1, fill: graphColors.orange, opacity: 0.5}),\n Plot.ruleX([cutoff], {strokeWidth: 3, stroke: graphColors.red}),\n Plot.ruleY([truecut], {strokeWidth: 3, stroke: graphColors.red} )\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 | ${rateAll.fn.toFixed(2)} | ${rateAll.tp.toFixed(2)} | |\n| | Negative | ${rateAll.tn.toFixed(2)} | ${rateAll.fp.toFixed(2)} | |","pinned":false,"hidden":false,"mode":"md"},{"id":63,"value":"## Setup","pinned":false,"hidden":false,"mode":"md"},{"id":350,"value":"br =[ {all: (rateAll.fn + rateAll.tp).toFixed(3), g1: (rate1.fn + rate1.tp).toFixed(3), g2: (rate2.fn + rate2.tp).toFixed(2) }]","pinned":false,"hidden":false,"mode":"ojs"},{"id":353,"value":"sr = [{all: (rateAll.tp + rateAll.fp).toFixed(3), g1: (rate1.tp + rate1.fp).toFixed(3), g2: (rate2.tp + rate2.fp).toFixed(3)}]","pinned":false,"hidden":false,"mode":"ojs"},{"id":230,"value":"cutoff = -jStat.normal.inv(.5*(1-jStat.normal.cdf(truecut,0,1) + 1 - jStat.normal.cdf(truecut,mu2,1)),0,1)","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(m1, m2) {let z1 = m1 + s1 * jStat.normal.sample(0,1);\n let z2 = m2 + lambda*(z1 - m1) + 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":"data1 = {let dd = [];\n for(let i = 1; i <= n; i++)\n dd.push(biRand(0,0))\n return dd;\n }\n ","pinned":false,"hidden":false,"mode":"ojs"},{"id":246,"value":"data2 = {let dd = [];\n for(let i = 1; i <= n; i++)\n dd.push(biRand(0,mu2))\n return dd;\n }","pinned":false,"hidden":false,"mode":"ojs"},{"id":277,"value":"dataAll = data1.concat(data2)","pinned":true,"hidden":false,"mode":"ojs"},{"id":74,"value":"rate1 = {let fp = 0; let tp = 0; let fn = 0; let tn = 0;\n for(let i = 0; i <= data1.length-1; i++)\n {let x1 = data1[i].z1;\n let x2 = data1[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":266,"value":"rate2 = {let fp = 0; let tp = 0; let fn = 0; let tn = 0;\n for(let i = 0; i <= data2.length-1; i++)\n {let x1 = data2[i].z1;\n let x2 = data2[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":279,"value":"rateAll = {let fp = 0; let tp = 0; let fn = 0; let tn = 0;\n for(let i = 0; i <= dataAll.length-1; i++)\n {let x1 = dataAll[i].z1;\n let x2 = dataAll[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/(2*n), tp: tp/(2*n), fn: fn/(2*n), tn: tn/(2*n)})}","pinned":false,"hidden":false,"mode":"ojs"},{"id":6,"value":"jStat = require('jStat')","pinned":false,"hidden":false,"mode":"ojs"},{"id":237,"value":"graphColors = ({\n blue: \"#4e79a7\",\n orange: \"#f2832d\",\n red: \"#e15759\"\n})","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W30","stdlib":"1","titling":"auto"}