繪制基因結(jié)構(gòu)圖的方法有很多路鹰，TBtools,GSDS網(wǎng)站等;今天介紹如何通過R繪制基因結(jié)構(gòu)圖贷洲，廢話不多說直接看代碼。
library(pacman)
pacman::p_load(genoPlotR, dplyr, magrittr,
stringi);options(warn = -1)
data <- "data.xls" %>% read.delim() %>%
arrange(molecule)
list <- data %>% mutate(strand = 1) %>%
select(1, 3, 4, 5) %>%
set_colnames(c("name", "start", "end", "strand")) %>%
group_split(name) %>% as.list() %>%
set_names(unique(data[, 1]))
for (i in seq_along(list)) { list[[i]] %<>%
as.dna_seg(col="black",fill="steelblue") }
list[data[1, 1] %>% stri_sub(length = 5) %>%
paste0(., (list %>% names() %>%
length() %>% seq_len())) %>%
match(names(list))] %>% plot_gene_map()
gene.jpeg
只需要安裝加載所需要的R包晋柱，讀入如下示例格式數(shù)據(jù)即可,
file.png
