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()

你是否还在寻找稳定的海外服务器提供商？创新互联www.cdcxhl.cn海外机房具备T级流量清洗系统配攻击溯源，准确流量调度确保服务器高可用性，企业级服务器适合批量采购，新人活动首月15元起，快前往官网查看详情吧

成都创新互联主打移动网站、网站制作、成都网站设计、网站改版、网络推广、网站维护、域名注册、等互联网信息服务，为各行业提供服务。在技术实力的保障下，我们为客户承诺稳定，放心的服务，根据网站的内容与功能再决定采用什么样的设计。最后，要实现符合网站需求的内容、功能与设计，我们还会规划稳定安全的技术方案做保障。
文章标题：R曼哈顿图QQ图-创新互联
标题路径：http://chengdu.cdxwcx.cn/article/iehjd.html