## Not run:
# # Load the XGR package and specify the location of built-in data
# library(XGR)
# RData.location <- "http://galahad.well.ox.ac.uk/bigdata_dev"
#
# # a) provide the genomic regions
# ## load ImmunoBase
# ImmunoBase <- xRDataLoader(RData.customised='ImmunoBase',
# RData.location=RData.location)
# ## get lead SNPs reported in AS GWAS and their significance info (p-values)
# gr <- ImmunoBase$AS$variant
# df <- as.data.frame(gr, row.names=NULL)
# chr <- df$seqnames
# start <- df$start
# end <- df$end
# data <- paste(chr,':',start,'-',end, sep='')
#
# # b) create a GRanges object
# GR <- xGR(data=data, format="chr:start-end",
# RData.location=RData.location)
# ## End(Not run)
Run the code above in your browser using DataLab