data(geneByTrait, traitDef, geneDef, package="PCAN")
omim <- "612285"
traitDef[which(traitDef$id==omim),]
# Gene associated to this disease
entrez <- geneByTrait[which(geneByTrait$id==omim), "entrez"]
geneDef[which(geneDef$entrez %in% entrez),]
# All diseases associated to this gene
traitDef[
which(
traitDef$id %in%
geneByTrait[which(geneByTrait$entrez==entrez), "id"]
),
]
Run the code above in your browser using DataLab