R Notebook

library('readr')
library('dplyr')
library('ggplot2')
library('ggthemes')

read_data <- function(fn)
{
    read_csv(fn, col_types = 'iniin')
}
#data1834_0_4 <- read_data('decoder-out-hpar_new_1834_0-4.csv')

make_plots <- function(alldata, modelname=waiver()) 
{
    alldata <- filter(alldata, abs(coef) < 4)
    components <- unique(alldata$component)
    lapply(components,
           function(comp) {
               pltdata <- filter(alldata, component==comp) %>% mutate(coef=factor(coef), rgn=paste('rgn', rgn))
               ggplot(data=pltdata, aes(x=hour, y=value, color=coef)) + 
                   geom_line(size=0.5) +
                   facet_wrap(~rgn) +
                   ggtitle(paste('Component',comp), subtitle=modelname) +
                   theme_solarized(light=FALSE) +
                   scale_color_solarized()
           })
}

drawplots <- function(plotlist) {
    for (p in plotlist) {
        print(p)
    }
}

top5nets <- c(net1='1813_0-1', net2='1519_0-7', net3='1400_0-1', net4='1880_0-6', net5='117_0-2')
filenames <- paste0('encoding_',top5nets, '.csv')
modelnames <- paste('model', top5nets)

names(filenames) <- names(top5nets)
datasets <- lapply(filenames, read_data)

all_model_plts <- mapply(make_plots, datasets, modelnames, SIMPLIFY=FALSE)

Model 1 (1813_0-1)

drawplots(all_model_plts$net1)

Model 2 (1519_0-7)

drawplots(all_model_plts$net2)

Model 3 (1400_0-1)

drawplots(all_model_plts$net3)

Model 4 (1880_0-6)

drawplots(all_model_plts$net4)

Model 5 (117_0-2)

drawplots(all_model_plts$net5)