28:null 29:["$","title",null,{"children":"Taylor-Russell Diagram and Computations for Two Groups | Gary McClelland | Observable"}] 2e: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]) ] })2f:[] 2c:{"id":"c6a01374acb3791e","slug":null,"title":"Taylor-Russell Diagram and Computations for Two Groups","thumbnail":"f56196b811cdb964e138864515dc0a202efc11ecbb44fa13f2ffdf0b9d2b78f3","access_level":"public","create_time":"$D2023-12-14T05:17:11.258Z","update_time":"$D2023-12-22T18:11:26.496Z","view_time":null,"first_public_time":"$D2023-12-14T05:17:11.258Z","has_unpublished_edits":null,"fork":{"id":"62644f3923f087e1","slug":null,"title":"Taylor-Russell Diagram and Computations","owner":{"id":"8e902c6c168a5327","login":"mcclella","name":"Gary McClelland"}},"stars":1,"starred":false,"owner":{"type":"team","id":"8e902c6c168a5327","login":"mcclella","name":"Gary McClelland","avatar_url":"https://avatars-next.observableusercontent.com/5fb0f9ce97ef2ff380c45d797b84b995cfeccb679fa1fcc25ae40f1a1b81414b","tier":"starter_2024"},"version":404,"publish_version":null,"role":null,"authored":false,"authors":[{"type":"individual","id":"bb02701f3bca15dc","login":"mcclella","name":"Gary McClelland","avatar_url":"https://avatars-next.observableusercontent.com/5fb0f9ce97ef2ff380c45d797b84b995cfeccb679fa1fcc25ae40f1a1b81414b","tier":null}],"members":[{"type":"individual","id":"bb02701f3bca15dc","login":"mcclella","name":"Gary McClelland","avatar_url":"https://avatars-next.observableusercontent.com/5fb0f9ce97ef2ff380c45d797b84b995cfeccb679fa1fcc25ae40f1a1b81414b","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":0,"collection":null,"type":"notebook","body":{"titling":"auto","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","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":38,"value":"viewof rho = Inputs.range([0,1], {label: \"Rho\", step: 0.01, start: 0})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":55,"value":"viewof truecut = Inputs.range([-3,3], {label: \"truecut\", step: 0.01, start: 0})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":383,"value":"viewof mu2 = Inputs.range([0,1.5], {label: \"Relative Prevalence\", step: 0.01, start: 0})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":214,"value":"\n \n \n \n

\n ${trPlot}\n

\n ${trReport}\n

\n \n","pinned":false,"hidden":false,"mode":"html","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":161,"value":"$2e","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":63,"value":"## Setup","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":8,"value":"s1=1","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":10,"value":"s2 = 1","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":12,"value":"lambda = s2/s1*rho","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":17,"value":"nu = Math.sqrt(1-rho*rho)*s1*s2","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":28,"value":"n = 10000","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":277,"value":"dataAll = data1.concat(data2)","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":6,"value":"jStat = require('jStat')","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":237,"value":"graphColors = ({\n blue: \"#4e79a7\",\n orange: \"#f2832d\",\n red: \"#e15759\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W2f","stdlib":"1"}} 2d:"$2c:body"