【问题标题】:Excel: Ensembl continue, how to save to CSV file after change RExcel:Ensembl 继续,更改 R 后如何保存到 CSV 文件
【发布时间】: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 列旁边创建一个列,两者都很好。

标签: r biomart


【解决方案1】:

你正在取回 ensembl 中的所有条目,你需要合并:

my_ids = data.frame(ids=c("ENSG00000130203","ENSG00000136717","ENSG00000186868"))

df = merge(my_ids,results_end_1,by.x="ids",by.y="ensembl_gene_id")
df
              ids external_gene_name
1 ENSG00000130203               APOE
2 ENSG00000136717               BIN1
3 ENSG00000186868               MAPT

或者您需要设置“过滤器”选项以获取 my_ids 中的条目,:

out <- getBM(attributes = c("ensembl_gene_id","external_gene_name"), 
values = my_ids, mart = ensembl,filters="ensembl_gene_id")

  ensembl_gene_id external_gene_name
1 ENSG00000130203               APOE
2 ENSG00000136717               BIN1
3 ENSG00000186868               MAPT

你可以写:

write.csv(out,"out.csv",row.names=FALSE,quote=FALSE)

【讨论】:

    猜你喜欢
    • 1970-01-01
    • 1970-01-01
    • 1970-01-01
    • 2016-03-20
    • 2012-12-31
    • 1970-01-01
    • 2017-03-17
    • 2011-04-03
    • 1970-01-01
    相关资源
    最近更新 更多