results_log = fread("result.assoc.txt")
dim(results_log)
head(results_log)
summary(results_log)

select = dplyr::select
d1 = results_log %>% select(Chr = chr,Marker = rs,Pos = ps,p=p_wald)
d1 = d1 %>% drop_na(p)
summary(d1)

manhattan(d1,chr="Chr",bp="Pos",p="p",snp="Marker", main = "Manhattan plot: logistic")
tiff("y1-曼哈頓圖.tiff")
manhattan(d2,chr="Chr",bp="Pos",p="p",snp="Marker", main = "Manhattan plot",suggestiveline = 5.237,genomewideline = 6.538,annotatePval = 0.00000029, col = c("#F07E4A", "#3985F0"))

qq(d1$p, main = "Q-Q plot of GWAS p-values : log")
tiff("y1-QQ圖.tiff")
qq(d1$p, main = "Q-Q plot of GWAS p-values : log")
dev.off()

你是否還在尋找穩(wěn)定的海外服務器提供商？創(chuàng)新互聯(lián)www.cdcxhl.cn海外機房具備T級流量清洗系統(tǒng)配攻擊溯源，準確流量調(diào)度確保服務器高可用性，企業(yè)級服務器適合批量采購，新人活動首月15元起，快前往官網(wǎng)查看詳情吧

成都創(chuàng)新互聯(lián)主打移動網(wǎng)站、網(wǎng)站制作、成都網(wǎng)站設計、網(wǎng)站改版、網(wǎng)絡推廣、網(wǎng)站維護、域名注冊、等互聯(lián)網(wǎng)信息服務，為各行業(yè)提供服務。在技術實力的保障下，我們?yōu)榭蛻舫兄Z穩(wěn)定，放心的服務，根據(jù)網(wǎng)站的內(nèi)容與功能再決定采用什么樣的設計。最后，要實現(xiàn)符合網(wǎng)站需求的內(nèi)容、功能與設計，我們還會規(guī)劃穩(wěn)定安全的技術方案做保障。
當前文章：R曼哈頓圖QQ圖-創(chuàng)新互聯(lián)
本文URL：http://weahome.cn/article/iehjd.html