# 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/"
# }
# NOT RUN {
# a) provide the seed nodes/genes with the weight info
## load ImmunoBase
ImmunoBase <- xRDataLoader(RData.customised='ImmunoBase',
RData.location=RData.location)
## get genes within 500kb away from AS GWAS lead SNPs
seeds.genes <- ImmunoBase$AS$genes_variants
## seeds weighted according to distance away from lead SNPs
data <- 1- seeds.genes/500000
# b) prepare a data frame
df <- data.frame(name=names(data), value=data, stringsAsFactors=FALSE)
# c) do correlation
gp <- xPolarDot(df[1:50,])
gp
# }
Run the code above in your browser using DataLab