30:T183a,genome_view = vegalite({ //"$schema": "https://vega.github.io/schema/vega-lite/v5.json", "config": { "view": { "stroke": null }, }, "vconcat": [ { "layer": [ { "height": 125, "width": 750, "data" : {"values" : _.filter (fragments, (d)=>!(d[1]<=view_this_range[0] && d[2] <= view_this_range[1] || d[2]>=view_this_range[0] && d[2] >= view_this_range[1]))}, "mark" : {"type" : "bar"}, "transform": [ { "calculate": "(datum['7'] % 2 == 0) ? -30 : -40", "as": "y" }, { "calculate": "(datum['7'] % 2 == 0) ? -40 : -50", "as": "y2" } ], "encoding" : { "tooltip": [{"field": "0", "type": "ordinal", "title" : "Fragment"}, {"field": "3", "type": "quantitative", title : "BUSTED p-value"}, {"field": "1", "type": "quantitative", title : "Start coordinate"}, {"field": "2", "type": "quantitative", title : "End coordinate"}, ], "x" : {"field" : "1", "type": "quantitative","axis": null, "scale" : {domain : view_this_range}}, "x2" : {"field" : "2", "type": "quantitative","axis": null}, "y" : {"field" : "y", "type": "quantitative","axis": null}, "y2" : {"field" : "y2", "type": "quantitative","axis": null}, "stroke" : { "condition": {"test": "datum['3'] <= 0.05", "value": "black"}, "value" : null }, "fill" : {"field" : "6", //"condition": {"test": "datum['3'] > 0.05", "value": "white"}, "legend" : {"title" : "Mean dN/dS internal branches", "orient" : "bottom"}, "type" : "quantitative", "scale": {"scheme": "lighttealblue", "type" : "sqrt", "domain" : [0,1], "clamp" : true}}, } }, //================================ { "data" : {"values" : _.filter (conserved_sites, (d)=>d[0] >= view_this_range[0] && d[0] <= view_this_range[1])}, "transform": [{ "calculate": "3", "as": "y" }, { "calculate": "6", "as": "y2" }], "mark": { type: "rule", tooltip : true, }, "encoding": { "x": { "field": "0", "type": "quantitative", "axis": null, }, "y": { "field": "y", "type": "quantitative", "axis": null, }, "y2": { "field": "y2", "type": "quantitative", "axis": null, }, "size" : {"value" : 0.05}, "color" : {"value" : "black"} } }, { "data" : {"values" : _.filter (negative_sites, (d)=>d[0] >= view_this_range[0] && d[0] <= view_this_range[1])}, "transform": [{ "calculate": "23", "as": "y" }, { "calculate": "28", "as": "y2" }], "mark": { type: "rule", tooltip : true, }, "encoding": { "x": { "field": "0", "type": "quantitative", "axis": null, }, "y": { "field": "y", "type": "quantitative", "axis": null, }, "y2": { "field": "y2", "type": "quantitative", "axis": null, }, "size" : {"value" : 0.5}, "color" : {"value" : "grey"}, "opacity" : {"value" : 0.8} } }, { "data" : {"values" : _.filter (selected_sites, (d)=>d.coordinate >= view_this_range[0] && d.coordinate <= view_this_range[1])}, "transform": [{ "calculate": "38", "as": "y" }, { "calculate": "43", "as": "y2" }], "mark": { type: "rule", tooltip : true, }, "encoding": { "x": { "field": "coordinate", "type": "quantitative", "axis": null, }, "y": { "field": "y", "type": "quantitative", "axis": null, }, "y2": { "field": "y2", "type": "quantitative", "axis": null, }, "size" : {"value" : 2}, "color" : {"field" : "branches", "type": "quantitative", "legend" : null, "scale": {"scheme": "lightgreyred"}} } }, { "data": { //"values": gene_coordinates "values" : _.filter (gene_coordinates, (d)=>!(d[0]<=view_this_range[0] && d[2] <= view_this_range[1] || d[0]>=view_this_range[0] && d[2] >= view_this_range[1])) }, "transform": [{ "calculate": "(datum['3'] % 2) ? -1 : 1", "as": "coloring" }, { "calculate": "(datum['3'] % 2) ? -5 : 0", "as": "y" }, { "calculate": "(datum['3'] % 2) ? 10 : 15", "as": "y2" } ], "mark": { "type": "bar", "stroke": "#000", }, "encoding": { "tooltip": [{"field": "0", "type": "quantitative", "title" : "Start"}, {"field": "2", "type": "quantitative", title : "End"}, {"field": "1", "type": "nominal", title : "Gene"}, ], "x": { "field": "0", "type": "quantitative", "axis": null }, "x2": { "field": "2" }, "y": { "field": "y", "type": "quantitative", "axis": null }, "y2": { "field": "y2" }, "fill": { "value" : null }, "opacity": { "value": "0.5" } } }, ] } ], //"resolve": {"scale": {"y": "shared"}}, })31:T593,table (busted_distribution_summary, { limit: 5000, enableCSVDownload: true, sort: true, enableFilter: false, columns: [ { key: 'file', name: 'Segment' }, { key: 'omega1', name: '⍵1', render : (val)=>d3.format (".3g") (val) }, { key: 'proportion1', name: 'p1', render : (val)=>d3.format (".3p") (val) }, { key: 'omega2', name: '⍵2', render : (val)=>d3.format (".3g") (val) }, { key: 'proportion2', name: 'p2', render : (val)=>d3.format (".3p") (val) }, { key: 'omega3', name: '⍵3', render : (val)=>d3.format (".3g") (val) }, { key: 'proportion3', name: 'p3', render : (val)=>d3.format (".3p") (val) }, { key : 'p', name : 'p-value', render : (val)=>val <= 0.05 ? '' + d3.format (".3g") (val) + '' : d3.format (".3g") (val) }, { key : 'q', name : 'q-value', render : (val)=>val <= 0.2 ? '' + d3.format (".3g") (val) + '' : d3.format (".3g") (val) }, ] });32:T10f9,{ let tree = tree_object[0]; var t = tree.render("#tree_container_0", { //height:1024, width:800, 'show-scale' : true, 'left-right-spacing': 'fit-to-size', //'top-bottom-spacing': 'fit-to-size', //'internal-names' : true, 'is-radial' : false, 'node_circle_size' : (n)=>0 }); const bs = (n) => { if (tree.is_leafnode (n)) { if (n.data.name == sc2_seq) { return 6; } return 3; } return 1; }; let siteID = summary_json[site2segment [site2view]? site2segment [site2view] : segment_id].subs[site2view]; let nodeStates = null; if (siteID) { nodeStates = {}; let node_labels = siteID; function recurse_node (node) { if (node.data.name in node_labels) { nodeStates[node.data.name] = node_labels[node.data.name]; } else { nodeStates[node.data.name] = nodeStates[node.parent.data.name]; } _.each (node.children, (c)=>recurse_node(c)); } recurse_node(tree.nodes); if (doAA) { for (let node in nodeStates) { nodeStates[node] = nodeStates[node] in translation_table ? translation_table[nodeStates[node]] : '?'; } } } let c10 = d3.scaleOrdinal(d3.schemeDark2); //var c10 = d3.scaleOrdinal(d3.schemeDark2); //t.node_label = n=>n.data.name.substr (0,3); t.node_bubble_size = bs; t.phylotree.node_bubble_size = bs; t.options["draw-size-bubbles"] = nodeStates ? true : false; t.font_size = 12; t.spacing_x (12); t.phylotree.branch_length_accessor = n=> { //console.log (summary_json[segment_id]); let bl = tree_object[1][n.data.name][2]; if (bl > 0.25) { n.data.censored_bl = true; return 0.25; } else { n.data.censored_bl = false; } return bl; } t.node_label ((n)=>{ if (nodeStates) { return nodeStates[n.data.name];// + " (" + n.data.name + ")"; } return n.data.name; //return ""; }); function sort_nodes (asc) { tree.traverse_and_compute (function (n) { var d = 1; if (n.children && n.children.length) { d += d3.max (n.children, function (d) { return d["count_depth"];}); } n["count_depth"] = d; }); tree.resort_children (function (a,b) { return (a["count_depth"] - b["count_depth"]) * (asc ? 1 : -1); }); } sort_nodes (true); tree.traverse_and_compute ((n) => { //if (!tree.is_leafnode (n) && !n.tagged_view && n.parent && n.parent.tagged_view) { //n.notshown = true; //} }); const is_focus = (n,d)=> { if (n == sc2_seq) { return 'Reference'; } return "Isolate"; } t.style_nodes ((e,n) => { if (n.data.name in tree_object[1]) { let binfo = tree_object[1][n.data.name]; if (nodeStates) { e.style ('fill', c10 (nodeStates[n.data.name])); if (binfo[0].length) { e.selectAll ("circle").style ('fill', species_colors (is_focus(n.data.name,binfo))); } } else { if (binfo[0].length) { e.style ('fill', species_colors (is_focus(n.data.name,binfo))); } } } //} }); t.style_edges ((e,d) => { //console.log (site); if (nodeStates) { if (nodeStates[d.target.data.name] != nodeStates[d.source.data.name]) { e.style ('stroke-width',5); } } else { if (d.target.data.name in tree_object[1]) { let binfo = tree_object[1][d.target.data.name]; if (binfo[0].length) { e.style ('stroke-color', species_colors (species_colors (is_focus(d.target.data.name,binfo)))); } if (tree_object[2][d.source.data.name]) { e.style ('stroke-width',5); } } //e.selectAll("title").text ("p = " + p); } if (d.target.data.censored_bl) { e.style ('stroke-dasharray',10); } }); t.placenodes(); t.update(); return t; }33:T5a6,fragments = { let result = _.sortBy(_.flatten (_.map (summary_json, (d,k)=> { let ranges = _.filter (d.map, (i)=>i>=0); let my_range = []; ranges.push (1e100); let last_check = ranges[0]; _.reduce (ranges, (accumulator, value)=> { if ((value - accumulator) % 3 == 0) { return value; } else { let tag = k.split ("_"); if (tag.length > 1) { tag = tag[0]+"/"+tag[1]; } else { tag = tag[0].split (".")[0]; } my_range.push ([k, last_check, accumulator + 2, d["busted"]["p"], 0, tag, d["rates"]["mean-omega"]["test"] ]); last_check = value; return value; } }, ranges[0]-3); return my_range; } ), true), (d)=>d[1]); let index = {}; let count = 0; _.each (result, (d) => { d.push (count++); }); return result; }34:T576,plotPie = (data)=>vegalite ({ "data": { "values": data }, "encoding": { "theta": {"field": "proportion", "type": "quantitative", "stack": true}, "color": {"field": "omega", "type": "quantitative", //"sort": "descending", "title" : "dN/dS", "scale": { "domain" : [0,4], "type" : "sqrt", "clamp" : true, //"domainMid": 1, "scheme": "lightmulti" } }}, "layer" : [ { "mark": {"type": "arc", "innerRadius": 40, "outerRadius" : 60, tooltip : true, padAngle : 0.0, strokeWidth : 1, strokeOpacity: 1}, "encoding" : { "stroke" : {value : "black"}, } }, { "mark": {"type": "text", "outerRadius" : 70, "fontSize" : 18}, "encoding": { "color" : {"value" : "black"}, "text": { "field": "omega", "type": "nominal", "format" : ".2f", "condition": {"test": "datum['omega'] <= 1.", "value": ""}, } } }, { "mark": {"type": "text", "outerRadius" : 20, "fontSize" : 18}, "encoding": { "color" : {"value" : "black"}, "text": { "field": "proportion", "type": "nominal", "format" : ".2p", "condition": {"test": "datum['omega'] <= 1.", "value": ""}, } } } ], "view": {"stroke": null} })35:[] 2e:{"title":"Visualizing selection analysis results for evolution of WNV","cells":[{"id":0,"value":"md`# Visualizing selection analysis results for evolution of WNV`","pinned":false,"hidden":false,"mode":"ojs"},{"id":2156,"value":"md`**Overview of selection regimes on the WNV isolates** \n\n- Sites found to be under positive selection by MEME on internal branches (colored by the number of branches showing evidence of diversifying selection).\n- Sites found to be under negative selection by FEL on internal branches. \n- Ticks in th middle of the bars show sites that are **conserved** on the amino-acid level among **all** sequences. \n- Non-recombinant fragments inferred by GARD; colors: mean dN/dS on the fragment; black outline -- segment is under positive selection.`","pinned":false,"hidden":false,"mode":"ojs"},{"id":2584,"value":"viewof coord_range = text ({\n title: \"View this coordinate range\",\n value: \"1-15000\",\n submit: true\n});","pinned":false,"hidden":false,"mode":"ojs"},{"id":627,"value":"$30","pinned":false,"hidden":false,"mode":"ojs"},{"id":2167,"value":"md`Mean dN/dS and total branch length for each of internal (test) and terminal (background) branches`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1807,"value":"segment_profile = vegalite({\n\t\"$schema\": \"https://vega.github.io/schema/vega-lite/v5.json\",\n\t\"config\": {\n\t\t\"view\": {\n\t\t\t\"stroke\": null\n\t\t},\n\t},\n \"data\" : {\"values\" : selection_profile},\n \"mark\" : {\"type\" : \"line\", \"clip\": true, \"tooltip\": {\"content\": \"data\"}},\n \"width\" : 300,\n \"height\" : 50,\n \"encoding\": {\n \"facet\": {\n \"field\": \"s\",\n \"type\": \"ordinal\",\n \"columns\": 2,\n \"title\" : \"Branch set\",\n }, \n \"x\" : {\"field\" : \"x\", \"type\": \"quantitative\", \"axis\": {\"title\" : \"Coordinate\"}}, \n \"y\" : {\"field\" : \"omega\", \"type\": \"quantitative\", \"scale\": {\"clip\": true}, \"axis\": {\"title\" : \"Mean dN/dS\"}}, \n //\"color\" : {\"field\" : \"s\", \"legend\" : {\"title\" : \"Species\"}}\n }\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1826,"value":"segment_profile_T = vegalite({\n\t\"$schema\": \"https://vega.github.io/schema/vega-lite/v5.json\",\n\t\"config\": {\n\t\t\"view\": {\n\t\t\t\"stroke\": null\n\t\t},\n\t},\n \"data\" : {\"values\" : selection_profile},\n \"mark\" : {\"type\" : \"line\", \"clip\": true},\n \"width\" : 300,\n \"height\" : 50,\n \"encoding\": {\n \"facet\": {\n \"field\": \"s\",\n \"type\": \"ordinal\",\n \"columns\": 2,\n \"title\" : \"Branch set\"\n }, \n \"x\" : {\"field\" : \"x\", \"type\": \"quantitative\", \"axis\": {\"title\" : \"Coordinate\"}}, \n \"y\" : {\"field\" : \"T\", \"type\": \"quantitative\", \"scale\": { \"clip\": true}, \"axis\": {\"title\" : \"Tree length\"}}, \n //\"color\" : {\"field\" : \"s\", \"legend\" : {\"title\" : \"Species\"}}\n }\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2233,"value":"md`### Segment information`","pinned":false,"hidden":false,"mode":"ojs"},{"id":966,"value":"viewof segment_id = {\n const dd3 = select({\n title: \"View this segment\",\n options: _.map (_.uniqBy (fragments,\"0\"), (g)=>g[0])\n \n });\n dd3.input.style.fontSize = \"15px\";\n dd3.input.style.marginTop = \"8px\";\n return dd3;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2137,"value":"md`ω distribution inferred for internal branches on this segment by BUSTED`","pinned":false,"hidden":false,"mode":"ojs"},{"id":2134,"value":"plotPie (_.filter (busted_distributions, (d)=>d.file == segment_id))","pinned":false,"hidden":false,"mode":"ojs"},{"id":2236,"value":"md`### Site-level information`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1972,"value":"\nviewof site2viewExpanded = autoSelect({\n options: _.flatten ([\"Overall for \" + segment_id,_.map (site2segment, (s,k)=>k + \" \" + site2gene[k])]),\n title: `View this site (in NCBI reference coordinates)`,\n placeholder: \"Select a site (or overall)\"\n})\n","pinned":false,"hidden":false,"mode":"ojs"},{"id":2067,"value":"viewof doAA = checkbox({\n options: [{ value: \"toggle\", label: \"Amino-acids\" }]\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2301,"value":"site2view && site2view != \"Overall\" ? md`Site **${site2viewExpanded}** is in segment *${site2segment[site2view]}* and is codon ${site2alignment[site2view]} in the overall alignment. \n
FEL p-value ${float_format (annotation_json[site2view].bFEL.p)}, ω = ${report_ratio (annotation_json[site2view].bFEL.a, annotation_json[site2view].bFEL.b)}\n
MEME p-value ${float_format (annotation_json[site2view].bMEME.p)}, ω + = ${report_ratio (annotation_json[site2view].bMEME[\"a\"], annotation_json[site2view].bMEME[\"b+\"])}, p = ${float_format (annotation_json[site2view].bMEME[\"w+\"])}, **${annotation_json[site2view].bMEME.br}** selected branches\n
Evolutionary credibility report: ${evo_credibility_report(site2view)}\n
**Amino-acid composition** \n${freq_report (annotation_json [site2view].baa[\"background\"])} \n` : md`Thicker branches are those that have experienced episodic selection according to the aBSREL method`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1941,"value":"trendLegend = swatches({\n color: species_colors,\n title: \"Viral species\",\n //tickFormat: \"%\"\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":15,"value":"html`\n \n

\n`\n","pinned":false,"hidden":false,"mode":"ojs"},{"id":870,"value":"html`\n \n

`","pinned":false,"hidden":false,"mode":"ojs"},{"id":2239,"value":"md`### Summary tables`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1693,"value":"md`**Table 1** List of **${selected_sites.length}** sites found to be under selection by MEME (p ≤ ${pv_cut}) along internal branches, as well as biochemical properties that are _important_ at this site (via the PRIME method)`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1028,"value":"table (selected_sites,\n {\n limit: 5000,\n enableCSVDownload: true,\n sort: true,\n enableFilter: false,\n columns: [\n {\n key: 'coordinate',\n name: 'Coordinate (reference)'\n },\n {\n key: 'gene',\n name: 'Gene/ORF'\n \n },\n {\n key: 'site',\n name: 'Codon (in alignment)'\n },\n {\n key : 'meme_weight',\n name : '% of branches with ω>1',\n render : (val)=>d3.format (\".3p\") (val)\n },\n {\n key : 'p',\n name : 'p-value',\n render : (val)=>d3.format (\".3g\") (val)\n },\n {\n key : 'prime',\n name : 'Properties',\n render : (val)=>_.map (val, (d)=>{if (d[2]) return d[0] + (d[2]>0 ? \" (constrained) \" : \" (changing)\"); else return 'overall';}).join (\", \")\n }\n ]\n});","pinned":false,"hidden":false,"mode":"ojs"},{"id":2244,"value":"md`**Table 2** List of ${negative_sites.length} sites found to be under negative selection by FEL (p≤${pv_cut}) along internal branches clade`","pinned":false,"hidden":false,"mode":"ojs"},{"id":2249,"value":"table (negative_sites,\n {\n limit: 5000,\n enableCSVDownload: true,\n sort: true,\n enableFilter: false,\n columns: [\n {\n key: '0',\n name: 'Coordinate (reference)',\n render : (val)=>val + \" \" + site2gene [val]\n },\n {\n key: '3',\n name: 'Codon (in alignment)'\n },\n {\n key : '1',\n name : 'dN/dS',\n render : (val)=>d3.format (\".3g\") (val)\n },\n {\n key : '2',\n name : 'p',\n render : (val)=>d3.format (\".3g\") (val)\n }\n ]\n});","pinned":false,"hidden":false,"mode":"ojs"},{"id":2173,"value":"md`**Table 3** BUSTED selection results on internal branches across segments`","pinned":false,"hidden":false,"mode":"ojs"},{"id":1663,"value":"$31","pinned":false,"hidden":false,"mode":"ojs"},{"id":2563,"value":"md`**Table 4** aBSREL selection analysis summary on internal branches (p≤${pv_cut}, following multiple testing correction)`","pinned":false,"hidden":false,"mode":"ojs"},{"id":2633,"value":"table (absrel_summary,\n {\n limit: 5000,\n enableCSVDownload: true,\n sort: true,\n enableFilter: false,\n columns: [\n {\n key: 'segment',\n name: 'Segment'\n \n },\n {\n key: 'selected',\n name: '# of selected branches'\n },\n {\n key: 'mean',\n name: 'Mean fraction of sites under selection',\n render : (val)=>val ? d3.format (\".3p\") (val) : \"N/A\"\n },\n ]\n});","pinned":false,"hidden":false,"mode":"ojs"},{"id":2863,"value":"md`**Table 5** FADE selection analysis summary (directional selection) on internal branches (BF≥${BF_cut})`","pinned":false,"hidden":false,"mode":"ojs"},{"id":2871,"value":"table (fade_sites,\n {\n limit: 5000,\n enableCSVDownload: true,\n sort: true,\n enableFilter: false,\n columns: [\n {\n key: 'coord',\n name: 'Coordinate (reference)'\n },\n {\n key: 'site',\n name: 'Gene/ORF'\n \n },\n {\n key : 'Res',\n name : 'Target'\n },\n {\n key : 'BF',\n name : 'Bayes Factor',\n render : (val)=>d3.format (\".3g\") (val)\n },\n {\n key : 'Bias',\n name : 'Rate bias',\n render : (val)=>d3.format (\".3g\") (val)\n }\n ]\n});","pinned":false,"hidden":false,"mode":"ojs"},{"id":5,"value":"md`---`","pinned":false,"hidden":false,"mode":"ojs"},{"id":2809,"value":"md`### Data and libraries`","pinned":false,"hidden":false,"mode":"ojs"},{"id":2588,"value":"view_this_range = {\n let range = coord_range ? _.map (coord_range.split (\"-\"), (d)=>+d) : [];\n if (range.length == 2 && range[0] < range[1] && range[0] > 0 && range[1] < 29903) {\n return range;\n }\n return [1,29903];\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":9,"value":"d3 = require (\"d3\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":11,"value":"phylotree=require (\"https://dl.dropboxusercontent.com/s/i528j2lqnyzknit/phylotree.js\")\n//require (\"phylotree/build/phylotree.js\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":13,"value":"_=require(\"lodash\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":78,"value":"//resolve(\"phylotree@0.1/phylotree.css\")\n//require (\"https://dl.dropbox.com/s/mjhii51tjy1o5gt/phylotree.js\")\nhtml`\n\n`","pinned":false,"hidden":false,"mode":"ojs"},{"id":2524,"value":"site2view = site2viewExpanded.split (\" \")[0]","pinned":false,"hidden":false,"mode":"ojs"},{"id":1905,"value":"tree_object = [new phylotree.phylotree(summary_json[site2segment [site2view]? site2segment [site2view] : segment_id][\"tree\"]), summary_json[site2segment [site2view]? site2segment [site2view] : segment_id][\"tree_tags\"], summary_json[segment_id][\"aBSREL\"] || {}]","pinned":false,"hidden":false,"mode":"ojs"},{"id":127,"value":"$32","pinned":false,"hidden":false,"mode":"ojs"},{"id":874,"value":"import { table } from \"@gampleman/table\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":174,"value":"import {select,checkbox, autoSelect, text} from \"@jashkenas/inputs\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":488,"value":"vegalite = require(\"vega-embed\")","pinned":false,"hidden":false,"mode":"ojs"},{"id":1793,"value":"selection_profile = {\n let summary = [];\n _.each (fragments, (d,i)=> {\n _.each (species, (s)=> {\n let T = summary_json[d[0]].rates[\"T\"][s];\n if (T > 0) {\n summary.push ({'segment' : i+1, 'x': d[1], 'omega' : summary_json[d[0]].rates[\"mean-omega\"][s], 'T' : T, 's': s});\n summary.push ({'segment' : i+1, 'x': d[2], 'omega' : summary_json[d[0]].rates[\"mean-omega\"][s], 'T' : T, 's': s});\n\n } else {\n summary.push ({'segment' : i+1, 'x': d[1], 'T' : T, 's': s});\n summary.push ({'segment' : i+1, 'x': d[2], 'T' : T, 's': s});\n\n }\n })});\n return summary;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":684,"value":"$33","pinned":false,"hidden":false,"mode":"ojs"},{"id":2288,"value":"float_format = d3.format (\".3g\");","pinned":false,"hidden":false,"mode":"ojs"},{"id":2217,"value":"$34","pinned":false,"hidden":false,"mode":"ojs"},{"id":1876,"value":"prime_annotation = _.map (summary_json[_.first (_.keys(summary_json))][\"prime-properties\"], (d)=>d.split (\" \")[2])","pinned":false,"hidden":false,"mode":"ojs"},{"id":2321,"value":"report_ratio = (a,b)=> {\n if (a == 0 && b == 0) return 'N.D.';\n if (a == 0) return 'Infinite';\n return float_format (b/a);\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":1929,"value":"species = _.sortBy (_.filter (_.uniq (_.map (summary_json[_.first (_.keys(summary_json))][\"tree_tags\"], (d)=>d[0])), (d)=>d.length))","pinned":false,"hidden":false,"mode":"ojs"},{"id":1937,"value":"species_colors = d3.scaleOrdinal(d3.schemeCategory10).domain ([\"Isolate\",\"Reference\"]);","pinned":false,"hidden":false,"mode":"ojs"},{"id":1152,"value":"busted_distributions = _.flatten (_.map (summary_json, (d,k)=>{\n return _.map (d[\"rates\"][\"busted\"][\"Test\"], (dd)=>{\n return {\n file: k,\n omega: dd.omega,\n proportion: dd.proportion,\n p: d[\"busted\"][\"p\"]\n };\n });\n})); ","pinned":false,"hidden":false,"mode":"ojs"},{"id":2230,"value":"relax_distributions = _.flatten (_.map (summary_json, (d,k)=>{\n return _.flatten (_.map (d[\"rates\"][\"relax\"], (dd2,kk)=>{\n \n return _.map (dd2, (dd)=>{return {\n file: k,\n species: kk,\n omega: dd.omega,\n proportion: dd.proportion,\n p: d[\"relax\"][\"p\"]\n }\n });\n }))})); ","pinned":false,"hidden":false,"mode":"ojs"},{"id":1666,"value":"busted_distribution_summary = {\n let bds = _.sortBy (_.map (summary_json, (d,k)=>{\n let dd = d[\"rates\"][\"busted\"][\"Test\"];\n return {\n file: k,\n omega1: dd[0].omega,\n proportion1: dd[0].proportion,\n omega2: dd[1].omega,\n proportion2: dd[1].proportion,\n omega3: dd[2].omega,\n proportion3: dd[2].proportion,\n p: d[\"busted\"][\"p\"],\n };\n }), (d)=>d.p); \n return _.map (bds, (record, i)=> {\n record['q'] = record['p'] * bds.length / (i+1); \n return record;\n });\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":764,"value":"conserved_sites = _.filter (_.map (annotation_json, (v,k)=>[k,v]), (s)=>_.uniq (_.flatten (_.map (species, (sp)=>_.uniq (_.filter (_.keys (s[1].baa[sp]), (k)=>k!='-'))))).length == 1, (d)=>[+d[0]])","pinned":false,"hidden":false,"mode":"ojs"},{"id":2195,"value":"site2segment = {\n let r = {};\n _.each (summary_json, (d,s)=> {\n _.each (d.map, (c)=> {\n if (c >= 0) {\n r[c] = s; \n }\n });\n });\n \n return r;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2407,"value":"site2alignment = {\n let r = {};\n let offset = 1;\n let last_segment = \"\";\n let only_sites = _.sortBy (_.map (annotation_json, (d,k)=>[k,d.index]), (s)=>+s.k);\n _.each (only_sites, (d)=> {\n if (site2segment[d[0]] != last_segment) {\n if (last_segment && last_segment.length) {\n offset += summary_json [last_segment].map.length;\n }\n last_segment = site2segment[d[0]];\n }\n r [d[0]] = offset + d[1];\n });\n \n return r;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":780,"value":"selected_sites = _.map (_.filter (_.map (annotation_json, (v,k)=>[k,v]), (s)=>s[1].bMEME.p<=pv_cut), (d)=>{return {'coordinate' : +d[0],'p' : d[1].bMEME.p,'meme_weight' : d[1].bMEME[\"w+\"],'gene': d[1].G + \"/\" + d[1].S, 'site' : site2alignment[d[0]], 'fel_p' : d[1].bFEL.p, 'fel_n' : d[1].bFEL.b < d[1].bFEL.a, 'branches' : d[1].bMEME.br, 'prime' : _.map(_.filter (_.map (d[1].bPRIME ? d[1].bPRIME.p : [], (p,i)=>[i,p,i>0 ? d[1].bPRIME['lambda'][i-1] : null]), (pr)=>pr[1]{if (d[0]>0) d[0] = prime_annotation [d[0]-1]; else d [0] = \"Overall\"; return d;})};})","pinned":false,"hidden":false,"mode":"ojs"},{"id":2816,"value":"fade_sites = _.chain (annotation_json).toPairs().filter (d=>'bFADE' in d[1]).map (d=>[d[0], d[1].G+\"/\"+d[1].S,_.chain(d[1].bFADE).toPairs().map (d=> ({'Res' : d[0], 'Bias' : d[1].r, 'BF' : d[1].bf})).value()]).map (d=>_.map (d[2], r=>_.extend ({'coord' : d[0], 'site' : d[1],}, r))).flatten().value()","pinned":false,"hidden":false,"mode":"ojs"},{"id":2776,"value":"pv_cut = 0.05","pinned":false,"hidden":false,"mode":"ojs"},{"id":2817,"value":"BF_cut = 100","pinned":false,"hidden":false,"mode":"ojs"},{"id":798,"value":"negative_sites = _.map (_.filter (_.map (annotation_json, (v,k)=>[k,v]), (s)=>_.size (s[1].bcdn[species[0]]) > 1 && s[1].bFEL.p<=pv_cut && s[1].bFEL.a > s[1].bFEL.b), (d)=>[+d[0], d[1].bFEL.b / d[1].bFEL.a, d[1].bFEL.p, site2alignment[d[0]]])","pinned":false,"hidden":false,"mode":"ojs"},{"id":2574,"value":"credibility_heatmap = {\n return _.map (_.filter (_.map (annotation_json, (d,s)=> {\n if (\"evo\" in d) {\n return [s, \n d3.sum (_.map (d[\"evo\"], (d)=>-d * Math.log2 (d)))\n \n ];\n }\n return [s, null];\n }), (d)=>d[1]), (d)=>{return {'coordinate' : d[0], 'entropy' : d[1]}});\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2663,"value":"majority_discordant = {\n return _.map (_.filter (_.map (annotation_json, (d,s)=> {\n if (\"evo\" in d) {\n let most_common = _.chain (d[\"evo\"]).toPairs().sortBy (d=>-d[1]).value()[0];\n if (most_common[0] != d[\"bSC2\"] && most_common[1] > 0.9) {\n return [s, \n most_common,\n d[\"bSC2\"]\n ];\n }\n }\n return [s, null];\n }), (d)=>d[1]), (d)=>{return {'coordinate' : d[0], 'discordant' : [d[1],d[2]]}});\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2078,"value":"site_report = annotation_json [site2view]\n ","pinned":false,"hidden":false,"mode":"ojs"},{"id":2339,"value":"freq_report = (freqs)=> {\n return _.map (_.sortBy (_.map (freqs, (f,a)=>[a,f]), (d)=>-d[1]), (d)=> d[0] + \"_{\" + d[1] + \"}\").join (\" \");\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2517,"value":"site2gene = {\n let r = {};\n _.each (annotation_json, (d,s)=> {\n r[s] = d.G + \"/\" + d.S;\n });\n \n return r;\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2354,"value":"cfel_site_report = (d)=> {\n const site_data = annotation_json[d][\"bCFEL\"];\n return \"\" + _.map (species, (sp)=>\"\"). join (\"\") + \"

Species	ω	Substitutions
\" + sp + \"	\" + report_ratio (site_data.a, site_data.b[sp]) + \"	\" + site_data.s[sp] +\"

\" + \n \"\" +\n _.map (site_data.p, (p,k)=> {return \"\"}).join (\"\")\n +\"

Pair	p-value
\" + k + \"	\" + float_format(p) + \"

\";\n}\n\n","pinned":false,"hidden":false,"mode":"ojs"},{"id":2534,"value":"evo_credibility_report = (d)=> {\n let site_data = annotation_json[d][\"evo\"];\n if (!site_data) return \"No data (invariable site)\";\n site_data = _.chain (site_data).toPairs().sortBy ((d)=>-d[1]).value();\n return \"\" + _.map (site_data, (d)=>\"\"). join (\"\") + \"

Codon	AA	Predicted probability in SARS-CoV-2
\" + d[0] + \"	\" + translation_table[d[0]] + \"	\" + float_format (d[1]) +\"

\";\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":2621,"value":"absrel_summary = _.chain (summary_json).toPairs ().map ((d)=>{return {'segment':d[0], 'selected' :_.size(d[1].aBSREL), 'mean' :d3.mean (_.map (d[1].aBSREL, (r)=>r.r[1][1]))}}).sortBy (d=>-d.selected).value()","pinned":false,"hidden":false,"mode":"ojs"},{"id":2497,"value":"summary_json = summary_json_2[1]","pinned":false,"hidden":false,"mode":"ojs"},{"id":2500,"value":"sc2_seq = summary_json_2[0]","pinned":false,"hidden":false,"mode":"ojs"},{"id":681,"value":"summary_json_2 = {\n let raw_res = await d3.json (\"https://dl.dropbox.com/s/pvwuytsyw8ne89o/summary.json\");\n let sc2_seq = raw_res['SC2'];\n delete raw_res['SC2'];\n let gene_coordinates = raw_res ['gene-coordinates'];\n delete raw_res['gene-coordinates'];\n raw_res = _.chain (raw_res).toPairs ().filter (d=>d[1].busted).fromPairs().value();\n return [sc2_seq,raw_res,gene_coordinates];\n}\n ","pinned":false,"hidden":false,"mode":"ojs"},{"id":754,"value":"annotation_json = {\n let r = await d3.json (\"https://dl.dropbox.com/s/w1ty3qc0zbv7cvm/annotation.json\");\n return _.chain (r).toPairs ().filter ((d)=>\"baa\" in d[1]).fromPairs().value();\n}","pinned":false,"hidden":false,"mode":"ojs"},{"id":629,"value":"gene_coordinates = _.map (summary_json_2[2], (d,i)=>{\n let dd = [d[0],d[2],d[1],i];\n return dd;\n})","pinned":false,"hidden":false,"mode":"ojs"},{"id":1598,"value":"name_fragment = (n)=>{ \n let pieces = n.split ('_');\n if (pieces.length > 3) {\n return pieces[1]+\"/\"+pieces[2];\n }\n return pieces[1];\n};","pinned":false,"hidden":false,"mode":"ojs"},{"id":1943,"value":"import { legend, swatches } from \"@d3/color-legend\"","pinned":false,"hidden":false,"mode":"ojs"},{"id":2060,"value":"import {translation_table} from \"@spond/usual-vs-unusual-mutation-analyses\"","pinned":false,"hidden":false,"mode":"ojs"}],"files":"$W35","stdlib":"1","titling":"auto"}