workflows/210127_counts_of_completely_ko_genes.R

tab <- read.delim('http://compbio.charite.de/jenkins/job/hpo.annotations/lastSuccessfulBuild/artifact/util/annotation/genes_to_phenotype.txt', skip = 1)
cols = unlist(strsplit("entrez-gene-id<tab>entrez-gene-symbol<tab>HPO-Term-Name<tab>HPO-Term-ID<tab>Frequency-Raw<tab>Frequency-HPO<tab>Additional Info from G-D source<tab>G-D source<tab>disease-ID for link", "<tab>"))
colnames(tab) <- cols

tab[tab$`HPO-Term-ID` == 'Autosomal recessive inheritance',]

genes = unique(tab[tab$`HPO-Term-ID` == 'Autosomal recessive inheritance',]$`entrez-gene-symbol`)
frhl/our documentation built on Feb. 5, 2021, 7:30 p.m.