plot script fixes

This commit is contained in:
Alwin Berger 2025-05-13 12:09:58 +00:00
parent a0ab58a6c9
commit 12396f7104

View File

@ -1,3 +1,4 @@
# install.packages(c("mosaic", "dplyr", "DBI", "tikzDevice", "colorspace", "heatmaply", "RColorBrewer", "RSQLite"))
library("mosaic") library("mosaic")
library("dplyr") library("dplyr")
library("DBI") library("DBI")
@ -103,10 +104,13 @@ BREW=RdYlGn(4)
draw_plot <- function(data, casename) { draw_plot <- function(data, casename) {
# evo, cov, random, fret # evo, cov, random, fret
MY_COLORS <- c(BREW[[3]], BREW[[2]], BREW[[1]], BREW[[4]], "cyan", "pink", "gray", "orange", "black", "yellow","brown") MY_COLORS <- c(BREW[[4]], BREW[[3]], BREW[[2]], BREW[[1]], "cyan", "pink", "gray", "orange", "black", "yellow","brown")
# MY_COLORS <- c("orange", "blue", "red", "green", "orange", "cyan", "pink", "gray", "orange", "black", "yellow","brown") # MY_COLORS <- c("orange", "blue", "red", "green", "orange", "cyan", "pink", "gray", "orange", "black", "yellow","brown")
# MY_COLORS <- c("green", "blue", "red", "magenta", "orange", "cyan", "pink", "gray", "orange", "black", "yellow","brown") # MY_COLORS <- c("green", "blue", "red", "magenta", "orange", "cyan", "pink", "gray", "orange", "black", "yellow","brown")
LEGEND_POS=LEG_POS[[casename]] LEGEND_POS=LEG_POS[[casename]]
if (is.null(LEGEND_POS)) {
LEGEND_POS = "bottomright"
}
ISNS_PER_US = (10**3)/(2**5) ISNS_PER_US = (10**3)/(2**5)
@ -120,6 +124,9 @@ draw_plot <- function(data, casename) {
data[[n]]$sdiv <- data[[n]]$sdiv / ISNS_PER_US data[[n]]$sdiv <- data[[n]]$sdiv / ISNS_PER_US
} }
data <- data[c('stgwoet', 'feedgeneration100', 'frafl', 'random')] # manual re-order
data <- data[!sapply(data, is.null)] # remove NULL entries
wcrt = KNOWN_WCRT[[casename]] wcrt = KNOWN_WCRT[[casename]]
if (!is.null(wcrt)) { if (!is.null(wcrt)) {
wcrt = wcrt / ISNS_PER_US wcrt = wcrt / ISNS_PER_US
@ -127,7 +134,7 @@ draw_plot <- function(data, casename) {
wcrt = 0 wcrt = 0
} }
static_wcrt = STATIC_WCRT[[casename]] static_wcrt = STATIC_WCRT[[casename]]
if (!is.null(wcrt)) { if (!is.null(static_wcrt)) {
static_wcrt = static_wcrt / ISNS_PER_US static_wcrt = static_wcrt / ISNS_PER_US
} else { } else {
static_wcrt = 0 static_wcrt = 0