28:null 29:["$","title",null,{"children":"Taylor-Russell Diagram and Computations | Gary McClelland | Observable"}] 2e: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: [" "]}) ] })2f:[] 2c:{"id":"62644f3923f087e1","slug":null,"title":"Taylor-Russell Diagram and Computations","thumbnail":"03cae96c008416827435332776c516eacb855ebcb595305786d7773657161c32","access_level":"public","create_time":"$D2023-12-08T23:23:38.295Z","update_time":"$D2023-12-11T23:09:09.891Z","view_time":null,"first_public_time":"$D2023-12-08T23:23:38.295Z","has_unpublished_edits":null,"fork":{"id":"75a03d47c3d920de","slug":null,"title":"Untitled","owner":{"id":"8e902c6c168a5327","login":"mcclella","name":"Gary McClelland"}},"stars":0,"starred":false,"owner":{"type":"team","id":"8e902c6c168a5327","login":"mcclella","name":"Gary McClelland","avatar_url":"https://avatars-next.observableusercontent.com/5fb0f9ce97ef2ff380c45d797b84b995cfeccb679fa1fcc25ae40f1a1b81414b","tier":"starter_2024"},"version":227,"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":1,"collection":null,"type":"notebook","body":{"titling":"auto","title":"Taylor-Russell Diagram and Computations","cells":[{"id":0,"value":"# Taylor-Russell Diagram and Computations","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"","database":"$undefined"},{"id":38,"value":"viewof rho = Inputs.range([-1,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":49,"value":"viewof cutoff = Inputs.range([-3,3], {label: \"Cutoff\", 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":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,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","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 | ${rate.fn.toFixed(2)} | ${rate.tp.toFixed(2)} | |\n| | Negative | ${rate.tn.toFixed(2)} | ${rate.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":21,"value":"mu1 = 0","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"id":23,"value":"mu2 = 0","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() {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","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":"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","author":"$undefined","summary":"$undefined","providerMetadata":"$undefined","attachments":"$undefined","output":"$undefined","database":"$undefined"},{"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","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"}],"files":"$W2f","stdlib":"1"}} 2d:"$2c:body"