BioMart: Is there a way to easily change the species for all of my code?

后端 未结 1 1490
孤独总比滥情好
孤独总比滥情好 2020-12-22 02:40

Below is a small fraction of my code:

library(biomaRt)

ensembl_hsapiens <- useMart(\"ensembl\", 
                        dataset = \"hsapiens_gene_ensem         


        
相关标签:
1条回答
  • 2020-12-22 02:52

    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)

    0 讨论(0)
提交回复
热议问题