Below is a small fraction of my code:
library(biomaRt)
ensembl_hsapiens <- useMart(\"ensembl\",
dataset = \"hsapiens_gene_ensem
Easy way is to just wrap it all in a for
loop:
library(biomaRt)
species <- c("hsapiens_gene_ensembl", "mmusculus_gene_ensembl", "ggallus_gene_ensembl")
for (s in species) {
tmp <- useMart("ensembl", dataset = paste0(s))
hsapien_PC_genes <- getBM(attributes = c("ensembl_gene_id", "external_gene_name"),
filters = "biotype",
values = "protein_coding",
mart = tmp)
paralogues[[s]] <- getBM(attributes = c("external_gene_name",
"hsapiens_paralog_associated_gene_name"),
filters = "ensembl_gene_id",
values = c(ensembl_gene_ID) , mart = tmp)
}
This should work, I haven't tested it because I don't have those packages installed. I've changed the names of some vars to make more sense (eg tmp
)