29:null 2a:["$","title",null,{"children":"Visualizing BUSTED-PH analysis results | Avery Selberg | Observable"}] 2f:T4ad,md` 1. ** ID ** : gene / alignment ID 2. ** Sequences ** : number of sequences in the aligment 3. ** Sites ** : number of sites in the alignment 3. ** L (FG) / L (BG) ** : tree length under the unconstrained model (FG branches / BG branches) 4. ** FG p-value ** : p-value of the BUSTED test on foreground branches; if significant (bolded), implies that there are detectable differences in selection **intensity** between ciliate groups. 5. ** BG p-value ** : p-value of the BUSTED test on background branches; if significant (bolded), implies that there are detectable differences in selection **intensity** between ciliate groups. 6. ** Sel. sites (FG/BG) ** : Sites that are showing support for positive selection in the Foreground/Background branches with the log-evidence ratio of 5 or higher. 7. ** DIFF p-value ** : p-value of the tests for differences in selection between FG and BG; if significant (bolded), implies that there are detectable differences in selection **regimes**. 8 ** Selection associated with phenotype ** : if Yes, then there is selection on **foreground**, no selection on **background**, and **differences** in selection between the two groups. `30:T8ed,displayed_tree = { let tree = tree_object[0]; var t = tree.render("#tree_container", { //height:1024, width:400, '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)) { return 3; } return 1; }; let nodeStates = null; t.node_bubble_size = bs; t.phylotree.node_bubble_size = bs; t.options["draw-size-bubbles"] = true; t.font_size = 10; t.spacing_x (10); t.phylotree.branch_length_accessor = n=> { //console.log (summary_json[segment_id]); let bl = tree_object[2][n.data.name]; /*if (bl > 0.25) { n.data.censored_bl = true; return 0.25; } else { n.data.censored_bl = false; }*/ return bl; } t.node_label ((n)=>{ return n.data.name.split ("_").slice (0,3).join (" "); //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); t.style_nodes ((e,n) => { e.selectAll ("title").data ([n.data.name]).join ("title").text ((d)=>d); if (n.data.name in tree_object[1]) { let binfo = tree_object[1][n.data.name]; e.selectAll ("circle").style ('fill', species_colors (binfo)); e.selectAll ("text").style ('fill', species_colors (binfo)); } //} }); t.style_edges ((e,d) => { if (d.target.data.name in tree_object[1]) { let binfo = tree_object[1][d.target.data.name]; console.log (e); e.style ('stroke', species_colors (binfo)); e.style ('stroke-width', 2); } }); t.placenodes(); t.update(); return t; }31:T710,plotPie = (values, id, skipK)=>{ skipK = true; return vegalite ({ "columns": 3, "hconcat": _.map (_.groupBy (values, (d)=>d.type), (d, type) => { let is_dnds = "omega" in d[0]; let title = (is_dnds ? "ω for " : "") + type + (type == baseline_by_segment[id] ? " (reference)" : ""); let rate = is_dnds ? "omega" : "rate"; return { "title" : title, "data" : {"values" : d}, "encoding": { "theta": {"field": "proportion", "type": "quantitative", "stack": true}, "color": {"field": rate, "type": "quantitative", "title" : "ω", "legend" : {"titleFontSize" : 14, "orient" : "top"}, "scale": { "domain" : [0,1,4], "type" : "sqrt", "clamp" : true, "scheme": "sinebow" } }}, "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": rate, "type": "nominal", "format" : ".2f", "condition": {"test": "datum['" + rate +"'] <= 1.", "value": ""}, } } }, { "mark": {"type": "text", "outerRadius" : 20, "fontSize" : 18}, "encoding": { "color" : {"value" : "black"}, "text": { "field": "proportion", "type": "nominal", "format" : ".2p", "condition": {"test": "datum['" + rate +"'] <= 1.", "value": ""} } } } ]} }), "view": {"stroke": null} }) }32:T400,decision_tree = { if (_.get (result_json,[gene_id,"tests","FG","p-value"]) > .05) { return "no evidence of selection associated with the " + phenotype + " phenotype" } else { if (_.get (result_json,[gene_id,"tests","BG","p-value"]) > 0.05) { if (_.get (result_json,[gene_id,"tests","DIFF","p-value"]) > .05) { return "evidence of selection associated with the " + phenotype + " phenotype, but no statistically significant difference between the test and background selective regimes" } else { return "strong evidence of selection is associated with the " + phenotype + " phenotype" } } else { if (_.get (result_json,[gene_id,"tests","DIFF","p-value"]) > .05) { return "selection is occurring with or without the " + phenotype + " phenotype, but with similar selective regimes" } else { return "selection is occurring with or without the " + phenotype + " phenotype, but with different selective regimes" } } } }33:T4e4,averys_table = { const new_table = []; const keys = ["ID", "Sequences", "Sites", "FG p-value", "BG p-value", "DIFF p-value", "PH"]; for (let j = 0; j < table_data.length; j++){ for (let i = 0; i < keys.length; i++){ new_table[j] = {"ID": table_data[j]["ID"], "Sequences": table_data[j]["Sequences"], "Sites": table_data[j]["Sites"], "L BG": table_data[j]["L"]["background"], "L FG": table_data[j]["L"]["test"], "FG p-value": table_data[j]["FG p-value"], "BG p-value": table_data[j]["BG p-value"], "SelSites FG": table_data[j]["SelSites"]["0"], "SelSites BG": table_data[j]["SelSites"]["1"], "DIFF p-value": table_data[j]["DIFF p-value"], "Selection associated with phenotype": table_data[j]["PH"] }; } } for (let k = 0; k < new_table.length; k++){ if(new_table[k]["Selection associated with phenotype"] == false){ new_table[k]["Selection associated with phenotype"] = "No"} else { new_table[k]["Selection associated with phenotype"] = "Yes"} } return new_table } 34:T41eb,marine_genes = ['GNAT3','OR6B1','CALHM1','SSTR4','COL9A2','FGF11','HMGCS2','CLSTN2','PERP','S100A5','OR51M1','OR51I2','KRTAP10-3','OR51I1','TRIM29','CNGA4','PAPL','GRIN2C','GRIN3B','FER1L4','KRT80','GLRA4','RAG1','DSP','CAB39L','OR52D1','STOML3','PIK3R5','RNF222','OR10Z1','TGM3','SFTPB','CAPN12','OR6A2','OR5C1','OR51A7','CAPN8','OR2A25','ANKDD1A','OR10P1','PDE1C','METTL11B','TACR2','ANO2','KRT17','OTOF','HECTD1','S100Z','OR8S1','TCRBV9S1A1T','OR4E2','EPHA8','OR5A1','NEB','OR51G2','PNPLA4','KRT6C','TPRG1','PKNOX2','MYOF','OR51G1','CNGA2','ALDH3B2','LHFP','OR10A4','TMEM235','HNRNPF','MYH7','ASIC4','TRPV3','SERPINB12','NT5DC3','OR6K6','CYP3A43','PRSS3','LGALS7B','FMO1','OR4S1','ANXA6','SVOPL','OR51L1','OR51V1','UNC5C','TKT','GSN','CDH12','SBSN','TMEM25','OR52K1','C1QTNF8','OR10G3','ULK2','CHAD','SAP18','CSDE1','TUBA1B','GSTM4','PGAM2','VSIG8','KRTAP12-3','SAP30BP','OMP','AQP10','FHL5','TAGLN2','CSMD3','CALML3','OR1I1','SLC22A20','COL9A3','OR9I1','ZDHHC23','OR9Q1','FEZ1','MADD','PRSS1','SLC47A1','MFRP','OR10S1','IGSF22','HTR1D','MARCH9','NPSR1','SLC2A4','CAPNS2','KRTAP8-1','MIP','OR2F2','NT5C2','SLC12A9','PLEC','FAM179A','OLFR959','KRT35','PCDHB1','MYO5A','CALHM3','TAAR3','HAO1','ANKRD23','APBB1','JUP','CLK4','OR2W3','VTI1B','PIGU','OR2F1','KRT40','FMO2','CMTM1','CADPS2','OR5AU1','VPS36','ZEB1','PTPRM','TMPRSS13','HTR6','KRTAP10-12','NT5DC4','ZBED2','HIGD1B','CHI3L2','MAGEA4','GSTM1','SLITRK5','PAH','CPNE2','OR6X1','STRC','KRT85','YRDC','AMY1A','OR6N1','KCNK18','PLEKHD1','EPHA4','GABRP','KRTAP10-8','POR','ACTR1B','HTR3E','ANO4','PRM1','ALOX12B','TRIM54','COL11A2','KRT84','GRAMD1A','PCDH9','OR2B11','AVIL','OR52M1','OR1J4','ABCC8','ENO3','OR7C2','TCRBV22S1A2N1T','OR5E1P','OR6K2','OR10Q1','BAC05914.1','PLD2','GRM8','FER1L6','ACSL3','ANKFN1','KRT28','OR5H2','OR8D4','PAMR1','NDRG4','WIPI1','PRG4','OR5L2','KRTAP6-2','PTK2B','PDK2','FSD1','PKP1','MNT','KRT10','MCHR2','OR6B2','ASB15','SGSM2','BPIFA1','KLK5','CRP','KRTAP4-5','SLC6A14','OR56A5','COX17','DUSP22','PTAFR','ANXA11','MYH4','KRTAP6-1','OR2C3','ARSI','FUS','OR10AD1','IL1F10','PFN1','OTUD6A','TESC','CKM','MKS1','MYH1','MAMDC2','MYO5C','MYH7B','ATP2A1','GRAMD1B','ABHD16A','AQP5','CELA1','OR10J3','SLC6A18','ZFHX4','ABHD1','NKD2','MYH2','SCN1A','RTP2','MUC13','OR10A6','PPP2R1B','OR10A2','TAAR2','MYBPC3','PTK6','WBP2','UNC13C','SCN5A','KCNIP1','HSPA6','PRSS3P2','TRIM63','INSIG2','OR8D2','ATP2B3','NRAP','NOL12','CKMT1A','LCE5A','KLK11','OR52Z1','CKMT1B','GRIA4','NPY6R','MYBPHL','CACNA1E','KRT15','IST1','OR4D1','TCRBV1S1A1N1','OR3A1','ANXA7','PADI4','PLXNA4','CCDC85A','GUCY1B2','OR5R1','OR5K2','DYRK1B','CCT3','ZBTB8OS','OR9Q2','KRTAP9-2','LCE2A','OR10D3','CASQ1','ALOX5','TAAR9','KIAA1147','KRT32','OR4K14','DTNB','KRTAP10-9','FLOT2','TRIM47','BPIFB4','MB','OR10G2','LCE3B','MTNR1B','TGM5','OR51B2','CDH23','DFNB31','ESYT3','OR7D2','ATP5B','FAM229B','MUC2','VANGL1','PPARGC1A','OR56A3','OPN1SW','ABHD17A','CLDN11','ARNT2','YBX2','TMEM132E','OR52E6','AGO1','MYH13','OR5AP2','OR8B3','SYT7','FMO9P','RAPSN','PAPSS2','CALB2','BOLL','ACTBL2','SCIN','KRTAP20-1','DPH3','ACSL1','TUFM','OR2V1','OR52B4','LCE1D','KREMEN1','MYH8','OR9K2','GRIK1','VIM','OR10T2','OR51Q1','OR6S1','VARS','MEAF6','UBB','AK074396','PANX3','RIT2','IGLON5','OR6Y1','DSC1','FAM89A','RYR1','SLC35F4','LYRM5','CAV1','SIDT1','MYBPC2','OR6N2','OR2B6','TSACC','OR5M9','KRTAP10-1','CCL14','PPIP5K1','KLHDC8A','KRT36','KCNA1','CYP2S1','CARTPT','IPO4','SLC25A13','OR52E8','UCMA','TIAM1','OR51A4','TMEM145','LHX2','PYGM','MEGF11','COL9A1','KRT13','CHRNA9','OR56B1','LCE3C','CIB2','CELA2A','ATP13A5','NDRG1','OR4F17','RNF167','MYO1D','SV2C','HIST1H2AA','PSAPL1','KRT31','CDSN','COX15','VIL1','NCKAP1','ARPP21','IFT88','KCNJ11','KCNJ10','TGM6','SLURP1','ZNF593','HIP1R','GNAT1','GBA3','CAMSAP3','SMPX','ARSH','CORO6','SPOCK1','GLRA3','LMBR1L','KRT27','SPEF1','OR51B4','TNNI1','RYR3','OR4L1','SDCBP2','UNC80','GRIK2','CALCOCO1','SLC12A2','SHANK1','OR52E2','CHCHD10','GRASP','KLK12','MAPK13','RBM22','GOLT1A','HTRA4','RNF182','SMYD1','ALDOB','ARPC1A','LINGO4','OR1K1','ENDOG','SULT6B1','IL20','MYL2','LACTB2','INMT','ABCF3','OR52N4','GSTA3','OR12D3','FLNC','PPP1R3E','TCRA','DIP2C','FEM1C','GHITM','UNC13D','LIMD2','DDX1','KRT16','RUSC2','TENM3','RNASE2','OR5H1','OR4F4','OR1J1','OR51B6','PGM1','C5orf46','PLN','SLC15A1','TMEM222','CITED2','HNRPDL','PKP3','HSPA1A','UGT3A2','OR7A10','CYP3A4','OR2A14','AQP3','STK39','IPO5','SDHB','PHKB','HEPHL1','DEGS2','OR52W1','SLC2A7','KCNAB3','ASCC2','PTPRG','TMEM200A','LRRC2','OCLN','OR6W1P','POPDC2','RBP4','OR5D18','SEC14L3','MYBPC1','RASGRP2','HOXD10','CDH6','SDR9C7','PGA3','ATP5A1','VWA2','CSRP1','C15orf62','HOXB6','OR52N2','TAAR6','SLC5A4','OR9A4','CSMD2','MTNR1A','C17orf89','HTR3C','PC','ACSS2','PTPRT','IRGC','ANO7','OR1D2','OR5AR1','GSTT1','CADM2','OR1L3','OR52N5','GFAP','NR4A1','PLEKHB1','ARMC1','POFUT2','XB','PRHOXNB','PRDM6','LMOD2','OACYLP','OR52B2','LCE1F','PHLDB3','SULT1C2','KRT76','TRPM5','ATL1','LDB2','CLTA','CAP1','PTPN6','AMY2B','NAT8','OR5B3','LGALS7','CA5B','LOC100093631','AKR7L','OR2A12','KRTAP9-3','OR5V1','OR6P1','VPS52','MYO6','MICU3','SLC37A4','CFTR','TIMM23','GSDMA','FAM46B','TRANK1','TMOD4','SNTA1','FAM83G','CCDC70','ALDH1B1','KRT34','TAF13','TTC32','CDH4','CEL','OR2Z1','CLCN3','FABP4','KRT5','EDAR','KCNS3','OR2A2','ASTN1','NUDT13','FXYD6-FXYD2','EXOG','SLITRK3','GABRR2','MSRA','DDX47','OR52J3','OR6C76','TGM7','ACOT7','FOXF2','OR1F1','LRRC16B','NPY2R','GCC1','ADAD1','MPP7','GPR108','KRTDAP','OR5K1','C11orf53','OR2K2','UPK1B','TAAR5','SURF4','EXOC3','DOK3','DLD','TRPM8','FLRT1','EPC2','AK302694','TAS1R1','LY6G6C','SCGB1D1','OPCML','RAB3C','GPD1','CNN1','CLDN19','PGK1','POU6F1','OR1G1','PRKCDBP','KRT86','OTUB2','KRT24','PMCH','THBS4','PLA2G2E','PARP3','CYP2C19','ITGB4','KLK8','KRT20','CELP','AIF1L','GRIP1','RGS21','MORN4','POLR1C','TMSB4X','SLC9A3','SEC14L6','TRPV6','ADORA2B','REG4','FRRS1L','TMPRSS11F','FAM188A','LEPREL1','RNF39','TAS2R41','DZIP1L','TAC1','GCOM1','FAM117B','KCNAB2','PTPRE','HYI','GNB5','KRTAP1-1','PROSER1','ACSL5','CHPF2','MYLIP','LMNA','NAT10','SPTBN1','YTHDC2','GPR4','GLRX5','TMEM74B','KRT12','TSSK2','RAP2A','USP2','CDH24','KCNA10','RTP1','ASTN2','TAL1','OR5B21','DES','IGSF11','MYO7A','INTS4','HSF2BP','CACNG1','RGS5','AANAT','OR56B4','HRG','RELL2','SYPL2','CNNM1','GPX7','GPR155','ANKRD35','RAB3B','CALHM2','SERPINB1','OR4X2','LRRC10B','OR5B2','TMEM189-UBE2V1','FAM24B-CUZD1','POLDIP2','FCN3','ASB17','SYNJ2BP-COX16','EML2','YIF1A','DBNL','ITM2B','LRRC15','PLP1','HSD17B11','LDLRAD2','UBALD2','KRTAP26-1','CDKN3','IGSF21','DNAJC18','NT5C1A','TFF2','LECT2','CHRNA3','A3GALT2','HIST1H2AM','KLC3','MAP3K8','RASGEF1A','CHI3L1','DNAJB12','AIFM2','PRKAA2','EFR3A','OTOP2','SLC25A30','LRFN5','LCE3A','PRRT1','LGI3','ZNF259','TECTA','CYP2B7P1','BSX','OR7A17','LYG2','KIAA1143','TAF7','SPTBN4','TFAP2E','OR10R2','OR51S1','SERPINA12','TMEM215','PRKCA','LRRC8C','TFEB','AARS','MDGA1','IRX3','DKFZP434H168','PHC1','RBMXL2','OR14J1','XKR7','MANEAL','HARS','TMEM150C','SEZ6','KRT4','ARNTL','LNX2','FMO6P','PPIE','GTF2IP1','ETHE1','ABCC12','RPL3L','GSTM3','SLMO2','OR5D13','FCGBP','PON1','CNFN','GBX1','PAX8','OR4D2','LCE3E','OPN1MW','TECTB','COL3A1','BAI3','ITGBL1','SRL','PITPNB','ENTPD2','REEP6','KDM4E','ATP13A1','DSG4','IFNW1','KRTAP4-6','OR51D1','COA3','ANKRD13D','TMCC3','STX3','PODNL1','CA10','SAMD8','LRRC39','ADRA2B','SPOCK3','THEM5','KRT81','SPEG','HIST3H3','C1QL2','HIST2H2AA4','AX746670','PCMTD1','SPOCK2','TEC','RIMKLB','WIPF1','HRH3','OPRK1','TARS','PANX2','OR6K3','RNF175','TENM1','OR51E1','CREM','FZD3','KCND3','UGT2B15','OR6M1','OR8D1','FUT4','ADRB3','PABPC4','SEC24C','SLC24A4','TMC1','CILP','PADI1','ABHD16B','OR2A5','DCTN6','RHOH','PDK4','LAMB1','OR52A5','PFN2','PCDH15','NEUROD6','AIP','ELN','CNGB1','ENHO','TNNT3','TAS2R50','GBAP1','DDR1','HRASLS5','LOR','IPO9','FITM1','LCE1B','DPCR1','FANCD2OS','CAMK2A','PIR','CHIT1','NKPD1','UGT2A1','OR13F1','FAM189B','FBP2','SDAD1','PSD','MEX3C','WIPF3','OR52R1','GGN','UCKL1','NR1I2','NMBR','KIF5A','RBP2','SMYD5','ALDH1A2','KCNQ3','TENM2','EPS8L3','PPFIBP2','OR8G2','ZNF155','CYP2B6','OR2AG1','ULK1','ATP13A2','RHOBTB2','DENND6A','GRIN2A','COQ9','SPON2','BC043578','EPX','CASKIN1','USH1C','GJC1','AP2B1','GRM1','POLR3D','C2orf40','PHOX2A','IQGAP1','CLDN17','CAPN1','SFN','IRX6','OR51F1','FAAH','NCOA2','TENM4','TMEM117','LRFN4','TSPEAR','OR51A2','KRTAP7-1','RGS8','PRKAG3','LARGE','KRT2','PFKFB1','GAS2L1','MROH6','OR4K15','ALG3','PARP6','APOH','ABCC4','CACNG5','SCAND1','OVOL3','ATF6B','TGM1','HRAS','DNAJB5','COL1A2','WAC','KCNMA1','GPBAR1','MMP10','FLII','RPL32','OTUD5','LOC149086','GPR162','LCE1E','SPINK5','SORL1','OR51E2','DOCK7','RAP1GDS1','TCRBV3S1','KRT14','TRPV4','NPBWR2','TPM3','MYL1','DDX25','TMEM86A','TAS1R2','OR51T1','SNX6','GPR157','TRIM7','LRTOMT','KRT74','PDE1B','QKI','PPL','MYO5B','RFX3','DNAH5','OR10A5','EIF3G','LAMTOR1','CKMT2','TMEM38A','KRTAP12-1','RGS16','DMKN','EXO5','MOAP1','SLC2A14','CACNB2','KRT73','SLC25A45','KRT82','SLC17A4','PYGL','RYK','DAPK2','AKT2','S100A2','EVPL','OR8U8','OR52H1','SGMS1','DDX3X','ATP5I','NAT8B','SPRR2F','TCRGV','OR3A2','OR52K2','NR2F1','OR5K3','HSPB7','SPRR2E','VPS51','PGLYRP3','OR51F2','SNRNP40','SNCG','RALGDS','OR9G9','ACY3','RHBDF2','TANC2','STT3A','SCN4A','KLHL15','SSTR3','MYOZ1','SSPO','CASP14','OR6C75','KRTAP9-7','OR2T2','KRTAP9-4','OR52I1','NOMO1','OR4B1','PRNP','KRT39','ASS1','CHMP1B','SLC34A1','OR2T35','IGFL1','DEFB4B','RPL23A','NDUFAB1','KRTAP21-2','VWA7','MCTP2','DENND4B','DSC3','MAP3K6','UNC45B','ELOVL1','SH3RF1','PTGS1','PWP2','TNK1','PKLR','RPS6KB1','KLKB1','S100A3','FOXB2','SERINC4','KRTAP2-3','OR4C13','MYL6B','SRPR','C7orf55-LUC7L2','FHOD1','HIST1H1C','COL10A1','SOCS7','BAAT','PAPD7','OR10V2P','WRB','FUNDC2','GPRC5B','HLA-DQB2','KIAA1239','VDAC2','ZBTB20','MYOT','ARL5A','LRRC20','GJA1','IDH1','KRT25','KRT79','KRTAP4-9','HS6ST3','OR9A2','ABHD14B','OR2C1','PRDM12','HIST1H2BE','KIF21B','SLC25A4','HAPLN1','PDE11A','PLEK','GTSF1','PFKFB4','TUBB3','CCRN4L','SHB','FBL','MS4A15','KRTAP4-1','OR2A9P','ERCC2','DSG2','ROBO3','GALR2','OR11H6','TRIM5','ACSM4','NSF','RCHY1','DEPDC1B','MAEA','TRIM17','KRTAP19-8','KRTAP19-1','NXNL2','FFAR3','CTNNA2','DOCK11','GJB5','LOXHD1','VAT1L','PRPF4','WASF1','NNT','TREX2','RGS9BP','TP53I3','STAC3','OR5K4','OR52E4','CLEC6A','CELA2B','BACE1','CNTFR','RAPGEFL1','KRTAP10-5','CELSR2','KRTAP10-2','GLOD5','SLC28A3','COL8A1','ATP8A1','SCGB2A2','HLA-DPB1','KLHL17','TMIGD2','SRM','CASD1','SLC10A7','TSSK3','TBC1D22B','RAB27B','RBM20','KRTAP9-9','PRRT4','BEND4','GPRC5C','TTC39A','OR10A3','TCERG1L','BTBD6','GCKR','ADH7','KIAA1024L','KRTAP4-2','PER2','CAPNS1','CCT2','SCARA3','MCM7','PRDX3','C12orf66','NEURL','OLFM4','EPHA2','LCE1C','MRPL23','EBF4','TMEM33','CHN2','CPNE6','PHKG1','GBE1','ADAP1','TRIM16L','GUCA1C','STUB1','FAM206A','CYP2D7P1','HLA-DPB2','BIN1','OAF','CD1D','PCDH1','LIPK','SRSF7','GYS2','DMRTC2','VPS4B','NXPH4','LUC7L','MTAP','GCLM','NCCRP1','CD46','OR4F15','CIDEC','SPARC','GSR','SLC1A1','KLHDC7A','AMMECR1L','LYPD6','DUOX2','XKR6','PPP1R2','CNIH3','LOXL4','RFX7','MYH9','SUFU','CYP26B1','SLC2A9','NLGN2','KRTAP4-7','TXNIP','KCNQ5','NABP2','KRTAP3-2','OR5H15','OR6C2','FOXJ3','ANKRD6','KIAA1467','WDFY2','PSMD3','AADACL4','C19orf71','DFNB59','TSC22D1','DSE','REXO2','PAX2','SLC25A11','FSTL5','POU2F3','OR6C65','AIM1L','CDR2L','LRFN2','ENDOU','C14orf159','TUBA1A','STMN3','SLC6A12','PPP6R3','ZFHX2','FAM5B','KRTAP13-1','NR1I3','SPRYD4','NPY5R','PEX5L','HYOU1','COL16A1','FCER1G','ZMYND15','DOCK10','NDUFV2','BAMBI','SGIP1','CA3','HSDL2','PSTPIP2','KRT38','CHRNA4','ASL','KLHDC4','SPTBN2','HTR1B','MFHAS1','GLRX','PGA5','TEX264','PRX','BEND6','SDK2','ADRA1B','ALOX15','SEPT6','DSG1','INO80C','GMPR','NOS1','OR1F2P','ATXN7','MEIS3','OR52I2','GNRHR2','KRTAP10-10','KRR1','NIPAL4','INPP4A','OR1M1','SLC25A41','CCT4','SLC6A7','HTR2B','RXFP3','MATN1','SLC17A8','PTPRU','VAMP8','PRPS1L1','RPTN','SCGB1C1','CIDEA','NRG2','CORO1A','FAH','GUCA1B','CRCT1','PLEKHG7','VASH2','OR13C9','INO80B','CLUAP1','COPB2','SEMA6B','FBXO27','ARL2BP','NXPH1','MPV17','HPX','OR3A3','GRIK4','CLSTN3','SEMA6D','FAT1','SYT9','KSR2','CYB561A3','MYOM2','RAB7A','INHBB','OR4M1','HIST1H2BA','TUBA1C','AGPAT2','ME3','TMEM242','SDC3','VTA1','PNPLA6','TM4SF1','ATP7B','TMEM178B','LCE6A','EIF2B4','CYR61','CTRB1','MTMR1','AGO4','THBS3','GRID1','DIXDC1','PPID','NKX6-1','DNASE1L2','GAL3ST2','SDHAP1','HSD17B7P2','QRFPR','KRTAP16-1','PLXNA1','CCL18','ANKDD1B','CPNE8','ACN9','SLC52A2','PKD1L3','ZNF225','CIB1','PHB','COL6A1','LPAR3','RNF183','IDH3B','CAPSL','HCK','GALNT13','CNOT2','AGPAT6','LPCAT4','PLEKHH3','GRIN3A','PODN','TBKBP1','TMEM160','KRTAP6-3','KRTAP10-11','NRCAM','SNTN','PFKL','EPHA10','NLRP10','LCE2B','MPZ','EPHB6','ASPA','VIPR2','TGIF2','FMN1','TMC5','TRIM45','CRTC3','DENND2D','FAM53C','VPS16','CTGF','OSBPL6','TMEM209','SLC2A1','TCRBV5S1A1T','OR4C3','TBR1','ANXA8','ABCB4','TKTL2','MEDAG','PAK6','TDRD12','AAMP','EIF2B5','ZAP70','ATP4A','SNX5','EIF4G1','TNP1','PITPNM3','FAM214B','BLMH','NMT1','LRRC73','FOXN1','STK25','SYCN','VTCN1','DDX53','TTN','GLYATL3','IAPP','PNMAL2','LMX1B','KRTAP17-1','IFIT3','TBCB','SCN7A','OR10H3','GATA2','PVALB','CAMK2B','KRT121P','TBX10','STRIP1','CELA3B','MAPK10','N4BP3','ZDHHC2','SLIT3','XKRX','CCL5','PLIN4','LCE3D','POMGNT1','RPS16','OR10H5','FAM219B','BID','SGTA','CYC1','SPIRE1','PRPH','TSGA10','AKAP10','GMPS','ADPRH','UCP3','ABCE1','ZNF664','SH3GLB2','METTL24','WDR75','NSMCE2','PDHB','EML6','NUDT12','PDCD6IP','UTP11L','ANKH','RALGPS1','TRIB3','TNKS','TRPC6','SPRR3','X06774','NKAIN4','CYP4X1','IL37','SFTPA1','CCM2L','PLG','GPR113','JAG2','OR1L1','KLK10','RBP5','KLHL10','TAS2R19','DNM3','FRK','MFF','RANBP9','DCPS','CCT8','C12orf44','PDCL3','EIF4E3','MT4','SGCZ','CCDC97','KATNB1','PGLYRP1','CHIA','SEC31B','KRTAP1-4','TNN','OR4F5','SLC28A1','ITPKA','LCE2D','CYP4B1','PTPRVP','HSPA8','TAS1R3','TRAPPC10','CD300LF','SS18L1','LHX1','MAL2','CDH22','ZNF362','DHRS7C','VSTM2A','LRRC8E','PDE5A','WDFY3','DNER','TRPM3','NOL10','SPNS2','SPRYD3','GNPDA1','C16orf70','PSMB7','SLC18B1','AMY2A','ARL5C','OR2H2','TIAF1','FPR2','MRGPRE','CDK18','SCN3A','MEP1A','LYN','SPTLC2','ISCU','FAM168A','42067','PLD1','TFDP2','NPC1','PDE12','SART1','LYZL4','DGAT2L6','NYAP1','TMA7','DEPTOR','DNAJC5','PPP3CA','VSIG10L','OR8B8','ZSWIM4','FGF22','MUC19','ULBP2','CLEC12B','GLRB','OR4M2','KRT75','CYP2F1','LRRC3B','EIF2S1','RAP1GAP2','INTS9','SERPINF1','SMARCD2','CAT','PKP4','NEDD4L','CTDSP2','IMMP1L','DOK4','F13A1','CD3E','GTF2H4','MYO15B','PLK3','PKN1','BLVRB','LOC284023','BDH2','BRD7','HSPD1','MEIG1','VPS37C','VAT1','AUTS2','PLAC9','CYP2C9','RPL15','LRRC29','ITLN2','OR8J1','OR6Q1','HSP90AB1','LILRB3','TAAR8','ITLN1','SIGLEC14','COL6A4P1','CYP4A11','HDHD2','DGKG','CTNNA1','CHRM5','OPN5','TMEM63B','DTNBP1','CNOT4','RELN','UNC45A','GRK6','SNX2','PCDHGC5','ASB16','MAG','DRG2','PROM2','TMEM181','SLC7A13','TRPM4','LAMB4','AQP12B','PIEZO1','PRSS57','RAPGEF2','ITIH1','TEX261','PHEX','USP21','GNAO1','RTN4RL1','ABCA12','AGR2','BPIFC','RASGRP1','CLDN9','ZG16','ZBTB7C','DOCK8','EIF5','SNAP23','SNAPC5','ARSF','AMPD2','MRP63','ADAL','SLC27A2','SIGLEC6','CYP3A7-CYP3AP1','CRHR2','LYSMD2','NUP93','FMOD','ANK1','DYNC1LI1','PRELP','LIN7C','KCNH8','C17orf105','TMOD1','MTMR7','VCPIP1','NCDN','POF1B','GAS7','TIMP2','EIF5A2','SEMA5A','BMP8A','CDH26','TMEM236','ADCY5','ONECUT2','CMTM5','OR2A1','COPS8','SV2B','USF2','ADAM11','RAB3GAP2','OAZ3','SLC41A1','DCN','FLRT3','TRHDE','CCS','DDX3Y','TMEM244','AWAT1','ITGA6','CTXN2','FABP5','DCDC2B','RBM42','EXOSC3','LTBP4','PECR','OLFM2','TRIML1','HBB','PGA4','CMTM3','C17orf64','KRT72','OR8B12','CLDN25','NDUFAF5','OR8B2','ZNF707','C8orf22','YIPF5','CLCN6','GAD2','POU5F2','MLLT3','RAB37','KARS','TMEM82','CACNG2','RERG','FMNL2','LMAN2L','RPL10L','C2CD5','MRPL18','PAPSS1','LHX4','KRTAP19-7','AVPR1A','FAM167A','SURF1','LPO','RASSF5','SLC44A4','METTL7A','OR52B6','ST6GALNAC4','CPLX1','LYZ','FLG2','SORT1','TLK1','LEP','SLC6A4','NETO2','CNIH','DNASE1L3','TAGLN','DHRS2','OR4X1','DSG3','CRYAA','ATF5','PCID2','OR5T2','CCL3','S100A13','RNASE3','KRTAP13-2','TNNT2','SEPN1','C1QBP','GPX3','SYNPR','ARHGAP44','ROBO2','ZAK','TRIM3','MFN2','STXBP5L','SLC6A5','EXTL3','ATG4C','CLPS','DUSP14','CASP16','LCE2C','OR5D14','CEBPG','EIF2S2','PTCHD2','PCDH7','PLCD1','SMG1','TMEM185B','C1QL3','MX1','DDAH2','ASPRV1','MURC','GJB4','SRMS','TMEM11','PRM3','FAM163B','TLX3','PXMP4','S100A14','PPP1R3D','MOG','AK098763','KCNE2','UBE2O','TGFBI','SNX25','VGLL3','CIT','CHST10','LSAMP','PSMB6','FECH','WNT10A','IRS2','RPL22L1','C15orf48','EIF4EBP3','ZC2HC1B','ATP6V1G2','GCH1','CRMP1','TRAPPC11','ZNF665','FAM78B','B9D1','RPL22','NEUROD1','HIPK4','UBAP2L','GRID2','CA12','ADC','VPS33A','VCL','MYL3','TUBA4A','HNRNPA1','CLPB','USF1','CYP2E1','ADCY2','TCRBV20S1','TRAPPC3L','FCGR1A','RIN3','TAS2R3','CIRBP','DNAJB8','SIPA1','TMIE','OS9','GARS','SAG','ALDH6A1','SPINK6','CCDC64B','CSNK1A1','ZNF423','PHLPP2','IDE','OR10G9','PRR7','MYH6','HOPX','CCNA1','PRMT5','PBX3','CAPN2','ASB13','TRIM2','TSHZ3','STX16','GRM2','PEX2','KCTD5','MXI1','DSTN','SMPDL3B','COL11A1','CLRN3','KCTD21','EXOSC4','RASL12','CPZ','HS3ST3A1','NCAM1','RBM7']35:[] 2d:{"id":"7aaf4617ccc54b7d","slug":"visualizing-busted-ph-analysis-results","title":"Visualizing BUSTED-PH analysis results","thumbnail":"6d83c410e4b3cc0e8735f5de904c1fe70caef39c7651a96fc97d3f84ad0f23af","access_level":"public","create_time":"$D2022-02-11T17:10:41.570Z","update_time":"$D2022-09-08T14:31:44.315Z","view_time":null,"first_public_time":"$D2022-07-12T23:29:42.627Z","has_unpublished_edits":null,"fork":{"id":"61ca4428b1c669bd","slug":"visualizing-busted-ph-analysis-results","title":"Visualizing BUSTED-PH analysis results","owner":{"id":"9bf992815a7b66f2","login":"spond","name":"Sergei Pond"}},"stars":0,"starred":false,"owner":{"type":"team","id":"c1465becab0bcdf5","login":"averyselberg","name":"Avery Selberg","avatar_url":"https://avatars-next.observableusercontent.com/7f856996a861d933b1c34b674aeea0e0ee9950febb5f9546c2f372fea88cc2a3","tier":"starter_2024"},"version":3491,"publish_version":3491,"role":null,"authored":false,"authors":[{"type":"individual","id":"bb48a9987020e1cf","login":"averyselberg","name":"Avery Selberg","avatar_url":"https://avatars-next.observableusercontent.com/7f856996a861d933b1c34b674aeea0e0ee9950febb5f9546c2f372fea88cc2a3","tier":null}],"members":[{"type":"individual","id":"bb48a9987020e1cf","login":"averyselberg","name":"Avery Selberg","avatar_url":"https://avatars-next.observableusercontent.com/7f856996a861d933b1c34b674aeea0e0ee9950febb5f9546c2f372fea88cc2a3","tier":null,"role":"editor"}],"comments":[],"custom_thumbnail":false,"trashed":false,"lapsed":null,"subscribed":null,"forks":1,"collection":null,"type":"notebook","body":{"title":"Visualizing BUSTED-PH analysis results","cells":[{"id":0,"value":"md`# Visualizing BUSTED-PH analysis results`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3176,"value":"viewof phenotype = select ({\n \"options\" : [\"marine\",\"masstop80\"],\n \"title\" : \"Phenotype\"\n});","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2156,"value":"md`### Summary of BUSTED-PH analyses\nBUSTED-PH results are available for **${_.size (_.filter (result_json, (d)=>\"fits\" in d))}** alignments. Of those, **${busted_significant.length}** yielded at least one sigfinicant result (p≤0.05) in BUSTED-PH test. Of those \n\n* **${_.filter (table_data, (d)=>d[\"FG p-value\"]<=0.05).length}** had evidence of selection on **Foreground**\n* **${_.filter (table_data, (d)=>d[\"BG p-value\"]<=0.05).length}** had evidence of selection on **Background**\n* **${_.filter (table_data, (d)=>d[\"DIFF p-value\"]<=0.05).length}** had evidence of differences between **Foreground** and **Background**\n`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2809,"value":"table1_title = md`**Table 1** BUSTED-PH results on individual alignments`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2817,"value":"table1 = Inputs.table(averys_table)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2838,"value":"$2f","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2912,"value":"fig1_label = md`**Figure 1**. Disributions of mean ω by group`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2920,"value":"fig1=vegalite ({\n \"$schema\": \"https://vega.github.io/schema/vega-lite/v4.json\",\n \"data\": {\n \"values\": _.filter (mean_dnds_plot, (r)=>r.omega < 2)\n },\n \"width\": 800,\n \"height\": 50,\n \"transform\":[{\n \"density\": \"omega\",\n \"groupby\": [\"clade\"],\n \"bandwidth\": 0.025\n }],\n \"mark\": {type: \"area\", \"opacity\" : 0.5, tooltip : true, line : true},\n \"encoding\": {\n \"x\": {\n \"field\": \"value\",\n \"grid\" : null,\n \"title\": \"Density\",\n \"type\": \"quantitative\",\n \"axis\": {\"grid\": false}\n },\n \"y\": {\n \"field\": \"density\",\n \"type\": \"quantitative\",\n \"title\" : \"\",\n \"axis\": {\"grid\": false}\n },\n \"color\" : {\"value\" : \"grey\"},\n \"row\" : {\"field\" : \"clade\", \"title\" : \"dN/dS\", \"sort\": \"ascending\"}\n }\n})\n\n","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2943,"value":"table2_label = md`**Table 2**. Disributions of mean ω by group`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2946,"value":"table2 = table (mean_dnds_groups,{\n limit: 5000,\n enableCSVDownload: true,\n sort: true,\n enableFilter: false,\n columns: [\n {\n key : 'group',\n name: \"Group\"\n }, \n {\n key : 'mean',\n name: \"Mean\",\n render: (val)=>float_format (val)\n }, \n {\n key : 'median',\n name: \"Median\",\n render: (val)=>float_format (val)\n }, \n {\n key : 'Q1',\n name: \"5%\",\n render: (val)=>float_format (val)\n },\n {\n key : 'Q3',\n name: \"95%\",\n render: (val)=>float_format (val)\n }\n ]\n})\n ","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2992,"value":"md`### Summary of BUSTED[S] analyses\nBUSTED[S] results are available for **${_.size (_.filter (result_json, (d)=>\"BUSTED\" in d))}** alignments. Of those, **${busted_significant.length}** yeilded a sigfinicant result (p≤0.05) for epidosic selection in at least one group.\n`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2233,"value":"md`### Gene-by-gene information`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":966,"value":"viewof gene_id = {\n const dd3 = autoSelect ({\n title: \"Gene ID\",\n options: _.keys (result_json)\n \n });\n dd3.input.style.fontSize = \"12px\";\n dd3.input.style.marginTop = \"8px\";\n return dd3;\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2643,"value":"md`#### Standard RELAX test`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2636,"value":"md`\n**${gene_id}** has \n${_.get (result_json,[gene_id,\"sequences\"])} sequences and ${_.get (result_json,[gene_id,\"sites\"])} sites.\n\n_Test summary._\n\n1. The p-value for selection of **foreground** branches **${float_format(_.get (result_json,[gene_id,\"tests\",\"FG\",\"p-value\"]))}**.\n2. The p-value for selection of **background** branches **${float_format(_.get (result_json,[gene_id,\"tests\",\"BG\",\"p-value\"]))}**.\n3. The p-value for differences in selection between the two sets of branches **${float_format(_.get (result_json,[gene_id,\"tests\",\"DIFF\",\"p-value\"]))}**.\n\n\n_Interpretation:_\n\n **${gene_id}** shows ${decision_tree}.\n`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2137,"value":"fig2_label = md`**Figure 2**. Rate distributions for this gene using the BUSTED-PH model`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2134,"value":"fig2 = plotPie (_.flatten (_.map (_.get (result_json, [gene_id,\"distributions\",\"unconstrained\",\"Rate Distributions\"]), (d,k)=> {\n return _.map (d, (r)=>{\n let rp = r;\n rp.type = k;\n return rp;});\n \n})) || [], gene_id)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2662,"value":"fig3_label = md`**Figure 3**. Rate distributions for this gene using the BUSTED-PH model, where background and foreground share the same distribution`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2665,"value":"fig3 = plotPie (_.flatten (_.map (_.get (result_json, [gene_id,\"distributions\",\"shared\",\"Rate Distributions\"]), (d,k)=> {\n return _.map (d, (r)=>{\n let rp = r;\n rp.type = k;\n return rp;});\n \n})) || [], gene_id)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2236,"value":"md`### Tree for ${gene_id}`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1941,"value":"trendLegend = swatches({\n color: species_colors,\n title: \"Viral species\",\n //tickFormat: \"%\"\n})","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":870,"value":"html`\n \n

`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":5,"value":"md`---`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2701,"value":"md`### Libraries and data`","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":9,"value":"d3 = require (\"d3\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":13,"value":"_=require(\"lodash\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1905,"value":"tree_object = [new phylotree.phylotree(_.get (result_json, [gene_id,\"tree\"]) || \"()\"), _.get (result_json, [gene_id,\"tested\",0]) || [], _.get (result_json, [gene_id,\"lengths\"])]","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":127,"value":"$30","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":874,"value":"import { table } from \"@gampleman/table\"","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":174,"value":"import {select,checkbox, autoSelect} from \"@jashkenas/inputs\"","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":488,"value":"vegalite = require(\"vega-embed\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2288,"value":"float_format = d3.format (\".3g\");","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2217,"value":"$31","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"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","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1929,"value":"species = new Set(_.flatten(_.map (has_relax, (f)=>_.map (_.uniqBy (_.toPairs (f[\"tested\"][0]), (d)=>d[1]), (d)=>d[1]))))","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3030,"value":"has_relax = _.chain (result_json).toPairs().filter ((d)=>\"tested\" in d[1]).fromPairs().value()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":1937,"value":"species_colors = d3.scaleOrdinal(d3.schemeCategory10).domain ([...species]);","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2785,"value":"table_data = _.map (_.toPairs (result_json), (d)=> {\n let res = {\n \"ID\" : d[0] \n };\n \n let pfg = (d)=>d.tests.FG[\"p-value\"];\n let pbg = (d)=>d.tests.BG[\"p-value\"];\n let pdiff = (d)=>d.tests.DIFF[\"p-value\"];\n let haz = (d)=>pfg(d)<=0.05 && pbg (d)>0.05 && pdiff (d)<=0.05;\n let sites = (d)=>[d.fg_selected, d.bg_selected];\n let l = (d)=>_.mapValues(_.groupBy (_.toPairs(d.lengths), (i)=>d.tested[0][i[0]]), (ll)=>d3.sum (ll, (i)=>i[1]));\n \n let inject = [\n [\"Sequences\", [\"sequences\"]],\n [\"Sites\", [\"sites\"]],\n [\"FG p-value\", pfg],\n [\"BG p-value\", pbg],\n [\"DIFF p-value\", pdiff],\n [\"SelSites\", sites],\n [\"PH\", haz],\n [\"L\", l]\n \n ];\n \n _.each (inject, (i)=> {\n res[i[0]] = _.isFunction (i[1]) ? i[1](d[1], d[0]) :_.get (d[1], i[1]);\n });\n \n return res;\n});","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3317,"value":"$32","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3276,"value":"$33","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2998,"value":"busted_significant = _.filter (_.toPairs (result_json), (f)=>f[1].tests && d3.min (_.map(f[1].tests, (test)=>test[\"p-value\"])) <= 0.05)","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2475,"value":"baseline_by_segment = _.chain (result_json).toPairs().filter ((d)=>\"fits\" in d[1]).map ((d)=> {\n d[1]=\"Reference\";\n return d;\n}).fromPairs().value();","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2952,"value":"mean_dnds_groups = _.chain (_.groupBy (mean_dnds_plot, (d)=>d.clade)).map ((d,k)=>{\n return {\n group : k,\n mean : d3.mean (d, (e)=>e.omega),\n median : d3.median (d, (e)=>e.omega),\n Q1 : d3.quantile (d, 0.05, (e)=>e.omega),\n Q3 : d3.quantile (d, 0.95, (e)=>e.omega)\n };\n}).value()","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2904,"value":"mean_dnds_plot = _.filter (_.flatten (_.map (mean_dnds, (d, f)=> {\n return _.map (d, (dnds, gr)=> { \n return {'clade' : gr, \"omega\" : dnds, \"gene\" : f};\n });\n})), (r)=>_.isNumber (r.omega) && r.clade != \"Unclassified\")","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":2865,"value":"mean_dnds = {\n return _.mapValues (result_json, (d,f)=> {\n if (\"fits\" in d == false) return; \n let test_groups = new Set (_.values (d.tested[0]));\n let mapper = {\"Unclassified\" : \"Unclassified\"};\n if (test_groups.has (\"Test\")) {\n _.some (d.tested[0], (g,s) => {\n if (g == \"Test\") {\n mapper[\"Test\"] = s.split (\"_\")[0]; \n } else {\n mapper[\"Reference\"] = s.split (\"_\")[0]; \n }\n return _.size (mapper) == 3;\n }\n );\n } else {\n mapper = _.fromPairs ( _.map ([...test_groups], (d)=>[d,d]));\n }\n return _.fromPairs ([['mapping', mapper]].concat(_.map ([...test_groups], (g)=> {\n return ([mapper[g],d.fits[\"MG94xREV with separate rates for branch sets\"][\"Rate Distributions\"][\"non-synonymous/synonymous rate ratio for *\" + g + \"*\"][0][0]]);\n })));\n })\n}","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":754,"value":"result_json = await d3.json (\"https://data.hyphy.org/web/busted-ph/\"+phenotype+\"_ags.json\");","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3466,"value":" codeml_missed = {\n let missed = []\n for (const gene in result_json){\n if(result_json[gene][\"tests\"][\"FG\"][\"p-value\"] < 0.05 && result_json[gene][\"tests\"][\"BG\"][\"p-value\"] < 0.05 && result_json[gene][\"tests\"][\"DIFF\"][\"p-value\"] > 0.05){\n missed.push(gene)\n }\n }\n return missed\n} \n/*\ncodeml_missed = {\n let missed = []\n for (const gene in result_json){\n if(result_json[gene][\"tests\"][\"FG\"][\"p-value\"] < 0.05 && result_json[gene][\"tests\"][\"BG\"][\"p-value\"] > 0.05 && result_json[gene][\"tests\"][\"DIFF\"][\"p-value\"] > 0.05){\n missed.push(gene)\n }\n }\n return missed\n} */","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3477,"value":"$34","pinned":false,"hidden":false,"mode":"ojs","author":"$undefined","summary":"$undefined","output":"$undefined","database":"$undefined"},{"id":3473,"value":"missed_marine = {\n const missed = []\n for (let i = 0; i