#这里是下载人类的已知基因的信息 (35.4 MB)
source("http://bioconductor.org/biocLite.R")
biocLite("TxDb.Hsapiens.UCSC.hg19.knownGene")
然后查看我们下载的这个包里面所包含的信息 Continue reading
#这里是下载人类的已知基因的信息 (35.4 MB)
source("http://bioconductor.org/biocLite.R")
biocLite("TxDb.Hsapiens.UCSC.hg19.knownGene")
然后查看我们下载的这个包里面所包含的信息 Continue reading
创建txdb对象的三种方法:
一是读取实例文件
library(GenomicFeatures)
samplefile <- system.file("extdata", "hg19_knownGene_sample.sqlite",
package="GenomicFeatures")
txdb <- loadDb(samplefile)
这个文件是hg19的一个子集,里面只有178个转录本,620个外显子 Continue reading