微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

biomaRt-getSequence既不返回序列也不返回错误消息

如何解决biomaRt-getSequence既不返回序列也不返回错误消息

我有一套Sus scrofa的ENSEMBLE标识符,我想使用biomart获得它们各自的序列(cDNA,如果可能的话)。这对人类标识符没有任何问题,但是即使我可以使用ENSEMBLE网站手动检索我的标识符序列,我的猪数据也无法获得任何结果。

一个例子:

# Human - working
library(biomart)
ensembl <- useEnsembl(biomart = "genes",dataset = "hsapiens_gene_ensembl")
gene <- "ENSG00000175206"

biomartCacheClear()
seq <- getSequence(id      = gene,type    = "ensembl_gene_id",seqType = "cdna",mart    = ensembl)

# Pig - not working
library(biomart)
ensembl <- useEnsembl(biomart = "genes",dataset = "sscrofa_gene_ensembl")
gene <- "ENSSSCG00045003854"

biomartCacheClear()
seq <- getSequence(id      = gene,mart    = ensembl)

没有错误消息,seq保留为空:

> seq
[1] cdna            ensembl_gene_id
<0 Zeilen> (oder row.names mit Länge 0)

有人知道为什么它不起作用吗?也许我缺少了什么?由于存有序列信息,我相信不是因为不太流行的猪基因组数据。非常感谢您的帮助。

最好, 马丁

P.S:一个有关的小问题:对于上面的人类例子,我得到了3个具有相同标识符的不同cDNA序列。由于我不确定,它们之间到底有什么区别?

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。

相关推荐


Selenium Web驱动程序和Java。元素在(x,y)点处不可单击。其他元素将获得点击?
Python-如何使用点“。” 访问字典成员?
Java 字符串是不可变的。到底是什么意思?
Java中的“ final”关键字如何工作?(我仍然可以修改对象。)
“loop:”在Java代码中。这是什么,为什么要编译?
java.lang.ClassNotFoundException:sun.jdbc.odbc.JdbcOdbcDriver发生异常。为什么?
这是用Java进行XML解析的最佳库。
Java的PriorityQueue的内置迭代器不会以任何特定顺序遍历数据结构。为什么?
如何在Java中聆听按键时移动图像。
Java“Program to an interface”。这是什么意思?