28:null 29:["$","title",null,{"children":"DRHIP result viewer | ACME Computational Molecular Evolution Group | Observable"}] 2e:Tdf0,vegalite({ "$schema": "https://vega.github.io/schema/vega-lite/v5.json", "data": { "values": combined_summary }, "width": 1000, // Change this value to increase/decrease facet width "height": 600, // Optional: Adjust height if needed "layer": [ { "mark": "bar", "transform": [ {"calculate": "log(datum.RELAX_K)", "as": "log_k"}, ], "encoding": { "x": { "field": "gene", "type": "nominal", "title": "Gene", "axis": {"labelAngle": 0, "labelFontSize": 14, "titleFontSize": 16, } }, "y": { "field": "log_k", "type": "quantitative", "title": "log(RELAX K)", "axis": { "labelFontSize": 14, "titleFontSize": 16, }, "scale": { "domain": [extent[0]-0.5, extent[1]+0.5], } }, "opacity": { "condition": {"test": "datum.RELAX_overall_pval < 0.01", "value": 0.95}, "value": 0.45 } } }, { "mark": { "type": "text", "dy": {"expr": "log(datum.RELAX_K) >= 0 ? -25 : 10"}, "fontSize": 14, "fontWeight": "bold" }, "transform": [{"calculate": "log(datum.RELAX_K)", "as": "log_k"}, {"calculate": "'K=' + format(datum.RELAX_K, '.3f')", "as": "formatted_Kval"}, ], "encoding": { "x": {"field": "gene", "type": "nominal"}, "y": {"field": "log_k", "type": "quantitative"}, "text": {"field": "formatted_Kval", "type": "nominal"}, // "color": { // "condition": {"test": "datum.RELAX_overall_pval < 0.01", "value": "darkorange"}, // "value": "black" // } } }, { "mark": { "type": "text", "dy": {"expr": "log(datum.RELAX_K) >= 0 ? -10 : 25"}, "fontSize": 14, "fontWeight": "bold" }, "transform": [{"calculate": "log(datum.RELAX_K)", "as": "log_k"}, {"calculate": "'p=' + (datum.RELAX_overall_pval >= 0.01 ? format(datum.RELAX_overall_pval, '.2f') : format(datum.RELAX_overall_pval, '.2e'))", "as": "formatted_pval"} ], "encoding": { "x": {"field": "gene", "type": "nominal"}, "y": {"field": "log_k", "type": "quantitative"}, "text": {"field": "formatted_pval", "type": "nominal"}, // "color": { // "condition": {"test": "datum.RELAX_overall_pval < 0.01", "value": "#b47846"}, // "value": "black" // } } }, { "mark": { "type": "rule", "color": "black", //"strokeDash": [6,4], "strokeWidth": 1.5 }, "encoding": { "y": { "datum": 0 } } } ], "config": { "font": "Helvetica", "axis": { "labelFontSize": 10, "titleFontSize": 11, "tickSize": 3, "domain": true, "grid": false }, "header": { "labelFontSize": 11, "titleFontSize": 12, "labelFontWeight": "bold" }, "view": { "stroke": null }, //"facet": { "spacing": 10 } } })2f:Ta7b,**FEL** infers nonsynoymous (dN) and synonymous (dS) substitution rates for each site in a gene alignment. FEL assumes that the selection pressure for each site is constant across the entire phylogeny; _i.e._, a site must have evidence of negative selection across _all branches_ to have dN/dS < 1. Here, we report 'positive' if dN/dS > 1 and p<0.05, 'negative' if dN/dS < 1 and p<0.05, and 'neutral' if p>=0.05. **MEME** indicates that a site has been subject to episodic positive selection at some point in its history. Unlike BUSTED, MEME is not a gene-wide test for positive selection and is not subject to multiple testing correction, so positive results may be false positives and should be treated with some caution if BUSTED didn't show a significant result for the gene in question. MEME is also a better choice than FEL when searching for evidence of site-specific positive selection, because MEME does not require that all branches have evidence of positive selection for a site to be significant. **Contrast-FEL** searches for evidence that a site may be associated with adaptation to a different environment. A significant result at a site means that dN/dS (β/α) is different between the two sets of branches, with either an increase or a decrease on the foreground branches relative to the reference branches. A significant finding does not make any claims about positive (dN/dS > 1) or negative (dN/dS < 1), just that dN/dS differ among sets of branches, _i.e._, a difference need not change the mode of selection. Non-significant results do not necessarily mean that there is no difference in selection pressure between foreground and reference branches; it is also possible that whatever difference there may be does not rise to the level of statistical significance. **PRIME** searches for evidence that a site is experiencing selective pressure to maintain a particular physicochemical property. We have chosen to focus on the set of five biochemical properties described by Atchley et. al. (2005) when assessing property conservation at each site. **Intensified positive selection:** A site with evidence of episodic positive selection (_i.e._, a significant MEME p-value) and a statistically significant increase in the nonsynonymous substitution rate in the foreground relative to the reference (Contrast-FEL significant p-value and Beta_foreground > Beta_reference) may be subject to intensified positive selection. Such sites may also have different majority amino acid residues in the foreground and reference branches. Filters with a * require that a foreground branch set was specified during the analysis, and that Contrast-FEL results are present.30:T476,vegalite({ $schema: 'https://vega.github.io/schema/vega-lite/v5.json', description: 'Treemap of Composition', title: "Site composition", width: 350, // Change this value to increase/decrease facet width height: 350, // Optional: Adjust height if needed data: { values: focal_site_composition }, mark: { type: 'rect'}, encoding: { x: { field: 'category', type: 'ordinal', title: "Amino acid", axis: {"labelAngle": 0, "labelFontSize": 15}, sort: {x: "y"} }, y: { field: 'size', type: 'quantitative', title: "# Branches" }, // size: { field: 'size', type: 'quantitative' }, text: { field: 'size', type: 'quantitative' } }, layer: [ { mark: 'rect', encoding: { color: { field: 'category', type: 'nominal', legend: false }, size: { field: 'size', type: 'quantitative' } } }, { mark: { type: 'text', align: 'center', baseline: 'bottom' }, encoding: { text: { field: 'size', type: 'quantitative' }, } }, ] })31:T49b,vegalite({ $schema: 'https://vega.github.io/schema/vega-lite/v5.json', description: 'most common substitutions at the current site', title: "Most Common Substitutions", width: 350, // Change this value to increase/decrease facet width height: 350, // Optional: Adjust height if needed data: { values: focal_site_substitutions }, mark: { type: 'rect'}, encoding: { x: { field: 'category', type: 'ordinal', title: "Amino acid", axis: {"labelAngle": 0, "labelFontSize": 15}, sort: {x: "y"} }, y: { field: 'size', type: 'quantitative', title: "Substitutions" }, // size: { field: 'size', type: 'quantitative' }, text: { field: 'size', type: 'quantitative' } }, layer: [ { mark: 'rect', encoding: { color: { field: 'category', type: 'nominal', legend: false }, size: { field: 'size', type: 'quantitative' } } }, { mark: { type: 'text', align: 'center', baseline: 'bottom' }, encoding: { text: { field: 'size', type: 'quantitative' }, } }, ] })32:[] 2c:{"id":"df8bfd20604373de","slug":"drhip-result-viewer","title":"DRHIP result viewer","thumbnail":"79a7a73680f7faa92a882ffc559ebd1ab80283e1d2539dd4114300a58552ce36","access_level":"public","create_time":"$D2025-09-18T20:55:40.290Z","update_time":"$D2026-07-06T20:10:33.223Z","view_time":null,"first_public_time":"$D2025-09-22T16:51:20.062Z","has_unpublished_edits":null,"fork":true,"stars":0,"starred":false,"owner":{"type":"team","id":"7ada8a89538f700a","login":"hyphy","name":"ACME Computational Molecular Evolution Group","avatar_url":"https://avatars-next.observableusercontent.com/c49d98d47337d40736a6fb1a598558be99c1e5fbeea48e7cac329854f99ccb98","tier":"education_2026"},"version":1016,"publish_version":null,"role":null,"authored":false,"authors":[{"type":"individual","id":"d7d30d996f330cf9","login":"hverdonk","name":"Hannah Verdonk","avatar_url":"https://avatars-next.observableusercontent.com/e3886da95f22b11c72bd727c03af4612d64ebac68ae9c2a7a87ae5122c0d27f3","tier":null}],"members":[{"type":"individual","id":"d7d30d996f330cf9","login":"hverdonk","name":"Hannah Verdonk","avatar_url":"https://avatars-next.observableusercontent.com/e3886da95f22b11c72bd727c03af4612d64ebac68ae9c2a7a87ae5122c0d27f3","tier":null,"role":"editor"},{"type":"team","id":"7ada8a89538f700a","login":"hyphy","name":"ACME Computational Molecular Evolution Group","avatar_url":"https://avatars-next.observableusercontent.com/c49d98d47337d40736a6fb1a598558be99c1e5fbeea48e7cac329854f99ccb98","tier":"education_2026","role":"editor"}],"comments":[{"id":"324cff49ca6f719f","reply_to_id":null,"content":"add majority residue per-branch-set","create_time":"$D2025-09-19T18:54:56.614Z","update_time":null,"resolved":true,"user":{"type":"individual","id":"d7d30d996f330cf9","login":"hverdonk","name":"Hannah Verdonk","avatar_url":"https://avatars-next.observableusercontent.com/e3886da95f22b11c72bd727c03af4612d64ebac68ae9c2a7a87ae5122c0d27f3","tier":null}}],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":2,"collection":null,"type":"notebook","body":{"title":"DRHIP result viewer","cells":[{"id":0,"value":"# DRHIP result viewer","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":541,"value":"// uncomment to view results reported in \"CAPHEINE or everything and the kitchen sink: a workflow for automated selection analyses using HyPhy\"\n//combined_summary = await d3.csv(\"https://raw.githubusercontent.com/hverdonk/CAPHEINE-H5N1/refs/heads/main/drhip/combined_summary.csv\", d3.autoType);\nviewof combined_summary_file = Inputs.file({label: \"Combined Summary\", accept: \".csv\", required: true})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":544,"value":"// combined_sites = await d3.csv(\"https://raw.githubusercontent.com/hverdonk/CAPHEINE-H5N1/refs/heads/main/drhip/combined_sites.csv\", d3.autoType);\nviewof combined_sites_file = Inputs.file({label: \"Combined Sites\", accept: \".csv\", required: true})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":548,"value":"//combined_comparison_summary = await d3.csv(\"https://raw.githubusercontent.com/hverdonk/CAPHEINE-H5N1/refs/heads/main/drhip/combined_comparison_summary.csv\", d3.autoType);\nviewof combined_comparison_summary_file = Inputs.file({label: \"Combined Comparison Summary (optional)\", accept: \".csv\", required: false})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":546,"value":"//combined_comparison_sites = await d3.csv(\"https://raw.githubusercontent.com/hverdonk/CAPHEINE-H5N1/refs/heads/main/drhip/combined_comparison_site.csv\", d3.autoType);\nviewof combined_comparison_sites_file = Inputs.file({label: \"Combined Comparison Sites (optional)\", accept: \".csv\", required: false})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":117,"value":"## Data summary","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":305,"value":"`positive_sites` indicates how many sites had evidence of positive selection (significant MEME results)\n\n`negative_sites` indicates how many sites had evidence of negative selection (significant FEL results with site-level dN/dS < 1)\n\n`diff_sites` indicates how many sites had significant test results for contrast-FEL, indicating that selection may be acting differently for those sites across clades.","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":42,"value":"//viewof selectedClade = Inputs.select([\"All\", \"DENV1\", \"DENV2\", \"DENV3\", \"DENV4\"], {label: \"Filter by clade\"})\nviewof sites_selections_0 = Inputs.form({\n option1: Inputs.select([\"All\", ...comparison_groups], {label: \"Filter by group\"}),\n option2: Inputs.select([\"All\", ...new Set(summary.map(d => d.gene))], {label: \"Filter by gene\"}),\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":29,"value":"viewof table = Inputs.table(filteredSummary, {\n columns: [\n \"gene\",\n \"sites\",\n \"N\",\n \"T\",\n \"positive_sites\",\n \"negative_sites\",\n \"diff_sites\",\n \"comparison_group\",\n \"group_N\",\n \"group_T\",\n \"group_dN/dS\",\n ],\n header: {\n \"N\": \"# branches\",\n \"T\": \"tree length\",\n \"positive_sites\": \"+ sites\",\n \"negative_sites\": \"- sites\",\n \"diff_sites\": \"diff sites\",\n \"comparison_group\": \"branch set\",\n \"group_N\": \"# branches (group)\",\n \"group_T\": \"tree length (group)\",\n \"group_dN/dS\": \"dN/dS (group)\",\n },\n // width: {\n // \"gene\": 200,\n // },\n format: {\n \"group_T\": (x) => d3.format(\".4f\")(x),\n \"group_dN/dS\": (x) => d3.format(\".4f\")(x),\n \"BUSTED_pval\": (x) => d3.format(\".4\")(x),\n \"BUSTED_omega3\": (x) => d3.format(\".4\")(x),\n \"BUSTED_prop_sites_in_omega3\": (x) => d3.format(\".4\")(x),\n \"RELAX_clade_K\": (x) => d3.format(\".4\")(x),\n \"RELAX_overall_pval\": (x) => d3.format(\".4\")(x),\n }\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":107,"value":"## Gene-wide evidence for positive selection","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":109,"value":"BUSTED (Murrell et al. 2015) was used to search for gene-wide evidence of positive selection for at least one site on at least one branch of the provided gene tree. BUSTED estimates the proportion of sites belonging to each of three rate classes, constrained as ω1 ≤ ω2 ≤ 1 ≤ ω3. The positive class is the rate class where we assign sites with evidence of positive selection, for which we estimate ω3.\n\nA **p-value** < 0.05 indicates that BUSTED detected statistically significant evidence of positive selection. **omega3** is the evolutionary rate ratio estimated from all the sites in the positive class. **% of sites in positive class** indicates how many sites in the original alignment contributed to the evidence for positive selection. BUSTED itself isn't designed to detect *which* sites may be under positive selection; see **[Per-site evidence of selection](https://observablehq.com/d/562fa8ed5fbda59f#cell-146)** instead.","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":120,"value":"viewof keep_significant = Inputs.toggle({label: html`Show only significant results`, value: true, default: true})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":63,"value":"viewof table2 = Inputs.table(BUSTED_results, {\n columns: [\n \"gene\",\n \"BUSTED_pval\",\n \"BUSTED_omega3\",\n \"BUSTED_prop_sites_in_omega3\",\n ],\n header: {\n \"N\": \"internal branches\",\n \"T\": \"tree length\",\n \"BUSTED_pval\": \"P-value\",\n \"BUSTED_omega3\": \"omega3\",\n \"BUSTED_prop_sites_in_omega3\": \"% of sites in positive class\",\n \n },\n format: {\n \"BUSTED_pval\": (x) => d3.format(\".3\")(x),\n \"BUSTED_omega3\": (x) => d3.format(\".4\")(x),\n //\"BUSTED_prop_sites_in_omega3\": (x) => d3.format(\".4\"), //(x * 100),\n }\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":92,"value":"## Relaxed / intensified selection across clades","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":94,"value":"The K-value is the relaxation/intensification parameter from RELAX, plotted here on a log scale. Lower K-values (log(K) < 0) indicate relaxed selection compared to the Reference branches, while higher K-values (log(K) > 0) indicate intensified - but not *necessarily* positive - selection. p < 0.05 indicates that the differences are significant (orange bars).","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":876,"value":"$2e","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":146,"value":"## Per-site evidence of selection","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":148,"value":"### Selection across all branches","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":295,"value":"$2f","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":1006,"value":"viewof q_cut = Inputs.range([0, 1], {value: 0.2, step: 0.05, label: \"Significance threshold (q-value):\"})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":168,"value":"viewof sites_selections_1 = Inputs.form({\n option1: Inputs.select([\"All\", ...comparison_groups], {label: \"Filter by group\"}),\n option2: Inputs.select([\"All\", ...new Set(sites.map(d => d.gene))], {label: \"Filter by gene\"}),\n option3: Inputs.toggle({label: \"Show only significant FEL results\", value: false, default: false}),\n option4: Inputs.toggle({label: \"Show only significant MEME results\", value: false, default: false}),\n option5: Inputs.toggle({label: \"*Show only significant Contrast-FEL results\", value: false, default: false}),\n option6: Inputs.toggle({label: \"Show only significant PRIME results\", value: false, default: false}),\n option7: Inputs.toggle({label: \"*Show intensified positive selection\", value: false, default: false}),\n option8: Inputs.toggle({label: \"*Keep only sites with different majority residues\", value: false, default: false}),\n // option3: Inputs.select([\"All\", ...new Set(sites.map(d => d.site))], {label: \"Filter by site\"}),\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":281,"value":"viewof table3 = Inputs.table(filteredSites, {\n columns: [\n \"gene\",\n \"comparison_group\",\n \"site\",\n \"majority_residue\",\n \"fel_selection\",\n \"meme_marker\",\n \"cfel_marker\",\n \"cfel_beta\",\n \"prime_marker\"\n ],\n header: {\n \"comparison_group\": \"branch set\",\n \"majority_residue\": \"majority residue\",\n \"fel_selection\": \"FEL results\",\n \"meme_marker\": \"MEME pval\",\n \"cfel_marker\": \"C-FEL pval\",\n \"cfel_beta\": \"C-FEL Beta\",\n \"prime_marker\": \"PRIME properties under selection\",\n },\n width: {\n \"gene\": 100,\n \"site\": 30,\n // \"cfel_marker\": 120,\n // \"cfel_beta\": 120,\n \"prime_marker\": 250\n },\n format: {\n \"meme_marker\": (x) => formatMarker(x),\n \"cfel_marker\": (x) => formatMarker(x),\n \"cfel_beta\": (x) => d3.format(\".4\")(x),\n }\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":256,"value":"Site composition indicates the prevalence of each amino acid at a given site. It does not include how many sites had missing data (Codon \"NNN\"), as branches with missing site data could not contribute to the evolutionary analyses and were ignored for the given site.","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":195,"value":"//viewof choose_site = Inputs.select([...new Set(filteredSites.map(d => d.site))], {label: \"Sites in current serotype + gene with significant results\"})\nviewof choose_gene = Inputs.select([...new Set(sites.map(d => d.gene))], {label: \"Filter by gene\"})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":468,"value":"// viewof choose_group = Inputs.select(\n// [...comparison_groups], \n// {label: \"Filter by group\"}\n// )\nviewof choose_group = comparison_groups.size > 0\n ? Inputs.select([...comparison_groups], {value: \"Foreground\", label: \"Filter by group\"})\n : Inputs.select([], {label: \"Filter by group\"})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":459,"value":"viewof choose_site = Inputs.select(\n [...new Set(sites\n .filter(d => (!choose_group || d.comparison_group === choose_group) &&\n (!choose_gene || d.gene === choose_gene))\n .map(d => d.site))], \n {value: 2, label: \"Filter by site\"}\n)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":155,"value":"$30","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":279,"value":"The following chart indicates the most common substitutions of one amino acid for another at a given site. Hyphy analyses use _time reversible_ codon evolutionary models, so these substitutions should NOT be viewed as directional. For example, a substitution listed as Q:R could equally be written as R:Q. The total number of substitutions would reflect the sum of both Q:R and R:Q substitutions. ","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":273,"value":"$31","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":126,"value":"## Data","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":837,"value":"combined_summary = combined_summary_file.csv({typed: true})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":839,"value":"combined_sites = combined_sites_file.csv({typed: true})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":841,"value":"combined_comparison_summary = (\n combined_comparison_summary_file && combined_comparison_summary_file.name.endsWith(\".csv\")\n)\n ? await combined_comparison_summary_file.csv({ typed: true })\n : []","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":843,"value":"combined_comparison_sites = (\n combined_comparison_sites_file && combined_comparison_sites_file.name.endsWith(\".csv\")\n)\n ? await combined_comparison_sites_file.csv({ typed: true })\n : []","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":621,"value":"comparison_groups = {\n const groups = new Set(summary.map(d => d.comparison_group));\n groups.delete(\"background\");\n return groups\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":20,"value":"//summary = FileAttachment(\"capsid_protein_C_summary@2.csv\").csv()\nsummary = mergePrimaryIntoSecondary(combined_summary, combined_comparison_summary) // combineJsonData(summary_urls)\n .filter(d => d.comparison_group != 'background') // remove nuisance background (leaf) branches\n .map(r => (r.comparison_group && r.comparison_group.toLowerCase() === \"reference\") \n ? { ...r, RELAX_K: 1 } // reset RELAX K of 'Reference' to 1 (RELAX default for reference branches)\n : r\n )","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":54,"value":"filteredSummary = summary\n .filter(d => sites_selections_0.option1 === \"All\" || d.comparison_group === sites_selections_0.option1)\n .filter(d => sites_selections_0.option2 === \"All\" || d.gene === sites_selections_0.option2);","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":130,"value":"BUSTED_results = keep_significant\n ? combined_summary.filter(d => d.BUSTED_pval < 0.05)\n : combined_summary;","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":17,"value":"// sites = combineJsonData(sites_urls)\n// merge on gene AND site number\nsites = mergePrimaryIntoSecondary(combined_sites, combined_comparison_sites, {key: siteKey})\n .filter(d => d.comparison_group != 'background')","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":182,"value":"filteredSites = sites\n .filter(d => sites_selections_1.option1 === \"All\" || d.comparison_group === sites_selections_1.option1)\n .filter(d => sites_selections_1.option2 === \"All\" || d.gene === sites_selections_1.option2)\n .filter(d => !sites_selections_1.option3 || d.fel_qval < q_cut)\n .filter(d => !sites_selections_1.option4 || d.meme_marker < 0.05) // change to meme_qval < q_cut\n .filter(d => !sites_selections_1.option5 || d.cfel_marker < q_cut)\n .filter(d => !sites_selections_1.option6 || d.prime_qval < q_cut)\n .filter(d => !sites_selections_1.option7 || intensifiedSites.has(siteKey(d)))\n .filter(d => !sites_selections_1.option8 || majorityDiffPairs.has(siteKey(d)));\n //.filter(d => sites_selections_1.option3 === \"All\" || d.site === sites_selections_1.option3);","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":723,"value":"// Set of gene+site pairs that (1) have a non-'-' meme_marker,\n// (2) have a non-'-' cfel_marker, and (3) Foreground cfel_beta > Reference cfel_beta\nintensifiedSites = {\n const byKey = d3.group(sites, siteKey); // Map key -> rows for that site (both groups)\n const out = new Set();\n\n for (const [key, rows] of byKey) {\n // must have at least one non-'-' meme marker at this site\n const hasMeme = rows.some(r => r.meme_marker !== \"-\");\n if (!hasMeme) continue;\n\n const hasCFEL = rows.some(r => r.cfel_marker !== \"-\");\n if (!hasCFEL) continue;\n\n // find Foreground and Reference rows\n const fg = rows.find(r => r.comparison_group === \"Foreground\");\n const bg = rows.find(r => r.comparison_group === \"Reference\");\n if (!fg || !bg) continue;\n\n // numeric compare (coerce just in case)\n const fgBeta = +fg.cfel_beta;\n const bgBeta = +bg.cfel_beta;\n if (Number.isFinite(fgBeta) && Number.isFinite(bgBeta) && fgBeta > bgBeta) {\n out.add(key);\n }\n }\n return out;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":787,"value":"// Set of (gene,site) pairs whose majority_residue differs FG vs Reference\nmajorityDiffPairs = {\n const byKey = d3.group(sites, siteKey); // Map\n const out = new Set();\n\n for (const [key, rows] of byKey) {\n const fg = rows.find(r => r.comparison_group === \"Foreground\");\n const ref = rows.find(r => r.comparison_group === \"Reference\" || r.comparison_group === \"Background\");\n if (!fg || !ref) continue;\n\n const a = residue(ref.majority_residue);\n const b = residue(fg.majority_residue);\n if (a && b && a !== b) out.add(key);\n }\n return out;\n};","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":207,"value":"//focal_site = sites.filter(d => d.site === choose_site && d.comparison_group === choose_group && d.gene === choose_gene)\nfocal_site = sites.filter(d => \n d.site === choose_site &&\n d.gene === choose_gene &&\n (!choose_group || d.comparison_group === choose_group)\n)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":213,"value":"focal_site_composition = {\n const raw = focal_site?.[0]?.composition;\n if (!raw || raw === \"-\") return [];\n\n return raw\n .split(/\\s*,\\s*/) // e.g., [\"E:9\", \"T:1\"]\n .filter(Boolean)\n .map(tok => {\n const i = tok.lastIndexOf(\":\");\n if (i <= 0) return null; // malformed\n const category = tok.slice(0, i).trim();\n const size = Number(tok.slice(i + 1));\n return Number.isFinite(size) ? { category, size } : null;\n })\n .filter(Boolean);\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":260,"value":"focal_site_substitutions = {\n const s = focal_site?.[0]?.substitutions;\n if (!s || s === \"-\") return [];\n\n return s\n .split(/\\s*,\\s*/) // [\"E:E:9\", \"E:T:1\"]\n .filter(Boolean)\n .map(entry => {\n const i = entry.lastIndexOf(\":\");\n const category = entry.slice(0, i).trim(); // \"E:E\"\n const size = Number(entry.slice(i + 1)); // 9\n return { category, size };\n })\n .filter(d => Number.isFinite(d.size));\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":7,"value":"### Functions & Imports","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","database":"$undefined"},{"id":9,"value":"vegalite = require(\"vega-embed\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":707,"value":"siteKey = d => `${d.gene}::${d.site}`;","pinned":true,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":789,"value":"// Helper to normalize residue letters (and ignore missing/unknown)\nresidue = x => {\n if (x == null) return null;\n const r = String(x).trim().toUpperCase();\n return (r === \"-\" || r === \"?\") ? null : r;\n};","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":739,"value":"formatMarker = {\n const fmt = d3.format(\".4~g\");\n return x => (x == null || x === \"-\" ? x : fmt(+x));\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":584,"value":"// Merge \"primary\" fields into each row of \"secondary\" (if it exists).\n// Matching key defaults to `gene`, but you can pass a custom key fn.\nfunction mergePrimaryIntoSecondary(primary, secondary, { key = d => d.gene } = {}) {\n // If there's no secondary, return primary as-is\n if (!Array.isArray(secondary) || secondary.length === 0) return primary;\n\n // Index primary by the chosen key\n const idx = new Map(primary.map(p => [key(p), p]));\n\n // For each secondary row, overlay its fields on top of the matching primary row\n return secondary.map(s => {\n const p = idx.get(key(s));\n // If no match, keep the secondary row; otherwise merge (secondary wins on conflicts)\n return p ? { ...p, ...s } : s;\n });\n}\n\n// --- Usage ---\n// Basic (match by gene)\n// const merged = mergePrimaryIntoSecondary(firstArray, secondArray);\n\n// // Stricter match (gene + fasta_id)\n// const mergedStrict = mergePrimaryIntoSecondary(\n// firstArray,\n// secondArray,\n// { key: d => `${d.gene}::${d.fasta_id}` }\n// );","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":403,"value":"async function combineJsonData(urls) {\n const combinedData = [];\n\n // Fetch and combine data from each URL\n for (const url of urls) {\n try {\n const data = await d3.csv(url, d3.autoType);\n combinedData.push(...data); // Combine arrays\n } catch (error) {\n console.error(`Failed to load data from ${url}:`, error);\n }\n }\n\n return combinedData;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":322,"value":"function pivotCfelMarker(data) {\n return data.map(item => {\n const newItem = { ...item };\n if (item.cfel_marker && item.cfel_marker !== '-') {\n const markers = item.cfel_marker.split(',').map(m => m.trim());\n markers.forEach(marker => {\n const [key, value] = marker.split(':').map(s => s.trim());\n newItem[key] = value;\n });\n }\n delete newItem.cfel_marker;\n return newItem;\n });\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":956,"value":"// determines the extent of the log(RELAX_K values) on the RELAX chart, so we can add a little padding to the extent. This way, K and p values aren't clipped by the x-axis\nextent = d3.extent(summary.map(d => Math.log(d.RELAX_K)))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W32","stdlib":"1","titling":"auto"}} 2d:"$2c:body"