31:T44d,Plot.plot({ inset: 50, x: {label: "Std BUSTED", type : "log"}, y: {label: "BUSTED-E", type : "log"}, color: {legend: true}, marks: [ Plot.dot (_.filter (data, (d)=>d.kind != '--'), {"x" : d => Math.min (d["omega3_BUSTED"], 1000), "y" : d => Math.max (d["p"], 1e-6), fill : "kind"}), //Plot.dot (_.filter (data, (d)=>d["omega3_BUSTED"] < 100 && d["omega3_BUSTED-E"] < 100 && d["E_better"] <= 0.05 && (d["kind"] == "+-" || d["kind"] == "++")), {x: "omega3_BUSTED", y: "omega3_BUSTED-E", fill : "kind", "opacity" : 0.5}), //Plot.density(_.filter (data, (d)=>d["omega3_BUSTED"] < 100 && d["omega3_BUSTED-E"] < 100 && d["E_better"] <= 0.05 && d["kind"] == "++"), {x: "omega3_BUSTED", y: "omega3_BUSTED-E", stroke: "darkblue", thresholds: 10, strokeWidth: 2, stroke : "kind", "opacity" : 0.8}), //Plot.density(_.filter (data, (d)=>d["omega3_BUSTED"] < 100 && d["omega3_BUSTED-E"] < 100 && d["E_better"] <= 0.05 && d["kind"] == "+-"), {x: "omega3_BUSTED", y: "omega3_BUSTED-E", stroke: "darkblue", thresholds: 10, strokeWidth: 2, stroke : "kind", "opacity" : 0.8}), ] })32:T548,Plot.plot({ style: {"fontSize" : 14, "fontSize-swatches" : 14, "fontFamily" : "Helvetica"}, title : md`ω₃ **magnitude** by selection detection type`, round: true, color: {legend: true, "domain" : ["BUSTED","BUSTED-E"], "range" : ["midnightblue", "darkorange"]}, //color: {legend: true}, height: 800, width : 800, marginTop: 40, marginBottom: 40, marginLeft: 50, paddingBottom : 10, facet: { data : plot_data, x: 'BUSTED', y: 'BUSTED-E', }, x: {type: "linear", ticks : 5, 'label' : 'log10 (ω3)'}, y: {type: "linear", ticks : 5, 'label' : 'log10 (p[ω3])'}, //x: {}, marks: [ //Plot.rectY(plot_data, //Plot.binX({y2: "proportion-facet"}, {x: "omega", fill: "model", inset: 0, thresholds: 100, mixBlendMode: "multiply"})), Plot.density(plot_data, {x: "omega", y: "weight", fill: "model", "strokeWidth" : 0.5, thresholds: 100, bandwidth: 1,mixBlendMode: "luminosity"}), //Plot.rect(_.filter (plot_data, (d)=>d.model == "BUSTED"), Plot.bin({fill: "count"}, {x: "omega", y: "weight", thresholds: 50, tip: true, })), Plot.frame({}), Plot.axisY({anchor: "left", labelAnchor: "top", ticks : 5}), Plot.rect([{'x1' : 2, 'y1': -2, 'x2' : 5.0, 'y2': -5}], {'x1' : 'x1', 'y1': 'y1', 'x2' : 'x2', 'y2': 'y2', opacity : 0.2, "fill" : "black"}), //Plot.ruleY([0]) ] })33:T571,{ let dd = _.filter (plot_data, (d)=>d.model == plot_model); return Plot.plot({ style: {"fontSize" : 14, "fontFamily" : "Helvetica"}, //title : md`ω₃ **magnitude** by selection detection type`, round: true, //color: {legend: true, "domain" : ["BUSTED","BUSTED-E"], "range" : ["red", "blue"]}, color: {legend: true, type : "log", width : 500, ticks: 5, scheme: colors, reverse: true}, height: 500, width : 500, marginTop: 40, marginBottom: 40, marginLeft: 50, paddingBottom : 10, facet: { data : dd, x: 'BUSTED', y: 'BUSTED-E', }, x: {type: "linear", ticks : 5, 'label' : 'log10 (ω3)'}, y: {type: "linear", ticks : 5, 'label' : 'log10 (p[ω3])'}, //x: {}, marks: [ //Plot.rectY(plot_data, //Plot.binX({y2: "proportion-facet"}, {x: "omega", fill: "model", inset: 0, thresholds: 100, mixBlendMode: "multiply"})), //Plot.density(plot_data, {x: "omega", y: "weight", fill: "model", "strokeWidth" : 0.5, thresholds: 100, bandwidth: 0.01,mixBlendMode: "color"}), Plot.rect(dd, Plot.bin({fill: "count"}, {x: "omega", y: "weight", interval : 1/6, inset: 0.5, tip: true, })), Plot.frame({}), Plot.axisY({anchor: "left", labelAnchor: "top", ticks : 5}), Plot.rect([{'x1' : 2, 'y1': -2, 'x2' : 5.5, 'y2': -5.5}], {'x1' : 'x1', 'y1': 'y1', 'x2' : 'x2', 'y2': 'y2', opacity : 0.2, "fill" : "black"}), ] })}34:T63e,data = { let gene_handler = null; let annotation = null; if (urls[dataset].genes) { if (_.isFunction (urls[dataset].genes)) { gene_handler = urls[dataset].genes; } else { annotation = _.chain (await obtainSourceData(null, {'url' : urls[dataset].genes, 'tsv' : true, 'gzip' : true})).map ((d)=>[d.ID,d.NAME]).fromPairs ().value(); } } return _.map (await d3.csv (urls[dataset].address,d3.autoType), (d)=> { d.File = "" + d.File; if (annotation) { d.gene = annotation[d.File.split ('.')[0]] || d.File; } else { d.gene = gene_handler ? gene_handler (d) : d.File.split ('.')[0]; } d.kind = (d.p <= 0.05 ? '+' : '-') + (d.p_e <= 0.05 ? '+' : '-'); //_.each (d, (f,k)=> { //d.omega3_BUSTED = Math.min (1000, d.omega3_BUSTED); //d.omega_e = Math.min (1000, d.omega_e); //d["omega3_BUSTED-E"] = Math.min (1000, d["omega3_BUSTED-E"]); d.T = Math.min (d.T,100); const deltaaic = d["BUSTED_E_AIC"]; let weights = [0,0]; if (deltaaic < 0) { // +E is better weights[0] = Math.exp (deltaaic/2); weights[1] = 1; } else { weights[1] = Math.exp (-deltaaic/2); weights[0] = 1; } d.p_MA = (d.p*weights[0] + d.p_e*weights[1])/(weights[0]+weights[1]); d.E_better = 1-chi2cdf.default (2*(d.logL_e-d.logL), 2) ; //}); d.kind = (d.p <= 0.05 ? '+' : '-') + (d.p_e <= 0.05 ? '+' : '-') + (d.p_MA <= 0.05 ? '+' : '-'); return d; }); }35:[] 2f:{"title":"BUSTED vs BUSTED-E","cells":[{"id":4,"value":"# BUSTED vs BUSTED-E\n\nExploring the effect of adding an error-sink component on BUSTED selection detection","pinned":false,"hidden":false,"mode":"md"},{"id":20,"value":"viewof dataset = Inputs.select (_.keys (urls), {label : \"Dataset\", value : \"PMID 20333182\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":6,"value":"viewof search = Inputs.search(data, {autocomplete : true})","pinned":false,"hidden":false,"mode":"ojs"},{"id":7,"value":"Inputs.table(search, {\n rows: 20,\n format: {\n p : (d,i)=>annotator (d,i, false),\n p_e : (d,i)=>annotator (d,i, true),\n prop_e : (d)=>d3.format (\"4p\")(d),\n gene: (d)=>html`${d}`\n }\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":40,"value":"Inputs.table(get2x2 (search))","pinned":false,"hidden":false,"mode":"ojs"},{"id":930,"value":"Plot.plot({\n style: {\"fontSize\" : 14, \"fontSize-swatches\" : 14, \"fontFamily\" : \"Helvetica\"},\n marginBottom: 40,\n marginTop: 30,\n marginLeft : 50,\n x: {label: \"Significance cutoff\", type : \"log\", tickFormat : \"g\"},\n y: {label: \"Detection rate\", type : \"linear\", zero: false, tickFormat : \"g\"},\n color: {legend: true, domain: [\"BUSTED\", \"BUSTED-E\", \"MA\"], range : [\"firebrick\", \"steelblue\", \"forestgreen\"]},\n marks: [\n Plot.dot (p_trends, {\"x\" : \"p\", \"y\" : \"detection\", \"fill\" : \"method\", \"r\" : 5, tip : true}),\n Plot.line (p_trends, {\"x\" : \"p\", \"y\" : \"detection\", \"stroke\" : \"method\"}),\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":919,"value":"p_trends = _.flatten(_.map ([0.1,0.05,0.02,0.01,0.005,0.001,0.0005,0.0001,1e-5,1e-6,1e-8], (d)=> {\n let tt = get2x2 (search, d);\n let tt2 = get2x2 (search, d, \"p_e\",\"p_MA\");\n return [{'p' : d, 'method' : 'BUSTED', 'detection' : tt[4].p/100},{'p' : d, 'method' : 'BUSTED-E', 'detection' : tt[5].p/100},{'p' : d, 'method' : 'MA', 'detection' : tt2[5].p/100}];\n}))","pinned":false,"hidden":false,"mode":"ojs"},{"id":76,"value":"viewof var1 = Inputs.select (_.keys (data[0]), {label : \"var1\", value : \"omega3_BUSTED\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":79,"value":"viewof var2 = Inputs.select (_.keys (data[0]), {label : \"var2\" , value : \"omega3_BUSTED-E\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":830,"value":"//_.sortBy (_.filter (_.map (data, (d)=>[d.File,d.logL-d.logL_e,d.logL,d.logL_e,d.p,d.p_e]), (d)=>d[1]>0.25), (d)=>-d[1]), (d)=>d[0] + \".*.gz\").join (\" \")","pinned":false,"hidden":false,"mode":"ojs"},{"id":886,"value":"_.map (_.filter (data, (d)=>(d.logL - d.logL_e >= 0.5) || (d.p >= 0.05 && d.p_e < 0.05)), (d)=>d.File + \".*.gz\").join (\" \")","pinned":false,"hidden":false,"mode":"ojs"},{"id":361,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":482,"value":"Plot.plot({\n inset: 50,\n x: {label: \"Std BUSTED\", type : \"log\"},\n y: {label: \"BUSTED-E\", type : \"linear\"},\n color: {legend: true},\n marks: [\n Plot.link(_.filter (data, (d)=>d.kind == \"++\" && d.BUSTED_E_AIC < -5), {\n x1: d=>Math.min (10, d[\"omega3_BUSTED\"]),\n y1: \"p3_BUSTED\",\n x2: d=>Math.min (10, d[\"omega3_BUSTED-E\"]),\n y2: \"p3_BUSTED-E\",\n //stroke: (d) => d.R90_10_2015 - d.R90_10_1980,\n markerEnd: \"arrow\"\n }),\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":64,"value":"pairplot(search, {'vars' : [var1, var2], kind: 'scatter', 'hue' : 'kind', diag_kind: \"hist\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":545,"value":"plot_data = _.flatten(_.map (data, (d)=> {\n return [\n {'omega' : Math.min (5, Math.log10 (d[\"omega3_BUSTED\"])), 'weight' : Math.max (-5, Math.log10 (Math.max (1e-6,d[\"p3_BUSTED\"]))), 'class' : d[\"kind\"], 'model' : 'BUSTED', 'BUSTED' : d[\"kind\"][0] == '-' ? \"No\" : \"Yes\", 'BUSTED-E' : d[\"kind\"][1] == '-' ? \"No\" : \"Yes\"},\n {'omega' : Math.min (5, Math.log10 (d[\"omega3_BUSTED-E\"])), 'weight' : Math.max (-5, Math.log10 (Math.max (1e-6,d[\"p3_BUSTED-E\"]))), 'class' : d[\"kind\"], 'model' : 'BUSTED-E', 'BUSTED' : d[\"kind\"][0] == '-' ? \"No\" : \"Yes\", 'BUSTED-E' : d[\"kind\"][1] == '-' ? \"No\" : \"Yes\"} \n ]\n })) ","pinned":false,"hidden":false,"mode":"ojs"},{"id":1167,"value":"d3.max (plot_data, (d)=>d.weight)","pinned":true,"hidden":false,"mode":"ojs"},{"id":619,"value":"plot_data_w = _.flatten(_.map (data, (d)=> {\n return [\n {'omega' : Math.min (5, Math.log10 (Math.max (1e-6,d[\"p3_BUSTED\"]))), 'model' : 'BUSTED', 'BUSTED' : d[\"kind\"][0] == '-' ? \"No\" : \"Yes\", 'BUSTED-E' : d[\"kind\"][1] == '-' ? \"No\" : \"Yes\", 'class' : d[\"kind\"]},\n {'omega' : Math.min (5, Math.log10 (Math.max (1e-6,d[\"p3_BUSTED-E\"]))), 'model' : 'BUSTED-E', 'BUSTED' : d[\"kind\"][0] == '-' ? \"No\" : \"Yes\", 'BUSTED-E' : d[\"kind\"][1] == '-' ? \"No\" : \"Yes\", 'class' : d[\"kind\"]} \n ]\n })) ","pinned":false,"hidden":false,"mode":"ojs"},{"id":680,"value":"Plot.plot({\n style: {\"fontSize\" : 14, \"fontSize-swatches\" : 14, \"fontFamily\" : \"Helvetica\"},\n title : md`ω₃ **magnitude** by selection detection type`,\n round: true,\n color: {legend: true, \"domain\" : [\"BUSTED\",\"BUSTED-E\"], \"range\" : [\"darkgrey\",\"#ee0000\"]},\n height: 400,\n width : 600,\n marginTop: 40,\n marginBottom: 40,\n marginLeft: 50,\n paddingBottom : 10,\n facet: {\n data : plot_data,\n x: 'BUSTED',\n y: 'BUSTED-E',\n },\n x: {type: \"linear\", ticks : 5, 'label' : 'log10 (ω3)'},\n y: {type: \"sqrt\", ticks : 5, label : 'Proportion'},\n //x: {},\n marks: [\n Plot.rectY(plot_data, \n Plot.binX({y2: \"proportion-facet\"}, {x: \"omega\", fill: \"model\", inset: 0, thresholds: 100, mixBlendMode: \"multiply\"})),\n //Plot.density(plot_data, {x: \"omega\", fill: \"model\", opacity: 0.5, strokeWidth: 0.25, bandwidth: 10,mixBlendMode: \"multiply\"}),\n Plot.ruleY([0])\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1065,"value":"$32","pinned":false,"hidden":false,"mode":"ojs"},{"id":1230,"value":"viewof plot_model = Inputs.select ([\"BUSTED\", \"BUSTED-E\"], {\"label\" : \"Model to plot\"})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1237,"value":"viewof colors = Inputs.select ([\"turbo\", \"brbg\", \"prgn\", \"piyg\", \"puor\", \"rdbu\", \"rdgy\", \"rdylgn\", 'cividis', 'bugn', 'bupu', 'gnbu', 'orrd', 'pubu', 'purd', 'rdpu', 'ylgn', 'ylgnbu', 'ylorbr', 'ylorrd', 'magma'])","pinned":false,"hidden":false,"mode":"ojs"},{"id":1227,"value":"$33","pinned":false,"hidden":false,"mode":"ojs"},{"id":780,"value":"Plot.plot({\n style: {\"fontSize\" : 14, \"fontSize-swatches\" : 14, \"fontFamily\" : \"Helvetica\"},\n title : md`ω₃ **magnitude** by selection detection type`,\n round: true,\n //color: {legend: true, \"domain\" : [\"BUSTED\",\"BUSTED-E\"], \"range\" : [\"darkgrey\",\"#ee0000\"]},\n height: 400,\n width : 500,\n marginTop: 40,\n marginBottom: 100,\n marginLeft: 50,\n paddingBottom : 10,\n marginRight: 50,\n\n x: {'label' : '', tickRotate: 45},\n y: {type: \"sqrt\", ticks : 5, label : 'log10 (ω3)'},\n fx: {'label' : \"Detection class\"},\n //x: {},\n marks: [\n Plot.boxY(plot_data, {x: \"model\", y: \"omega\",\"fx\" : \"class\", \"fill\" : \"#DDD\", \"stroke\" : \"#444\"}), \n //Plot.density (plot_data, {x: \"model\", y: \"omega\",\"fx\" : \"class\"})\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1056,"value":"Plot.plot({\n style: {\"fontSize\" : 14, \"fontSize-swatches\" : 14, \"fontFamily\" : \"Helvetica\"},\n title : md`ω₃ **weight** by selection detection type`,\n round: true,\n //color: {legend: true, \"domain\" : [\"BUSTED\",\"BUSTED-E\"], \"range\" : [\"darkgrey\",\"#ee0000\"]},\n height: 400,\n width : 500,\n marginTop: 40,\n marginBottom: 100,\n marginLeft: 50,\n paddingBottom : 10,\n marginRight: 50,\n\n x: {'label' : '', tickRotate: 45},\n y: {type: \"sqrt\", ticks : 5, label : 'log10 (ω3)'},\n fx: {'label' : \"Detection class\"},\n //x: {},\n marks: [\n Plot.boxY(plot_data_w, {x: \"model\", y: \"omega\",\"fx\" : \"class\", \"fill\" : \"#DDD\", \"stroke\" : \"#444\"}), \n //Plot.density (plot_data, {x: \"model\", y: \"omega\",\"fx\" : \"class\"})\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":268,"value":"Plot.plot({\n style: {\"fontSize\" : 14, \"fontSize-swatches\" : 14, \"fontFamily\" : \"Helvetica\"},\n title : md`ω₃ **weight** by selection detection type`,\n round: true,\n color: {legend: true, \"domain\" : [\"BUSTED\",\"BUSTED-E\"], \"range\" : [\"darkgrey\",\"#ee0000\"]},\n height: 400,\n width : 600,\n marginTop: 40,\n marginBottom: 40,\n marginLeft: 50,\n paddingBottom : 10,\n facet: {\n data : plot_data_w,\n x: 'BUSTED',\n y: 'BUSTED-E',\n },\n x: {type: \"linear\", ticks : 5, 'label' : 'log10 P (ω3)'},\n y: {type: \"sqrt\", ticks : 5, label : 'Proportion'},\n //x: {},\n marks: [\n Plot.rectY(plot_data_w, \n Plot.binX({y2: \"proportion-facet\"}, {x: \"omega\", fill: \"model\", inset: 0, thresholds: \"freedman-diaconis\", mixBlendMode: \"multiply\"})),\n Plot.ruleY([0])\n ]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":694,"value":"annotator = (d, i, filtered) => {\n let ann = urls[dataset].annotate;\n if (ann) {\n let xt = urls[dataset][\"result-extension\"] || \"\";\n if (_.isString (ann)) {\n if (filtered) {\n return (html`${d}`);\n } else {\n return (html`${d}`);\n }\n } else {\n if (filtered) {\n return (html`${d}`);\n } else {\n return (html`${d}`);\n }\n }\n }\n return d;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":230,"value":"summary_stats = {\n const fmt = d3.format (\".2\");\n let get_stats = (vec)=> {\n return fmt (d3.mean (vec)) + \"/\" + fmt (d3.median (vec)) + \" [\" + fmt (d3.quantile (vec, .25)) + \",\" + fmt (d3.quantile (vec, .75)) + \"]\";}\n let res = {'N' : data.length};\n const sigl = 0.05;\n res[\"BUSTED+\"] = _.filter (data, (d)=>d.p <= sigl).length / res.N;\n res[\"BUSTED-E+\"] = _.filter (data, (d)=>d.p_e <= sigl).length / res.N;\n res[\"MA\"] = _.filter (data, (d)=>d.p_MA <= sigl).length / res.N;\n res[\"BUSTED-E preferred\"] = _.filter (data, (d)=>chi2cdf.default (2*(d.logL_e-d.logL), 2) >= 1- sigl).length / res.N;\n res['error'] = get_stats (_.map (data, (d)=>d.prop_e*100.));\n res['omega_BUSTED'] = get_stats (_.map (data, (d)=>d.omega3_BUSTED));\n res['omega_BUSTED-E'] = get_stats (_.map (data, (d)=>d[\"omega3_BUSTED-E\"]));\n res['p'] = get_stats (_.map (data, (d)=>100*d.p3_BUSTED));\n res['p-E'] = get_stats (_.map (data, (d)=>100*d[\"p3_BUSTED-E\"]));\n return res;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":236,"value":"null","pinned":false,"hidden":false,"mode":"ojs"},{"id":25,"value":"urls = ({\n//'PMID 30620335' : {'address' : \"https://data.hyphy.org/web/selection-data/30620335/30620335.csv\", 'annotate' : \"https://data.hyphy.org/web/selection-data/30620335/results\",'genes' : \"https://data.hyphy.org/web/selection-data/annotations/shultz-sackton.tsv.gz\"},\n//'zoonomia-32' : {'address' : \"https://data.hyphy.org/web/selection-data/zoonomia/32/32.csv\", 'annotate' : \"https://data.hyphy.org/web/selection-data/zoonomia/32/results\", 'result-extension' : '', 'genes' : (d)=>d.File.split ('.')[1]},\n'Endothermy' : {'address' : \"https://data.hyphy.org/web/busted-e/endothermy/busted-e_summary.csv\", 'annotate' : null, 'result-extension' : '', 'genes' : (d)=>d.File}})","pinned":false,"hidden":false,"mode":"ojs"},{"id":5,"value":"$34","pinned":false,"hidden":false,"mode":"ojs"},{"id":43,"value":"get2x2 = (data, pv, key1, key2)=> {\n let counts = [[0,0],[0,0]];\n let p = pv || 0.05;\n key1 = key1 || \"p\";\n key2 = key2 || \"p_e\";\n _.each (data, (d)=> {\n counts [d[key1] <= p ? 1 : 0][d[key2] <= p ? 1 : 0] += 1;\n });\n \n return [\n {'BUSTED' : 'No', 'BUSTED-E' : 'No', 'Count' : counts[0][0], 'p' : counts[0][0] / data.length * 100.},\n {'BUSTED' : 'Yes', 'BUSTED-E' : 'No', 'Count' : counts[1][0], 'p' : counts[1][0] / data.length * 100.},\n {'BUSTED' : 'No', 'BUSTED-E' : 'Yes', 'Count' : counts[0][1], 'p' : counts[0][1] / data.length * 100.},\n {'BUSTED' : 'Yes', 'BUSTED-E' : 'Yes', 'Count' : counts[1][1], 'p' : counts[1][1] / data.length * 100.},\n {'BUSTED' : 'Yes', 'BUSTED-E' : '*', 'Count' : counts[1][1] + counts[1][0], 'p' : (counts[1][1] + counts[1][0]) / data.length * 100.},\n {'BUSTED' : '*', 'BUSTED-E' : 'Yes', 'Count' : counts[0][1] + counts[1][1], 'p' : (counts[1][1] + counts[0][1]) / data.length * 100.},\n ];\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":60,"value":"import {pairplot} from '@chrispahm/plotting-pairwise-relationships-in-a-dataset'","pinned":false,"hidden":false,"mode":"ojs"},{"id":217,"value":"chi2cdf = import('https://cdn.skypack.dev/@stdlib/stats-base-dists-chisquare-cdf@0.1.0?min')","pinned":false,"hidden":false,"mode":"ojs"},{"id":718,"value":"import {obtainSourceData} from \"@spond/exploring-sars-cov-2-intrahost-variation-descriptive\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":270,"value":"vegalite = require (\"vega-embed\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":966,"value":"import {compute_corrections} from \"@spond/evaluate-different-forms-of-multiple-testing-correction\"","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W35","stdlib":"1","titling":"auto"}