安裝并加載所需R包
# install.packages("qqman")
library(qqman)
# 加載并查看示例數(shù)據(jù)
head(gwasResults)
SNP CHR BP P
1 rs1 1 1 0.9148060
2 rs2 1 2 0.9370754
3 rs3 1 3 0.2861395
4 rs4 1 4 0.8304476
5 rs5 1 5 0.6417455
6 rs6 1 6 0.5190959
基本數(shù)據(jù)格式：SNP名稱队伟，所在染色體抠艾，SNP位置低散，Pvalue值
繪制曼哈頓圖
基本使用
manhattan(x, chr = "CHR", bp = "BP", p = "P", snp = "SNP",
col = c("gray10", "gray60"), chrlabs = NULL,
suggestiveline = -log10(1e-05), genomewideline = -log10(5e-08),
highlight = NULL, logp = TRUE, annotatePval = NULL,
annotateTop = TRUE, ...)
常用參數(shù)
x 輸入數(shù)據(jù)文件
chr 設置染色體所在的列
bp 設置SNP位置所在的列
p 設置Pvalue值所在的列
snp 設置SNP名稱所在的列
col 設置不同染色體上點的顏色
chrlabs 設置染色體的標簽名
suggestiveline 設置"suggestive"線的閾值，默認為-log10(1e-5)
genomewideline 設置"genome-wide significant"線的閾值缎浇，默認為-log10(5e-8)
highlight 高亮感興趣的SNP位點
logp 設置是否對Pvalue去log10值
annotatePval 注釋低于指定p-value閾值的位點
annotateTop 是否只注釋低于指定p-value閾值位點中的top hit
默認繪圖
manhattan(gwasResults)
image.png
調(diào)整相關參數(shù)
manhattan(gwasResults, main = "Manhattan Plot", ylim = c(0, 10), cex = 0.6, cex.axis = 0.9,
col = c("blue4", "orange3"), suggestiveline = F, genomewideline = F,
chrlabs = c(1:20, "P", "Q"))
image.png
繪制單條染色體
manhattan(subset(gwasResults, CHR == 1))
image.png
高亮感興趣的SNP位點
head(snpsOfInterest)
## [1] "rs3001" "rs3002" "rs3003" "rs3004" "rs3005" "rs3006"
manhattan(gwasResults, highlight = snpsOfInterest)
image.png
只注釋低于指定p-value閾值位點中的top hit
manhattan(gwasResults, annotatePval = 0.01, annotateTop=TRUE)
image.png
注釋低于指定p-value閾值的所有位點
manhattan(gwasResults, annotatePval = 0.005, annotateTop = FALSE)
image.png
繪制QQ圖
默認繪圖
qq(gwasResults$P)
image.png
調(diào)整參數(shù)
qq(gwasResults$P, main = "Q-Q plot of GWAS p-values", xlim = c(0, 7),
ylim = c(0, 12), pch = 18, col = "blue4", cex = 1.5, las = 1)
image.png
sessionInfo()
## R version 3.5.1 (2018-07-02)
## Platform: x86_64-apple-darwin15.6.0 (64-bit)
## Running under: OS X El Capitan 10.11.3
##
## Matrix products: default
## BLAS: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRblas.0.dylib
## LAPACK: /Library/Frameworks/R.framework/Versions/3.5/Resources/lib/libRlapack.dylib
##
## locale:
## [1] zh_CN.UTF-8/zh_CN.UTF-8/zh_CN.UTF-8/C/zh_CN.UTF-8/zh_CN.UTF-8
##
## attached base packages:
## [1] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] qqman_0.1.4
##
## loaded via a namespace (and not attached):
## [1] compiler_3.5.1 backports_1.1.2 magrittr_1.5 rprojroot_1.3-2
## [5] tools_3.5.1 htmltools_0.3.6 yaml_2.2.0 Rcpp_0.12.18
## [9] calibrate_1.7.2 stringi_1.2.4 rmarkdown_1.10 knitr_1.20
## [13] stringr_1.3.1 digest_0.6.16 evaluate_0.11
