## ---- echo=FALSE, results="hide"，警告= FALSE ------------------------------- (“liftOver suppressPackageStartupMessages({库 ') }) ## ---- doit,回声= FALSE,结果=“隐藏 "------------------------------------------- 库(gwascat)库(GenomicRanges)图书馆(rtracklayer)图书馆(Homo.sapiens)图书馆(BiocGenerics)图书馆(liftOver) # #——lkOne eval = FALSE --------------------------------------------------------- # 库(gwascat) # cur = makeCurrentGwascat() # # #结果不同的一天——lkcur -------------------------------------------------------------------- 数据(坏蛋)cur # #——getch -------------------------------------------------------------------- 图书馆(rtracklayer)路径=系统。文件(包=“liftOver”、“extdata”、“hg38ToHg19.over.chain”)ch = import.chain(路径)ch str (ch [[1 ]]) ## ---- dolift ------------------------------------------------------------------- seqlevelsStyle(坏蛋)=“UCSC”#必要cur19 = liftOver (cur, ch)类(cur19) # # - ul ----------------------------------------------------------------------- cur19 = unlist (cur19)基因组(cur19) =“hg19”cur19 = new(“gwaswloc”,cur19) cur19 # #——lkloss ------------------------------------------------------------------- 长度(坏蛋)长度(cur19) setdiff (mcols(坏蛋)美元SNPS, mcols (cur19)美元SNPS)