# #——风格,呼应= FALSE,结果=“黑名单”- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - BiocStyle::减价()# #——lkd,缓存= TRUE - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - suppressPackageStartupMessages(库(ldblock)) = dir(道路系统。文件(“人类基因组单体型图”,包=“ldblock”),全= TRUE) ceu17 = hmld(路径,poptag =“CEU”,铬=“chr17”) ceu17 # #——abc,无花果= TRUE, fig.width = 7, fig.height = 7 - - - - - - - - - - - - - - - - - - - - - - - - - - - - -库(矩阵)图像(ceu17@ldmat [1:400 1:400], col.reg = heat.colors (120), colorkey = TRUE, useRaster = TRUE) # #——getg - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -库(gwascat)数据(ebicat37) seqlevelsStyle (ebicat37) =“NCBI”e17 = ebicat37 ((as.character (seqnames (ebicat37)) = =“17”)] # #——getrs - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - rsh17 =独特(e17 snp美元)头(rsh17) # #——doexpa - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -长度(rsh17) exset = expandSnpSet (rsh17, ldstruct = ceu17磅=。9)长度(exset)所有(rsh17 % % exset)