library(plyr)
library(tidyverse)
library(gridExtra)
library(cowplot)
library(data.table)
library(extrafont)
loadfonts(device = "win")
par(family = "LM Roman 10")

#source("CI-Functions.R")
source("CI-Functions-Bonferroni.R")

defaultpath <- dirname(rstudioapi::getActiveDocumentContext()$path)
setwd(defaultpath)


# ALL TASKS
path <- "plots/all_tasks_mapSize/raw/"
plots <- list()
vis_all <- c("GlyMa","CoordV", "LocSto")
vis_names <- c("Glyph Maps","Coordinated Views", "Location Storylines")
metrics <- c('time', 'error')
metrics_colors <- c('dodgerblue2', 'darkorange3')
task <- 'all'
task_name <- 'all'


## SETUP FILES
time_per_task <- file("./plots/all_tasks_mapSize/time_all_tasks_latex_table.txt")
error_per_task <- file("./plots/all_tasks_mapSize/error_all_tasks_latex_table.txt")

all_lines_time <- c("\\begin{table}[]","\\begin{tabular}{l|l|l}", " & Mean times & Mean Time diffs \\\\ \\hline")
all_lines_error <- c("\\begin{table}[]","\\begin{tabular}{l|l|l}", " & Mean Errors & Mean Error diffs \\\\ \\hline")
next_line <- 4



for (v in 1:3) {
  vis = vis_all[v]
  vis_name = vis_names[v]
  
  for (m in 1:2) {
    metric <- metrics[m]
    # READ METRIC
    filename_metric <- paste0(path,'printed_phase_', metric, '_means_task_', task, '_', vis,'.txt', sep='')
    printed_metric <- read.table(filename_metric, header = T, sep = ',')
    colnames(printed_metric)[colnames(printed_metric)=="upperBound_CI."] <- "upperBound_CI"
    #printed_metric$vis <- factor(printed_metric$phase, levels=rev(levels(printed_metric$phase)))
    
    filename_metricdiff <- paste0(path,'printed_phase_', metric,'_diffs_task_', task, '_', vis, '.txt', sep='')
    printed_metricdiff <- read.table(filename_metricdiff, header = T, sep = ',')
    #printed_metricdiff$phase <- factor(printed_metricdiff$phase, levels=c("Adv-Explo"))
    
    
    new_line <- paste0(vis_name, "&", "\\begin{tabular}[c]{@{}l@{}}")
    
    for (c in 1:2) {
      tech <- printed_metric[[1]][[c]]
      mean <- format(round(as.numeric(printed_metric[[2]][[c]]), 2), nsmall = 2)
      minCI_raw <- round(as.numeric(printed_metric[[3]][[c]]), 2)
      maxCI_raw <- round(as.numeric(printed_metric[[4]][[c]]), 2)
      minCI <- format(min(minCI_raw, maxCI_raw), nsmall = 2)
      maxCI <- format(max(minCI_raw, maxCI_raw), nsmall = 2)
      new_line <- paste0(new_line, tech, " ", mean, "[", minCI, ",", maxCI, "]")
      if(c != 3){
        new_line <- paste(new_line, "\\\\")
      }
    }
    new_line <- paste(new_line, "\\end{tabular}", "&")
    new_line <- paste(new_line, "\\begin{tabular}[c]{@{}l@{}}")
    
    for (c in 1:1) {
      tech <- printed_metricdiff[[1]][[c]]
      mean <- format(round(as.numeric(printed_metricdiff[[2]][[c]]), 2), nsmall = 2)
      minCI_raw <- round(as.numeric(printed_metricdiff[[3]][[c]]), 2)
      maxCI_raw <- round(as.numeric(printed_metricdiff[[4]][[c]]), 2)
      minCI <- format(min(minCI_raw, maxCI_raw), nsmall = 2)
      maxCI <- format(max(minCI_raw, maxCI_raw), nsmall = 2)
      new_line <- paste0(new_line, tech, " ", mean, "[", minCI, ", ", maxCI, "]")
      if(c != 3){
        new_line <- paste(new_line, "\\\\")
      }
    }
    new_line <- paste(new_line,"\\end{tabular}", "\\\\")
    
    if(m == 1){
      all_lines_time[[next_line]] <- new_line
    } else {
      all_lines_error[[next_line]] <- new_line
      next_line <- next_line + 1
    }
    
  }
}

all_lines_time[[next_line]] <- "\\end{tabular}"
all_lines_error[[next_line]] <- "\\end{tabular}"
next_line <- next_line + 1
all_lines_time[[next_line]] <- "\\end{table}"
all_lines_error[[next_line]] <- "\\end{table}"


writeLines(c(all_lines_time), time_per_task)
writeLines(c(all_lines_error), error_per_task)

close(time_per_task)
close(error_per_task)