28:null 29:["$","title",null,{"children":"CAPHEINE H5N1 Results | Hannah Verdonk | Observable"}] 2e:Tdeb,vegalite({ "$schema": "https://vega.github.io/schema/vega-lite/v5.json", "data": { "values": combined_summary }, "width": 900, // 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:T9f6,**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.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:T810,/** * Apply per-gene BH to MEME and FEL. * * Adds: * - meme_q_bh * - fel_q_bh * * Options: * - geneKey: field to group by (default "gene") * - methods: array of { pField, qField } definitions * - filter: optional predicate to select rows (e.g. Foreground only) */ function addPerGeneBH(sites, { geneKey = "gene", methods = [ { pField: "meme_marker", qField: "meme_qval" }, { pField: "fel_pval", qField: "fel_qval" } // <-- change if your FEL p-value field is named differently ], filter = null } = {}) { // Work on a shallow copy so we don't mutate caller data const out = sites.map(d => ({ ...d })); // Helper: parse p-values that might be "-", null, etc. const toP = (v) => { if (v === null || v === undefined) return null; if (typeof v === "string") { const t = v.trim(); if (t === "" || t === "-") return null; const n = Number(t); return Number.isFinite(n) ? n : null; } if (typeof v === "number") return Number.isFinite(v) ? v : null; const n = Number(v); return Number.isFinite(n) ? n : null; }; // Group row indices by gene const geneToIndices = new Map(); for (let i = 0; i < out.length; i++) { if (filter && !filter(out[i])) continue; const g = out[i][geneKey]; if (!g) continue; if (!geneToIndices.has(g)) geneToIndices.set(g, []); geneToIndices.get(g).push(i); } // For each gene, compute BH per method for (const [gene, indices] of geneToIndices.entries()) { for (const { pField, qField } of methods) { const indexedPvals = []; for (const idx of indices) { const p = toP(out[idx][pField]); // Accept only p in [0,1] if (p === null) continue; if (p < 0 || p > 1) continue; indexedPvals.push({ idx, p }); } const qMap = benjaminiHochbergForIndices(indexedPvals); // Write q-values back; leave missing as null for (const idx of indices) { out[idx][qField] = qMap.has(idx) ? qMap.get(idx) : null; } } } return out; }33:T665,// pre-filter sites for plotting filtered = sites .filter(d => d.comparison_group === "Foreground") //"Foreground") .map(d => { const site_num = +d.site; if (!Number.isFinite(site_num)) return null; //const gene_label = d.short_gene ?? d.gene; const cfel_raw = +d.cfel_marker; const fel_raw = d.fel_qval != null && d.fel_qval !== "" ? +d.fel_qval : +d.fel_pval; const meme_raw = d.meme_qval != null && d.meme_qval !== "" ? +d.meme_qval : +d.meme_marker; const prime_raw = +d.prime_qval; const fel_beta_raw = +d.fel_beta; const fel_alpha_raw = +d.fel_alpha; const fel_beta_eps = (Number.isFinite(fel_beta_raw) ? fel_beta_raw : 0) + epsilon; const fel_alpha_eps = (Number.isFinite(fel_alpha_raw) ? fel_alpha_raw : 0) + epsilon; const fel_effect = Math.log(fel_beta_eps / fel_alpha_eps); const method_rows = [ { method: "CFEL", raw_val: cfel_raw }, { method: "FEL", raw_val: fel_raw }, { method: "MEME", raw_val: meme_raw }, { method: "PRIME", raw_val: prime_raw } ]; return method_rows.map(m => ({ ...d, //gene_label, site_num, site_start: site_num - 0.5, // coords for plotting on graph site_end: site_num + 0.5, // coords for plotting on graph method: m.method, raw_val: m.raw_val, significant: Number.isFinite(m.raw_val) && m.raw_val <= q_cut, // ONLY WORKS IF ALL METHODS HAVE QVAL fel_alpha_raw, fel_beta_raw, fel_effect })); }) .filter(Boolean) .flat();34:T1b24,function geneSpec(gene, values, { showYAxis = true, showLegend = false } = {}) { const maxSite = d3.max(values, d => d.site_num); const width = panelWidth(maxSite); return { // "title": { // "text": gene, // "anchor": "end", // "fontSize": 16, // "fontWeight": "bold", // "offset": 8 // }, "width": width, "height": 160, "data": { "values": values }, "layer": [ // ------------------------------ Per-Gene Title Text ------------------------------ { "data": { "values": [{ "label": gene, "x": maxSite / 2 }] }, "mark": { "type": "text", "align": "center", "baseline": "bottom", "fontSize": 16, "fontWeight": "bold", "dy": -6 }, "encoding": { "x": { "field": "x", "type": "quantitative", "scale": { "domain": [0.5, maxSite + 0.5], "nice": false } }, "y": { "value": 0 }, "text": { "field": "label" } } }, // ---------------- non-FEL methods: colored by method if significant ---------------- { "transform": [ { "filter": "datum.significant && datum.method !== 'FEL'" } ], "mark": { "type": "tick", "orient": "vertical", "thickness": 2.2, "size": 40 }, "encoding": { "x": { "field": "site_num", "type": "quantitative", "scale": { "domain": [0.5, maxSite + 0.5], "nice": false }, "axis": { "title": "Site", "titleAnchor": "end", "titlePadding": 10 } }, "y": { "field": "method", "type": "nominal", "sort": ["FEL", "MEME", "CFEL", "PRIME",], "scale": { "domain": ["FEL", "MEME", "CFEL", "PRIME",], "paddingInner": 0.05, "paddingOuter": 0 }, "axis": showYAxis ? { "title": "Analysis" } : { "title": null, "labels": false, "ticks": false, "domain": true, "tickSize": 6 } }, "color": { "field": "method", "type": "nominal", "scale": { "domain": ["FEL", "MEME", "PRIME", "CFEL"], // "FEL+", "range": ["#56B4E9", "#E69F00", "#009E73", "#D55E00"] }, "legend": null }, "tooltip": [ { "field": "gene", "title": "Gene" }, { "field": "site_num", "title": "Site" }, { "field": "method", "title": "Method" }, { "field": "raw_val", "title": "q (or p fallback)" } ] } }, // ---------------- FEL: colored by effect size if significant ---------------- { "transform": [ { "filter": "datum.significant && datum.method === 'FEL' && isFinite(datum.fel_effect)" } ], "mark": { "type": "rect" }, "encoding": { "x": { "field": "site_start", "type": "quantitative", "scale": { "domain": [0.5, maxSite + 0.5], "nice": false }, "axis": { "title": "Site", "titleAnchor": "end", "titlePadding": 10 } }, "x2": { "field": "site_end" }, "y": { "field": "method", "type": "nominal", "sort": ["FEL", "MEME", "CFEL", "PRIME",], "scale": { "domain": ["FEL", "MEME", "CFEL", "PRIME",], "paddingInner": 0.05, "paddingOuter": 0 } }, "color": { "field": "fel_effect", "type": "quantitative", "scale": { "domain": [-maxAbsFelEffect, maxAbsFelEffect], "scheme": "redblue", "domainMid": 0, "reverse": true }, "legend": showLegend ? { "title": "FEL log(ω)", // "FEL log((β+ε)/(α+ε))", "orient": "left", "direction": "vertical", "gradientLength": 120 } : null // "legend": { // "title": "FEL log((β+ε)/(α+ε))", // "orient": "left", // "direction": "vertical", // "gradientLength": 120 // } }, "tooltip": [ { "field": "gene", "title": "Gene" }, { "field": "site_num", "title": "Site" }, { "field": "method", "title": "Method" }, { "field": "raw_val", "title": "q (or p fallback)" }, { "field": "fel_alpha_raw", "title": "FEL alpha" }, { "field": "fel_beta_raw", "title": "FEL beta" }, { "field": "fel_effect", "title": "log((β+ε)/(α+ε))", "format": ".3f" } ] } }, // { // "transform": [ // { "filter": "datum.significant && datum.method === 'FEL' && datum.fel_effect > 1" }, // { "calculate": "'FEL+'", "as": "method_plus" } // ], // "mark": { "type": "tick", "orient": "vertical", "thickness": 2.2, "size": 40 }, // "encoding": { // "x": { // "field": "site_start", // "type": "quantitative", // "scale": { "domain": [0.5, maxSite + 0.5], "nice": false } // }, // "x2": { "field": "site_end" }, // "y": { // "field": "method_plus", // "type": "nominal", // "sort": methodOrder, // "scale": { // "domain": methodOrder, // "paddingInner": 0.05, // "paddingOuter": 0 // }, // "axis": showYAxis ? { // "title": "Analysis", // "labels": true, // "ticks": true, // "domain": true // } : null // }, // "color": { "value": "#CC3311" }, // "tooltip": [ // { "field": "gene_label", "title": "Gene" }, // { "field": "site_num", "title": "Site" }, // { "value": "FEL+", "title": "Track" }, // { "field": "raw_val", "title": "FEL q" }, // { "field": "fel_effect", "title": "FEL log((β+ε)/(α+ε))", "format": ".3f" } // ] // } // } ] }; } 35:[] 2c:{"id":"09d06c638d806915","slug":"capheine-h5n1-results","title":"CAPHEINE H5N1 Results","thumbnail":"a44807cc25dc481c6ab1e39737e052c7c44dc05dc7f644831ff6f82c0bb5d1cc","access_level":"public","create_time":"$D2026-08-28T20:50:11.486Z","update_time":"$D2026-09-17T20:02:48.406Z","view_time":null,"first_public_time":"$D2026-08-28T20:50:11.486Z","has_unpublished_edits":null,"fork":{"id":"df8bfd20604373de","slug":"drhip-result-viewer","title":"DRHIP result viewer","owner":{"id":"7ada8a89538f700a","login":"hyphy","name":"ACME Computational Molecular Evolution Group"}},"stars":0,"starred":false,"owner":{"type":"team","id":"693c07e84b050d0c","login":"hverdonk","name":"Hannah Verdonk","avatar_url":"https://avatars-next.observableusercontent.com/e3886da95f22b11c72bd727c03af4612d64ebac68ae9c2a7a87ae5122c0d27f3","tier":"starter_2024"},"version":1108,"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":"team","id":"693c07e84b050d0c","login":"hverdonk","name":"Hannah Verdonk","avatar_url":"https://avatars-next.observableusercontent.com/e3886da95f22b11c72bd727c03af4612d64ebac68ae9c2a7a87ae5122c0d27f3","tier":"starter_2024","role":"editor"},{"type":"individual","id":"d7d30d996f330cf9","login":"hverdonk","name":"Hannah Verdonk","avatar_url":"https://avatars-next.observableusercontent.com/e3886da95f22b11c72bd727c03af4612d64ebac68ae9c2a7a87ae5122c0d27f3","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":0,"collection":null,"type":"notebook","body":{"title":"CAPHEINE H5N1 Results","cells":[{"id":0,"value":"# CAPHEINE H5N1 Results","pinned":false,"hidden":false,"mode":"md","author":"$undefined","summary":"$undefined","output":"","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":1044,"value":"vegalite({\n \"$schema\": \"https://vega.github.io/schema/vega-lite/v5.json\",\n \"vconcat\": rows,\n \"spacing\": 14,\n \"resolve\": {\n \"scale\": {\n \"color\": \"shared\"\n }\n },\n \"config\": {\n \"font\": \"Helvetica\",\n \"axis\": {\n \"labelFontSize\": 14,\n \"titleFontSize\": 13,\n \"tickSize\": 3,\n \"domain\": true,\n \"grid\": false\n },\n \"view\": { \"stroke\": null }\n }\n }\n)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","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 = {\n const f = await d3.csv(\"https://raw.githubusercontent.com/hverdonk/CAPHEINE-H5N1/refs/heads/main/drhip/combined_summary.csv\", d3.autoType);\n return f.filter(d => (d.gene != \"PA-X\" && d.gene != \"PB1-F2\"));\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":839,"value":"combined_sites = {\n const f = await d3.csv(\"https://raw.githubusercontent.com/hverdonk/CAPHEINE-H5N1/refs/heads/main/drhip/combined_sites.csv\", d3.autoType);\n return f.filter(d => (d.gene != \"PA-X\" && d.gene != \"PB1-F2\"));\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":841,"value":"combined_comparison_summary = {\n const f = await d3.csv(\"https://raw.githubusercontent.com/hverdonk/CAPHEINE-H5N1/refs/heads/main/drhip/combined_comparison_summary.csv\", d3.autoType);\n return f.filter(d => (d.gene != \"PA-X\" && d.gene != \"PB1-F2\"));\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":843,"value":"combined_comparison_sites = {\n const f = await d3.csv(\"https://raw.githubusercontent.com/hverdonk/CAPHEINE-H5N1/refs/heads/main/drhip/combined_comparison_site.csv\", d3.autoType);\n return f.filter(d => (d.gene != \"PA-X\" && d.gene != \"PB1-F2\"));\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":1073,"value":"sitesWithBH = addPerGeneBH(sites, {\n geneKey: \"gene\",\n methods: [\n { pField: \"meme_marker\", qField: \"meme_qval\" },\n { pField: \"fel_pval\", qField: \"fel_qval\" }\n ],\n});","pinned":false,"hidden":false,"mode":"ts","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_pval < 0.05)\n .filter(d => !sites_selections_1.option4 || d.meme_marker < 0.05)\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":1070,"value":"$32","pinned":false,"hidden":false,"mode":"ts","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1071,"value":"/**\n * Benjamini–Hochberg (BH) correction.\n * Input: array of objects that include a numeric p-value field.\n * Output: Map from row index -> q-value (BH-adjusted p).\n */\nfunction benjaminiHochbergForIndices(indexedPvals) {\n // indexedPvals: [{ idx, p }]\n const m = indexedPvals.length;\n if (m === 0) return new Map();\n\n // Sort ascending by p\n const sorted = [...indexedPvals].sort((a, b) => a.p - b.p);\n\n // Compute raw BH q = p * m / rank\n const qRaw = new Array(m);\n for (let i = 0; i < m; i++) {\n const rank = i + 1; // 1-based\n qRaw[i] = (sorted[i].p * m) / rank;\n }\n\n // Enforce monotonicity from the end: q_i = min(q_i, q_{i+1})\n const qMon = new Array(m);\n qMon[m - 1] = Math.min(qRaw[m - 1], 1);\n for (let i = m - 2; i >= 0; i--) {\n qMon[i] = Math.min(qRaw[i], qMon[i + 1], 1);\n }\n\n // Map back to original indices\n const out = new Map();\n for (let i = 0; i < m; i++) {\n out.set(sorted[i].idx, qMon[i]);\n }\n return out;\n}","pinned":false,"hidden":false,"mode":"ts","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"},{"id":1056,"value":"epsilon = 1e-12","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1054,"value":"$33","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1052,"value":"// group rows by gene\nbyGene = d3.group(filtered, d => d.gene);","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1050,"value":"customRows = [\n [\"PB1\", \"NA\"],\n [\"PB2\", \"NP\"],\n [\"PA\", \"HA\", ],\n [\"NEP\", \"M2\", \"M1\", \"NS1\"]\n];","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1046,"value":"rows = customRows.map((rowGenes, rowIndex) => ({\n hconcat: rowGenes\n .filter(gene => byGene.has(gene))\n .map((gene, colIndex) =>\n geneSpec(gene, byGene.get(gene), {\n showYAxis: colIndex === 0,\n showLegend: rowIndex === 0 && colIndex === 0\n })\n ),\n spacing: 24\n }))\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1064,"value":"// used to calculate global FEL color scale\nmaxAbsFelEffect = {\n const felEffects = filtered\n .filter(d => d.method === \"FEL\" && d.significant && Number.isFinite(d.fel_effect))\n .map(d => d.fel_effect);\n \n return d3.max(felEffects, d => Math.abs(d)) ?? 1;\n}\n\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1068,"value":"widthScale = {\n const geneSizes = Array.from(byGene, ([gene, values]) => ({\n gene,\n nSites: d3.max(values, d => d.site_num)\n }));\n\n const minSites = d3.min(geneSizes, d => d.nSites);\n const maxSites = d3.max(geneSizes, d => d.nSites);\n\n return d3.scaleSqrt()\n .domain([minSites, maxSites])\n .range([160, 520]); \n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1066,"value":"// Width function. Adjust the multiplier/caps to taste. nSites doesn't actually do anything here.\nfunction panelWidth(nSites) {\n //return Math.max(140, Math.min(700, nSites * 3.2));\n return Math.round(widthScale(nSites));\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1062,"value":"$34","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1048,"value":"genePanels = {\n const geneEntries = Array.from(byGene);\n\n return geneEntries.map(([gene, values], i) => {\n const col = i % 3;\n return geneSpec(gene, values, {\n showYAxis: col === 0,\n showLegend: i === 0\n });\n });\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"}],"files":"$W35","stdlib":"1","titling":"auto"}} 2d:"$2c:body"