plot_specs = (
{
"Evidence ratio for ω>1 (constrained)" : {
"width": 800, "height": 150,
"vconcat" : _.map (_.range (1, fig1data.length + 1, 70), (d)=> {
return ERPlot (fig1data, d, 70, "ER (constrained)")
})},
"Evidence ratio for ω>1 (optimized)" : {
"width": 800, "height": 150,
"vconcat" : _.map (_.range (1, fig1data.length + 1, 70), (d)=> {
return ERPlot (fig1data, d, 70, "ER (optimized null)")
})},
"Synonymous rates" : {
"width": 800, "height": 150,
"vconcat" : _.map (_.range (1, fig1data.length + 1, 70), (d)=> {
return SRVPlot (fig1data, d, 70, "SRV posterior mean", srv_hmm ? "SRV viterbi" : null)
})},
"Support for positive selection" : {
"vconcat" : _.map (_.range (1, results_json.input["number of sites"], 1000), (d)=> {
return BSPosteriorPlot (bsPositiveSelection, d, 1000)
})},
"Joy Plot" : {
"vconcat" : _.map (_.range (1, results_json.input["number of sites"], 7000), (d)=> {
return BSPosteriorPlot_Joy (bsPositiveSelection, d, 7000)
})},
"Support for 2H" : {
"vconcat" : _.map (_.range (1, results_json.input["number of sites"], 70), (d)=> {
return BSPosteriorPlot (twoHBranchSite, d, 70)
})},
"Support for 3H" : {
"vconcat" : _.map (_.range (1, results_json.input["number of sites"], 70), (d)=> {
return BSPosteriorPlot (threeHBranchSite, d, 70)
})},
"Support for 2H+3H" : {
"vconcat" : _.map (_.range (1, results_json.input["number of sites"], 70), (d)=> {
return BSPosteriorPlot (multiHBranchSite, d, 70)
})},
"Compare ER" : ERPlotComp (fig1data, scatter_plot[0], scatter_plot[1], scatter_plot[2]),
"Compare BS support" : BSPlotComp (compBSModels, scatter_plot[2]),
"LR difference CDS" : cdsQuant (fig1data, "Diff LR"),
"LogL difference" : {
"width": 800, "height": 150,
"vconcat" : _.map (_.range (1, fig1data.length + 1, 70), (d)=> {
return ERPlot (fig1data, d, 70, "Diff logl")
})},
"LogL difference scatter" : ScatterPlotComp (fig1data, "LogL ref", "LogL comp", scatter_plot[2]),
"LR difference scatter" : ScatterPlotComp (fig1data, "LR ref", "LR comp", scatter_plot[2]),
"LR difference" : {
"width": 800, "height": 150,
"vconcat" : _.map (_.range (1, fig1data.length + 1, 70), (d)=> {
return ERPlot (fig1data, d, 70, "Diff LR")
})},
}
)