###########################################
## Compute information content of each HP according to associated genes
data(geneByHp, hp_descendants, package="PCAN")
geneByHp <- unstack(geneByHp, entrez~hp)
ic <- computeHpIC(geneByHp, hp_descendants)
hist(
ic,
breaks=100, col="grey",
main="Distribution of Information Content",
xlab="IC base on genes associated to HP"
)
Run the code above in your browser using DataLab