仿GWAS曼哈顿图

1. 根据vcf文件获取一个类似曼哈顿图的input文件,只是第4列的p value 值变为了区间的snp number

如下:

SNP     CHR     BP      P
1       chr1    100000  1362
2       chr1    110000  1319
3       chr1    120000  1321
4       chr1    130000  1426
5       chr1    140000  1500
6       chr1    150000  1512
7       chr1    160000  1387
8       chr1    170000  1235
9       chr1    180000  1244
10      chr1    190000  1221
###########################

2. 使用R进行绘图,注意使用facet_grid 

library(ggplot2)
library("RColorBrewer")
#########################################
dat <- read.table("snp_manhattan.txt",header = TRUE, sep="\t")

pdf("snp.pdf")
ggplot(dat,aes(x=BP,y=P,colour=CHR)) + facet_grid(.~CHR,scales="free_x",space="free_x") +
  geom_point() + xlab("Position") + ylab("Numb

  • 0
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值