R/comet.R

In TiphaineCMartin/coMET: coMET: visualisation of regional epigenome-wide association scan (EWAS) results and DNA co-methylation patterns

##########################################################################
# This program is free software; you can redistribute it and/or modify
# it under the terms of the GNU General Public License as published by
# the Free Software Foundation; either version 2 of the License, or
# (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the Free Software

##########################################################################

##########################################################################
# File: comet.R
# Author: Tiphaine Martin
# Email: Tiphaine.Martin@kcl.ac.uk
# Purpose: coMET allows the display of p-values from association
#           with a correlation heatmap.
# Version : 0.99.9
###########################################################################

#coMET allows the display of p-values from association with a correlation heatmap.
#This version coMET.web is used for web site.

# The values given such as attributes of function erase ones of configuration file.
comet <- function(mydata.file = NULL,
                  mydata.format = "site",
                  mydata.type = "file",
                  mydata.large.file = NULL,
                  mydata.large.format = "site",
                  mydata.large.type = "listfile",
                  cormatrix.file = NULL,
                  cormatrix.method = "spearman",
                  cormatrix.format= "raw",
                  cormatrix.color.scheme = "bluewhitered",
                  cormatrix.conf.level=0.05,
                  cormatrix.sig.level= 1,
                  cormatrix.adjust="none",
                  cormatrix.type = "listfile",
                  mydata.ref = NULL,
                  start = NULL,
                  end = NULL,
                  zoom = FALSE,
                  lab.Y = "log",
                  pval.threshold = 10e-6,
                  pval.threshold.2 = 0,
                  disp.pval.threshold = 1 ,
                  disp.association = "FALSE",
                  disp.association.large = "FALSE",
                  disp.region = "FALSE",
                  disp.region.large = "FALSE",
                  disp.beta.association = "FALSE",
                  disp.beta.association.large = "FALSE",
                  factor.beta = 0.3,
                  symbols = "circle-fill",
                  symbols.large = NA,
                  sample.labels = NULL,
                  sample.labels.large = NULL,
                  use.colors = TRUE,
                  disp.color.ref = TRUE,
                  color.list = NULL,
                  color.list.large = NULL,
                  disp.mydata = TRUE,
                  biofeat.user.file = NULL,
                  biofeat.user.type = NULL,
                  biofeat.user.type.plot = NULL,
                  genome="hg19",
                  dataset.gene = "hsapiens_gene_ensembl",
                  tracks.gviz = NULL,
                  disp.mydata.names = TRUE,
                  disp.color.bar = TRUE,
                  disp.phys.dist = TRUE,
                  disp.legend = TRUE,
                  disp.marker.lines = TRUE,
                  disp.cormatrixmap = TRUE,
                  disp.pvalueplot = TRUE,
                  disp.type = "symbol",
                  disp.mult.lab.X = FALSE,
                  disp.connecting.lines = TRUE,
                  palette.file = NULL,
                  image.title = NULL,
                  image.name = "coMET",
                  image.type = NULL,
                  image.size = 3.5,
                  fontsize.gviz =5,
                  font.factor = 1,
                  symbol.factor = NULL,
                  print.image = TRUE,
                  connecting.lines.factor = 1.5,
                  connecting.lines.adj = 0.01,
                  connecting.lines.vert.adj = -1,
                  connecting.lines.flex = 0,
                  config.file = NULL,
                  verbose = FALSE) {

  #-------------------MAIN FUNCTION----------------------------

  #DEBUG STATEMENT
  if (verbose) cat("START COMET VERSION GENERIC\n")

  #-------------------GLOBAL VARIABLES BEGINS------------------

  #create three variables (global variables), which are lists of all other variables
  #used throughout the program to be passed to functions with values changed as
  #necessary.
  gbl.var <- NULL

  # remove trackview and ggbio in the possible visualisation
  tracks.ggbio <- NULL
  tracks.trackviewer <- NULL

  #DATA VARIABLES
  presence.mydata<- 0 # 0 presence of omic-WAS results (object mydata), 1 no presence of data
  general.data <- NULL
  split.biofeature.data.user.file <- NULL
  split.biofeature.data.user.type <- NULL
  split.biofeature.data.user.type.plot <- NULL

  mydata.data <- NULL
  mydata.num <- 0
  mydata.names <- NULL
  mydata.chr <- NULL
  mydata.gen <- NULL
  mydata.ref.pos <- NULL
  myfactor.beta <- 0

  sorted.mydata.names <- NULL
  sorted.mydata.pos <- NULL
  sorted.mydata.pos.zoom <- NULL

  cur.sample <- NULL
  mydata.samples <- NULL

  mydata.hash.names.pos <- new.env(hash=TRUE)
  mydata.hash.list.names.pos <- new.env(hash=TRUE)
  mydata.hash.names.start <- new.env(hash=TRUE)
  mydata.hash.names.end <- new.env(hash=TRUE)
  mydata.hash.pos.names <- new.env(hash=TRUE)
  mydata.hash.list.pos.names <- new.env(hash=TRUE)
  mydata.large.hash.pos.names <- new.env(hash=TRUE)
  mydata.large.hash.names.start <- new.env(hash=TRUE)
  mydata.large.hash.names.end <- new.env(hash=TRUE)
  mydata.large.hash.names.pos <- new.env(hash=TRUE)

  split.mydata.file <- NULL
  split.cormatrix.file <- NULL

  #GRAPH BOUNDARY VARIABLES

  pval.flag <- TRUE

  cur.exp <- 0

  total.dist <- NULL
  min.dist <- NULL
  max.dist <- NULL

  min.x <- NULL
  min.user.x <- NULL
  min.y <- NULL

  max.x <- NULL
  max.user.x <- NULL
  max.y <- NULL

  r.x <- NULL
  r.y <- NULL

  axis.x <- NULL
  axis.y <- NULL

  equidis.pos <- NULL

  mean.beta <- NULL
  sd.beta <- NULL

  #FORMATTING VARIABLES

  split.sample.labels <- NULL
  split.color.list <- NULL
  color.list <- NULL
  symbol.list <- NULL
  fill.list <- NULL
  split.format <- NULL
  split.association <- NULL
  split.beta.association <- NULL
  split.region <- NULL
  split.type <- NULL

  large.split.sample.labels <- NULL
  large.split.color.list <- NULL
  large.color.list <- NULL
  large.symbol.list <- NULL
  large.fill.list <- NULL
  large.split.format <- NULL
  large.split.association <- NULL
  large.split.beta.association <- NULL
  large.split.region <- NULL
  large.split.type <- NULL
  ref <- NULL
  samples <- NULL

  palette.size <- NULL
  cex.factor <- 0.5
  cex.factor.symbol <- 0.25
  font.size <- NULL
  line.width <- 0.5
  myfactor.beta <- 0.3

  cormatrix.data <- NULL
  split.cormatrix.type <- NULL
  matrix.data <- NULL
  cormatrix.data.full <- NULL
  cormatrix.pvalue.data <- NULL
  cormatrix.pvalue.data.full <- NULL
  cormatrix.CI.data <- NULL

  mySession <- NULL
  listtracks_user <- NULL
  listtracks_gviz <- NULL
  listtracks_ggbio <- NULL
  listtracks_trackviewer <- NULL

  gbl.var <- list(mydata.data = mydata.data,
                  general.data = general.data,
                  presence.mydata=presence.mydata,
                  split.mydata.file = split.mydata.file,
                  split.format = split.format,
                  split.association =  split.association,
                  split.beta.association =  split.beta.association,
                  split.region = split.region,
                  split.type = split.type,
                  split.biofeature.data.user.file = split.biofeature.data.user.file,
                  split.biofeature.data.user.type = split.biofeature.data.user.type,
                  large.split.format = large.split.format,
                  large.split.type = large.split.type,
                  large.split.association =  large.split.association,
                  large.split.beta.association =  large.split.beta.association,
                  large.split.region = large.split.region,
                  mySession = mySession,
                  listtracks_gviz = listtracks_gviz,
                  listtracks_ggbio = listtracks_ggbio,
                  listtracks_trackviewer = listtracks_trackviewer,
                  listtracks_user = listtracks_user,
                  mydata.num = mydata.num,
                  mydata.names = mydata.names,
                  mydata.chr = mydata.chr,
                  mydata.gen = mydata.gen,
                  sorted.mydata.names = sorted.mydata.names,
                  sorted.mydata.pos = sorted.mydata.pos,
                  sorted.mydata.pos.zoom = sorted.mydata.pos.zoom,
                  mydata.ref.pos = mydata.ref.pos,
                  cur.sample = cur.sample,
                  mydata.samples = mydata.samples,
                  mydata.hash.names.pos = mydata.hash.names.pos,
                  mydata.hash.list.names.pos = mydata.hash.list.names.pos,
                  mydata.hash.pos.names = mydata.hash.pos.names,
                  mydata.hash.list.pos.names = mydata.hash.list.pos.names,
                  mydata.hash.names.start = mydata.hash.names.start,
                  mydata.hash.names.end = mydata.hash.names.end,
                  mydata.large.hash.names.pos = mydata.large.hash.names.pos,
                  mydata.large.hash.names.start = mydata.large.hash.names.start,
                  mydata.large.hash.names.end = mydata.large.hash.names.end,
                  myfactor.beta = myfactor.beta,
                  mean.beta = mean.beta,
                  sd.beta = sd.beta,
                  cormatrix.data = cormatrix.data,
                  split.cormatrix.type = split.cormatrix.type,
                  split.cormatrix.file = split.cormatrix.file,
                  cormatrix.pvalue.data = cormatrix.pvalue.data,
                  cormatrix.pvalue.data.full = cormatrix.pvalue.data.full,
                  matrix.data = matrix.data,
                  cormatrix.data.full = cormatrix.data.full,
                  cormatrix.CI.data = cormatrix.CI.data,
                  samples = samples,
                  pval.flag = pval.flag,
                  cur.exp = cur.exp,
                  total.dist = total.dist,
                  min.dist = min.dist,
                  max.dist = max.dist,
                  min.x = min.x,
                  min.user.x = min.user.x,
                  min.y = min.y,
                  max.x = max.x,
                  max.user.x = max.user.x,
                  max.y = max.y,
                  r.x = r.x,
                  r.y = r.y,
                  axis.x = axis.x,
                  axis.y = axis.y,
                  equidis.pos = equidis.pos,
                  split.sample.labels = split.sample.labels,
                  split.color.list = split.color.list,
                  color.list = color.list,
                  symbol.list = symbol.list,
                  fill.list = fill.list,
                  large.split.sample.labels = large.split.sample.labels,
                  large.split.color.list = large.split.color.list,
                  large.color.list = large.color.list,
                  large.symbol.list = large.symbol.list,
                  large.fill.list = large.fill.list,
                  palette.size = palette.size,
                  cex.factor = cex.factor,
                  cex.factor.symbol = cex.factor.symbol,
                  font.size = font.size,
                  line.width = line.width,
                  verbose = verbose
  )

  #-------------------GLOBAL VARIABLES ENDS------------------

  #-------------------CONFIGURATION VARIABLES BEGINS---------
  list.tracks <- NULL
  DATASET.SNP="hsapiens_snp"
  VERSION.DBSNP="snp138"
  DATASET.SNP.STOMA="hsapiens_snp_som"
  DATASET.REGULATION="hsapiens_feature_set"
  DATASET.STRU="hsapiens_structvar"
  DATASET.STRU.STOMA="hsapiens_structvar_som"
  pattern.regulation="GM12878"
  BROWSER.SESSION="UCSC"

  #-------------------UPDATE CONFIGURATION VARIABLES---------
  config.var <- list(mydata.file = mydata.file,
                     mydata.format = mydata.format,
                     mydata.type = mydata.type,
                     mydata.large.file = mydata.large.file,
                     mydata.large.format = mydata.large.format,
                     mydata.large.type = mydata.large.type,
                     cormatrix.file = cormatrix.file,
                     cormatrix.method = cormatrix.method,
                     cormatrix.format = cormatrix.format,
                     cormatrix.color.scheme = cormatrix.color.scheme,
                     cormatrix.conf.level = cormatrix.conf.level,
                     cormatrix.sig.level = cormatrix.sig.level,
                     cormatrix.adjust = cormatrix.adjust,
                     cormatrix.type = cormatrix.type,
                     mydata.ref = mydata.ref,
                     color.list = color.list,
                     use.colors = use.colors,
                     disp.color.ref = disp.color.ref,
                     biofeat.user.file = biofeat.user.file,
                     biofeat.user.type = biofeat.user.type,
                     biofeat.user.type.plot = biofeat.user.type.plot,
                     start = start,
                     end = end,
                     zoom = zoom,
                     lab.Y = lab.Y,
                     fontsize.gviz = fontsize.gviz,
                     disp.pval.threshold = disp.pval.threshold,
                     pval.threshold = pval.threshold,
                     pval.threshold.2 = pval.threshold.2,
                     disp.type = disp.type,
                     disp.cormatrixmap = disp.cormatrixmap,
                     disp.pvalueplot = disp.pvalueplot,
                     disp.mydata = disp.mydata,
                     disp.association = disp.association,
                     disp.association.large = disp.association.large,
                     disp.beta.association = disp.beta.association,
                     disp.beta.association.large = disp.beta.association.large,
                     factor.beta = factor.beta,
                     disp.region = disp.region,
                     disp.region.large = disp.region.large,
                     sample.labels = sample.labels,
                     sample.labels.large = sample.labels.large,
                     symbols = symbols,
                     symbols.large = symbols.large,
                     genome = genome,
                     dataset.gene = dataset.gene,
                     DATASET.SNP = DATASET.SNP,
                     DATASET.SNP.STOMA = DATASET.SNP.STOMA,
                     DATASET.REGULATION = DATASET.REGULATION,
                     DATASET.STRU = DATASET.STRU,
                     DATASET.STRU.STOMA = DATASET.STRU.STOMA,
                     VERSION.DBSNP = VERSION.DBSNP,
                     pattern.regulation = pattern.regulation,
                     BROWSER.SESSION = BROWSER.SESSION,
                     tracks.gviz = tracks.gviz,
                     tracks.ggbio = tracks.ggbio,
                     tracks.trackviewer = tracks.trackviewer,
                     list.tracks = list.tracks,
                     disp.mydata.names = disp.mydata.names,
                     disp.marker.lines = disp.marker.lines,
                     disp.connecting.lines = disp.connecting.lines,
                     disp.color.bar = disp.color.bar,
                     disp.phys.dist = disp.phys.dist,
                     image.title = image.title,
                     disp.legend = disp.legend,
                     palette.file = palette.file,
                     image.type = image.type,
                     image.size = image.size,
                     disp.mult.lab.X = disp.mult.lab.X,
                     image.name = image.name,
                     print.image = print.image,
                     connecting.lines.factor = connecting.lines.factor,
                     connecting.lines.adj = connecting.lines.adj,
                     connecting.lines.vert.adj = connecting.lines.vert.adj,
                     connecting.lines.flex = connecting.lines.flex,
                     font.factor = font.factor,
                     symbol.factor = symbol.factor,
                     verbose = verbose)

  #-------------------CONFIGURATION VARIABLES BEGINS from config file---------
  gbl.var$verbose <- config.var$verbose

  if(!is.null(config.file)) {
    config.var <- read.config(config.file, config.var)
  }

  #-------------- CHECK if all parameters have values -----------------
  if(is.null(config.var$mydata.file)) {
    stop("Invalid MYDATA data file: ", config.var$mydata.file, "\n")
  }
  check.configVar(config.var)


  #------------- Track define outside
  if(!is.null(tracks.gviz) ){
    gbl.var$listtracks_gviz <- tracks.gviz
  }

  if(!is.null(tracks.ggbio) ){
    gbl.var$listtracks_ggbio <- tracks.ggbio
  }

  if(!is.null(tracks.trackviewer) ){
    gbl.var$listtracks_trackviewer <- tracks.trackviewer
  }

  #------------- connection to database
  #  if(!is.null(config.var$genome) & !is.null(config.var$BROWSER.SESSION)){
  #    mySession <- browserSession(config.var$BROWSER.SESSION)
  #    genome(mySession) <- config.var$genome
  #    gbl.var$mySession <- mySession
  #  }

  #------------- READ DATA for ANNOTATION TRACKS
  if(!is.null(config.var$biofeat.user.file)){
    split.biofeature.data.user.file <- strsplit(config.var$biofeat.user.file, ",")
    split.biofeature.data.user.type <- strsplit(config.var$biofeat.user.type, ",")

    gbl.var$split.biofeature.data.user.file <- split.biofeature.data.user.file
    gbl.var$split.biofeature.data.user.type <- split.biofeature.data.user.type

    if(!is.null(config.var$biofeat.user.type.plot)){
      split.biofeature.data.user.type.plot <- strsplit(config.var$biofeat.user.type.plot, ",")
      gbl.var$split.biofeature.data.user.type.plot <- split.biofeature.data.user.type.plot
    }
  }

  #------------- READ DATA and UPDATE VARIABLES
  if(!is.null(config.var$cormatrix.file)){
    if(!is.null(config.var$cormatrix.type)){
      split.cormatrix.type <- strsplit(config.var$cormatrix.type, ",")
      comatrix.type.length <- length(split.cormatrix.type[[1]])
      gbl.var$split.cormatrix.type <- split.cormatrix.type
    } else {
      stop("Need to define the format of your correlation matrix from FILE or from MATRIX")
    }
    if(gbl.var$split.cormatrix.type == "listfile"){
      split.cormatrix.file <- strsplit(config.var$cormatrix.file, ",")
      comatrix.file.length <- length(split.cormatrix.file[[1]])
      gbl.var$split.cormatrix.file <- split.cormatrix.file
    }else {
      split.cormatrix.file <- config.var$cormatrix.file
      comatrix.file.length <- length(split.cormatrix.file)
      gbl.var$split.cormatrix.file <- split.cormatrix.file
    }

  } else {
    warning("No visualisation of correlation matrice")
    disp.mydata <- FALSE
  }

  if(!is.null(config.var$mydata.file)){
    #--- TYPE of DATA
    if(!is.null(config.var$mydata.type)){
      split.type <- strsplit(config.var$mydata.type, ",")
      split.type.length <- length(split.type[[1]])
      gbl.var$split.type <- unlist(split.type)
      if(config.var$verbose){
        cat("TYPE DATA:",gbl.var$split.type,"\n")
      }
    } else{
      stop("Need to define the type of data: FILE, MATRIX\n")
    }

    mydata.file.length <- 0
    if(config.var$mydata.type == "file") {
      split.mydata.file <- strsplit(config.var$mydata.file, ",")
      mydata.file.length <- length(split.mydata.file[[1]])

      gbl.var$split.mydata.file <- split.mydata.file
      gbl.var$mydata.samples <- mydata.file.length
    } else {
      split.mydata.file <- list(config.var$mydata.file)
      mydata.file.length <- length(split.mydata.file)

      gbl.var$split.mydata.file <- split.mydata.file
      gbl.var$mydata.samples <- mydata.file.length
    }

    gbl.var$hap.samples <- 0
    samples <- mydata.file.length
    gbl.var$samples <- samples

    #REFERENCE DATA
    split.color.list <- create.color.list(config.var,gbl.var)
    gbl.var$color.list <- split.color.list[[1]]

    split.symbol.fill.lists <- create.symbol.list(config.var, split.color.list, gbl.var)

    gbl.var$symbol.list <- split.symbol.fill.lists$split.symbol.list
    gbl.var$fill.list <- split.symbol.fill.lists$split.fill.list[[1]]


    #create sample labels
    #split sample.labels on comma, then take the top 5 or less entries and create a list
    if(!is.null(config.var$sample.labels)) {
      split.tmp.sample.list <- strsplit(config.var$sample.labels, ",")
    } else {
      split.tmp.sample.list <- rep("sample", samples)
    }

    if(length(split.tmp.sample.list[[1]]) < samples) {
      warning("sample.labels contains less labels than the number of samples. Labels must be separated by commas without spaces.\n")
    }

    #--- FORMAT of DATA
    if(!is.null(config.var$mydata.format)){
      split.format <- strsplit(config.var$mydata.format, ",")
      split.format.length <- length(split.format[[1]])
      gbl.var$split.format <- split.format
    }else{
      stop("Need to define the format\n")
    }


    #-- Visualisation of ASSOCIATION
    if(!is.null(config.var$disp.association)){
      split.association <- strsplit(config.var$disp.association, ",")
      split.association.length <- length(split.association[[1]])
      gbl.var$split.association <- split.association
    }else{
      stop("Need to define if there is the visualisation of association\n")
    }

    if(!is.null(config.var$disp.beta.association)){
      split.beta.association <- strsplit(config.var$disp.beta.association, ",")
      split.beta.association.length <- length(split.beta.association[[1]])
      gbl.var$split.beta.association <- split.beta.association
      gbl.var$myfactor.beta <- config.var$factor.beta
    }else{
      stop("Need to define if there is the visualisation of association\n")
    }

    #---VISUALISATION of region
    if(!is.null(config.var$disp.region)){
      split.region <- strsplit(config.var$disp.region, ",")
      split.region.length <- length(split.region[[1]])
      gbl.var$split.region <- split.region
    }else{
      stop("Need to define if there is the visualisation of region\n")
    }

    split.tmp.sample.list.truncated <- substr(split.tmp.sample.list[[1]], 1, 12)

  } else {
    stop("No data visualise")
  }

  #DEBUG STATEMENT
  # if (config.var$verbose)  cat("samples ", samples, "\n")

  gbl.var$split.sample.labels <- list(head(split.tmp.sample.list.truncated, samples))

  #---LARGE REFERENCE DATA
  if(!is.null(config.var$mydata.large.file)){
    #--- TYPE of LARGE DATA
    if(!is.null(config.var$mydata.large.type)){
      large.split.type <- strsplit(config.var$mydata.large.type, ",")
      large.split.type.length <- length(large.split.type[[1]])
      gbl.var$large.split.type <- large.split.type
    }else{
      stop("Need to define the type of extra data (LISTFILE,listdataframe)\n")
    }
    large.mydata.file.length <- 0
    if(gbl.var$large.split.type == "listfile"){
      large.split.mydata.file <- strsplit(config.var$mydata.large.file, ",")
      large.mydata.file.length <- length(large.split.mydata.file[[1]])

      gbl.var$large.split.mydata.file <- large.split.mydata.file
      gbl.var$large.mydata.samples <- large.mydata.file.length
    } else {
      large.split.mydata.file <- config.var$mydata.large.file
      large.mydata.file.length <- length(large.split.mydata.file)

      gbl.var$large.split.mydata.file <- large.split.mydata.file
      gbl.var$large.mydata.samples <- large.mydata.file.length
    }

    gbl.var$large.hap.samples <- 0
    large.samples <- large.mydata.file.length
    gbl.var$large.samples <- large.samples

    large.split.color.list <- create.color.list.large(config.var,gbl.var)
    gbl.var$large.color.list <- large.split.color.list[[1]]

    large.split.symbol.fill.lists <- create.symbol.list.large(config.var, large.split.color.list, gbl.var)

    gbl.var$large.symbol.list <- large.split.symbol.fill.lists$large.split.symbol.list
    gbl.var$large.fill.list <- large.split.symbol.fill.lists$large.split.fill.list[[1]]


    #create sample labels
    #split sample.labels on comma, then take the top 5 or less entries and create a list


    if(!is.null(config.var$sample.labels.large)) {
      large.split.tmp.sample.list <- strsplit(config.var$sample.labels.large, ",")
    }

    if(length(large.split.tmp.sample.list[[1]]) < large.samples) {
      warning("LARGE.sample.labels contains less labels than the number of samples. Labels must be separated by commas without spaces.\n")
    }

    large.split.tmp.sample.list.truncated <- substr(large.split.tmp.sample.list[[1]], 1, 12)

    #DEBUG STATEMENT
    # if (config.var$verbose)  cat("large.samples ", large.samples, "\n")

    gbl.var$large.split.sample.labels <- list(head(large.split.tmp.sample.list.truncated, large.samples))

    #--- FORMAT of LARGE DATA
    if(!is.null(config.var$mydata.large.format)){
      large.split.format <- strsplit(config.var$mydata.large.format, ",")
      large.split.format.length <- length(large.split.format[[1]])
      gbl.var$large.split.format <- large.split.format
    }else{
      stop("Need to define the format of extra data\n")
    }

    #-- Visualisation of ASSOCIATION
    if(!is.null(config.var$disp.association.large)){
      large.split.association <- strsplit(config.var$disp.association.large, ",")
      large.split.association.length <- length(large.split.association[[1]])
      gbl.var$large.split.association <- large.split.association
    } else{
      stop("Need to define the visualisation of association of extra data\n")
    }

    # -- Visualisation of ASSOCIATION for beta
    if(!is.null(config.var$disp.beta.association.large)){
      large.split.beta.association <- strsplit(config.var$disp.beta.association.large, ",")
      large.split.beta.association.length <- length(large.split.beta.association[[1]])
      gbl.var$large.split.beta.association <- large.split.beta.association
      gbl.var$myfactor.beta <- config.var$factor.beta
    } else{
      stop("Need to define the visualisation of association of extra data\n")
    }

    #---VISUALISATION of region
    if(!is.null(config.var$disp.region.large)){
      large.split.region <- strsplit(config.var$disp.region.large, ",")
      large.split.region.length <- length(large.split.region[[1]])
      gbl.var$large.split.region <- large.split.region
    }else{
      stop("Need to define if there is the visualisation of region of extra data\n")
    }

  } else {
    gbl.var$large.split.mydata.file <- NULL
    gbl.var$large.mydata.samples <- NULL
    gbl.var$large.hap.samples <- 0
    gbl.var$large.samples <- NULL
    gbl.var$large.color.list <- NULL
    gbl.var$large.symbol.list <- NULL
    gbl.var$large.fill.list <- NULL
    gbl.var$large.split.sample.labels <- NULL
    gbl.var$large.split.sample.format <- NULL
    gbl.var$large.split.sample.type <- NULL
    gbl.var$large.split.association <- NULL
    gbl.var$large.split.beta.association <- NULL
    gbl.var$large.split.region <- NULL
  }

  #---- FIX GENERAL VARIABLES


  fix.var <- retrieve.data(config.var, gbl.var)
  config.var <- fix.var$config.var
  gbl.var <- fix.var$gbl.var


  #############DRAW Picture

  #---------------- DRAW DIFFERENT ANNOTATION TRACK ---------
  # gbl.var <- create.tracks.web(config.var,gbl.var)

  if( disp.pvalueplot == FALSE & disp.cormatrixmap == FALSE){
    stop("Impossible to have no draw pvalue plot and correlation matrix!\n
          Change the value of disp.pvalueplot and disp.cormatrixmap")
  }

  #------ DRAW the STRUCTURE COMET
  if(disp.pvalueplot) {
    if(print.image == FALSE || is.null(config.var$image.name)){
      gbl.var <- draw.plot.comet(config.var, gbl.var,newpage=TRUE)
    } else{
      printPlot.comet(config.var, gbl.var)
    }
  }else{
    if(print.image == FALSE || is.null(config.var$image.name)){
      gbl.var <- draw.plot.comet.nopval(config.var, gbl.var,newpage=TRUE)
    } else{
      printPlot.comet.nopval(config.var, gbl.var)
    }
  }
  #END FUNCTION LIST

  invisible(list(config.var, gbl.var))
  #DEBUG STATEMENT
  if (config.var$verbose)  cat("FINISH GENERIC COMET \n")
}