library(GenomicFeatures) metadata=data.frame(name="Resource URL",value=paste0("https://ftp.ensemblgenomes.ebi.ac.uk/pub/plants/release-56/gff3/triticum_aestivum/Triticum_aestivum.IWGSC.56.gff3.gz")) #从gtf文件中获取TxDb CS=makeTxDbFromGFF(file="~/genome/ChineseSpring/Triticum_aestivum.IWGSC.o...
value=paste0("https://ftp.ensemblgenomes.ebi.ac.uk/pub/plants/release-56/gff3/triticum_aestivum/Triticum_aestivum.IWGSC.56.gff3.gz"))#从gtf文件中获取TxDbCS=makeTxDbFromGFF(file="~/genome
使用gtf构建 library(GenomicFeatrues) spombe <- makeTxDbFromGFF("test.gff3") #从ucsc中下载,但是我没跑通 hg19.refseq.db <- makeTxDbFromUCSC(genome="hg19", table="refGene") image.png 方法2 从Bioconductor中直接下载,Bioconductor提供了30个TxDb包 if (!requireNamespace("BiocManager", quietly ...