R/mapa_chr.R

In charo1999/mapagenetica:

#' mapa_chr
#'
#'
#' @param cr chr vector
#' @param col_cr color
#' @param unidad cM
#'
#' @export
#'
mapa_chr <- function(cr= c(30.4,19.7,23.5,18.6,27), col_cr = "grey", unidad = "cM") {
  rango <- 1:length(cr)
  plot(x = rango, # cuantas cromosomas
       y = cr, # longitud de cada uno
       ylim = range(0, max(cr)),
       type = "h", # tipo: histograma
       xaxt = "n", # no poner las etiquetas en x
       yaxt = "n",
       bty = "n", # sin caja
       col = "white",
       xlab = "Cromosomas",
       ylab = unidad
  )
  axis(side = 2,
       at = pretty(0:max(cr)),
       labels = rev(pretty(0:max(cr))),
       lty = 1
  )
  ###ooiijmmnnnnnn
  axis(side = 3,
       at = rango,
       labels = paste0("LG", rango),
       lty = 0
  )
  segments(x0 = rango,
           y0 = abs(max(cr) - cr),
           x1 = rango,
           y1 = max(cr),
           lwd = 10,
           col = col_cr
  )
}

a_th <- list(crm = c(30.4, 19.7, 23.5, 18.6, 27),
             col = "grey",
             uni = "Mbp")

mapa_chr(cr = a_th$crm, col_cr = a_th$col, unidad = a_th$uni)