lapply和BioMart的另一个问题

Question

我最近问了一个关于lapply的问题，但现在我改变了方法，遇到了更多我无法解决的问题。

这是我遇到问题的代码：

加载biomaRt

library(biomaRt)

让hsapiens mart

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)

ensembl_gene_ID <- hsapien_PC_genes$ensembl_gene_id

设置物种向量

species <- c("mmusculus", "ggallus")

这一切都像预期的那样工作，我在使用lapply时遇到了问题

all_homologues <- lapply(species, function(s) getBM(attributes = c("ensembl_gene_id", 
                                                               "external_gene_name", 
                                                               get(paste0(s, "_homolog_ensembl_gene")), 
                                                               get(paste0(s, "_homolog_associated_ensembl_gene")), 
                                                               filters = "ensembl_gene_id", 
                                                               values = c(ensembl_gene_ID), 
                                                               mart = ensembl_hsapiens)))

我得到了错误消息：

 Error in martCheck(mart) : 
  You must provide a valid Mart object. To create a Mart object use the function: useMart.  Check ?useMart for more information.

代码的工作方式如下所示：

all_homologues <- getBM(attributes = c("ensembl_gene_id", "external_gene_name",  
                               "mmusculus_homolog_ensembl_gene", 
                               "mmusculus_homolog_associated_gene_name"),
                           filters = "ensembl_gene_id",
                           values = c(ensembl_gene_ID),
                           mart = ensembl_hsapiens)

但是我有许多物种，所以我希望能够使用一小段代码来提取物种的同源，而不是必须为每个物种重写所有同源。

我认为我用lapply做错了什么，因为在不使用lapply的情况下，mart可以很好地工作。

我之前的问题：Issue with lapply using biomart

Answer 1

在尝试过它之后，我通过使用以下命令设法让它正常工作：

all_homologues <- lapply(species, function(s) getBM(attributes = c("ensembl_gene_id", 
                                                               "external_gene_name", 
                                                               paste0(s, c("_homolog_ensembl_gene",
                                                                           "_homolog_associated_gene_name"))),
                                                filters = "ensembl_gene_id",
                                                values = c(ensembl_gene_ID),
                                                mart = ensembl_hsapiens))

虽然我不是很确定为什么？