前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >plasmapR优雅的可视化质粒图谱

plasmapR优雅的可视化质粒图谱

作者头像
R语言数据分析指南
发布2023-08-18 13:38:27
3210
发布2023-08-18 13:38:27
举报

欢迎关注R语言数据分析指南

❝本节来介绍一款R包「plasmapR」提供解析和绘制.gb质粒文件的功能,一旦质粒以Genbank格式导出就可以对其进行解析和绘图。感兴趣的朋友欢迎分享转发,「更多详细内容请参考作者官方文档」

官方文档

❝https://github.com/BradyAJohnston/plasmapR ❞

加载R包

代码语言:javascript
复制
devtools::install_github("bradyajohnston/plasmapr")  # 安装plasmapR包

library(plasmapR)  # 载入plasmapR包
library(tidyverse)  # 载入tidyverse包

导入系统数据

代码语言:javascript
复制
fl <- system.file('extdata', 'petm20.gb', package = "plasmapR")  

案例 1

代码语言:javascript
复制
fl %>% read_gb() %>% plot_plasmid(name = "pETM-20")  # 读取.gb文件并绘制质粒图,指定名称为"pETM-20"

案例 2

代码语言:javascript
复制
fl <- system.file('extdata', 'petm20.gb', package = "plasmapR")  # 获取文件路径
代码语言:javascript
复制
plasmid <- fl %>% read_gb()  # 读取.gb文件并将结果保存为plasmid对象
dat <- plasmid %>% as.data.frame()  # 将plasmid对象转换为数据框
代码语言:javascript
复制
dat[dat$type == "CDS", ] %>% 
  plot_plasmid(name = "pETM-20")  # 选择type为"CDS"的行数据并绘制质粒图,指定名称为"pETM-20"

案例3 修改主题

代码语言:javascript
复制
fl <- system.file('extdata', '20.gb', package = "plasmapR")  # 获取文件路径

plt <- fl %>% read_gb() %>% 
  plot_plasmid()  # 读取.gb文件并绘制质粒图,结果保存为plt对象
代码语言:javascript
复制
plt + ggplot2::theme_bw() +
  theme(axis.text=element_blank(),
        axis.ticks = element_blank(),
        axis.title.x = element_blank(),
        axis.title.y = element_blank())  # 修改图表主题,隐藏轴标签和刻度
本文参与 腾讯云自媒体同步曝光计划,分享自微信公众号。
原始发表:2023-06-22,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 R语言数据分析指南 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 欢迎关注R语言数据分析指南
  • 官方文档
  • 加载R包
  • 导入系统数据
  • 案例 1
  • 案例 2
  • 案例3 修改主题
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档