lapply和BioMart的另一个问题
我最近问了一个关于lapply的问题,但现在我改变了我的方法,遇到了更多我无法解决的问题 这是我遇到问题的代码: 装载生物零件lapply和BioMart的另一个问题,r,bioinformatics,lapply,bioconductor,biomart,R,Bioinformatics,Lapply,Bioconductor,Biomart,我最近问了一个关于lapply的问题,但现在我改变了我的方法,遇到了更多我无法解决的问题 这是我遇到问题的代码: 装载生物零件 library(biomaRt) 做海皮斯超市 ensembl_hsapiens <- useMart("ensembl", dataset = "hsapiens_gene_ensembl") 代码的工作原理如下: all_homologues <- getBM(attributes = c("ense
library(biomaRt)
做海皮斯超市
ensembl_hsapiens <- useMart("ensembl",
dataset = "hsapiens_gene_ensembl")
代码的工作原理如下:
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)
all_homologues在使用它之后,我通过使用以下工具使其正常工作:
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))
所有_同系物
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)
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))