library(pacman)pacman::p_load(genoPlotR,dplyr,magrittr,stringi);options(warn=-1)data-"data.xls"%%ad.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(iinseq_along(list)){list[]%%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()
数据
预览时标签不可点收录于话题#个上一篇下一篇