#' 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.