作者热门文章
- html - 出于某种原因,IE8 对我的 Sass 文件中继承的 html5 CSS 不友好?
- JMeter 在响应断言中使用 span 标签的问题
- html - 在 :hover and :active? 上具有不同效果的 CSS 动画
- html - 相对于居中的 html 内容固定的 CSS 重复背景?
假设我有一个名为“Tlr6”的基因(见下图),我想知道如何在 R 中检索染色体上该基因的起始值和终止值?例如图中,起始值为64952031,结束值为64960097。
图片网址为here 。这里的基因名称是 Tlr6
,Ensembl ID 是 ENSMUSG00000051498
。我可以使用这些信息将开始值和结束值读入 R 吗?
最佳答案
您可以使用 Bioconductor
中的 biomaRt
包来执行此操作:
#skip this if the package is already installed
source("http://bioconductor.org/biocLite.R")
biocLite("biomaRt")
library(biomaRt)
#select the ensembl mouse dataset
ensembl <- useMart("ensembl", dataset="mmusculus_gene_ensembl")
getBM(attributes=c("ensembl_gene_id","start_position","end_position"),
filters="ensembl_gene_id",values = "ENSMUSG00000051498",mart=ensembl)
有关该包的更多信息 here 。
关于r - 如何从 Ensembl 检索染色体上基因的起始值和终止值,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/31680148/
我已经开始尝试使用 Jenetics 库,但是我在尝试制作一组非常简单的“自定义”基因/染色体时遇到了一些问题。我试图做的是创建一个自定义染色体,其中包含不同(随机)数量的自定义基因。为了简单起见
我是一名优秀的程序员,十分优秀!