Showing content from https://www.bioconductor.org/packages/devel/workflows/html/../vignettes/liftOver/inst/doc/liftov.R below:
## ----echo=FALSE, results="hide", warning=FALSE-------------------------------- suppressPackageStartupMessages({ library('liftOver') }) ## ----doit,echo=FALSE,results="hide"------------------------------------------- library(gwascat) library(GenomicRanges) library(rtracklayer) library(Homo.sapiens) library(BiocGenerics) library(liftOver) ## ----lkOne,eval=FALSE--------------------------------------------------------- # library(gwascat) # cur = makeCurrentGwascat() # result varies by day ## ----lkcur-------------------------------------------------------------------- data(cur) cur ## ----getch-------------------------------------------------------------------- library(rtracklayer) path = system.file(package="liftOver", "extdata", "hg38ToHg19.over.chain") ch = import.chain(path) ch str(ch[[1]]) ## ----dolift------------------------------------------------------------------- seqlevelsStyle(cur) = "UCSC" # necessary cur19 = liftOver(cur, ch) class(cur19) ## ----ul----------------------------------------------------------------------- cur19 = unlist(cur19) genome(cur19) = "hg19" cur19 = new("gwaswloc", cur19) cur19 ## ----lkloss------------------------------------------------------------------- length(cur)-length(cur19) setdiff(mcols(cur)$SNPS, mcols(cur19)$SNPS)
RetroSearch is an open source project built by @garambo
| Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4