Below is a small fraction of my code:
library(biomaRt)
ensembl_hsapiens <- useMart("ensembl",
dataset = "hsapiens_gene_ensembl")
hsapien_PC_genes <- getBM(attributes = c("ensembl_gene_id", "external_gene_name"),
filters = "biotype",
values = "protein_coding",
mart = ensembl_hsapiens)
paralogues[["hsapiens"]] <- getBM(attributes = c("external_gene_name",
"hsapiens_paralog_associated_gene_name"),
filters = "ensembl_gene_id",
values = c(ensembl_gene_ID) , mart = ensembl_hsapiens)
This bit of code will only allow me to extract the paralogues for hsapiens, it there a way for me to easily get the same information for mmusculus (mouse) and ggallus (chicken) without having to rewrite the code by using something like Tapply? My code is much longer than the snippet provided, all I would need to do is swap the word hsapiens for mmusulus and ggallus.
See Question&Answers more detail:
os 与恶龙缠斗过久,自身亦成为恶龙;凝视深渊过久,深渊将回以凝视…