GCdepth散点图绘制

一 背景

基因组组装完成后,需要对组装结果进行评估,其中GC_depth图是一个比较重要的指标。该图的横坐标是GC含量,纵坐标是平均深度。以30k为窗口(window)无overlap计算其GC含量和平均深度。如果存在样品污染,通常能够从GC含量分析中呈现出来。
下文探索了几种绘制GC_depth散点图的方法思路。

二 具体方法

2.1 直接生成散点图

由于数据点太多,且相互重合,弱化了想要传达给读者的重要信息。

df <- read.table("result.txt", header = TRUE)
png("GC_depth.png", width = 20, height = 18, units = "cm", res = 300)
plot(df, ylab="Average depth (X)", xlab="GC content (%)", cex.lab = 1.4, cex.axis = 1.3, pch = 20, ylim = c(0,100), xlim = c(10,80))
dev.off()

2.2 通过透明度展示点密度

可能是点太多,调试了多次无法达到理想的效果。

df <- read.table("result.txt", header = TRUE)
MyGray <- rgb(t(col2rgb("black")), alpha=50, maxColorValue=255)
png("GC_depth_alpha.png", width = 20, height = 18, units = "cm", res = 300)
plot(df, ylab="Average depth (X)", xlab=<
  • 5
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 12
    评论
评论 12
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值