【发布时间】:2020-02-19 16:32:08
【问题描述】:
相关 Unable to use biomaRt package to get Gene Symbols from Entrez IDs
代码运行起来很吸引人,但我想将它保存到我的 当前 CSV 文件中,而不仅仅是“查看”它(添加一个额外的列也可以完成这项工作,我获得超过 65k 个 ID)
library(biomaRt)
marts <- listMarts()
ensembl <- useMart("ensembl")
datasets <- listDatasets(ensembl)
ensembl=useDataset("hsapiens_gene_ensembl",mart=ensembl)
attributes <- listAttributes(ensembl)
my_ids <- read.csv("LUADhtseq.csv")
ensembl = useEnsembl(biomart="ensembl", dataset="hsapiens_gene_ensembl")
results_end_1 <- getBM(attributes = c("ensembl_gene_id","external_gene_name"), values = my_ids, mart = ensembl )
View(results_end_1)
【问题讨论】:
-
所以您想在 my_ids data.frame 中添加另一列?
-
write.csv(results_end_1, "filename.csv") -
对不起,我不够清楚。我正在尝试更改我已经在使用“LUADhtseq”的文件中的 Ensembl ID,或者在 Ensembl_ID 列旁边创建一个列,两者都很好。