Question: how to create index
gravatar for Bioinformatics
8 months ago by
Bioinformatics30 wrote:

how can i use Rsubread​ to make the index based on human reference? i could not find the command 

In their documentation, I can see the following 

ref <- system.file("extdata","reference.fa",package="Rsubread")

but I don't know how to build it based on human 


ADD COMMENTlink modified 8 months ago by Mike Smith3.0k • written 8 months ago by Bioinformatics30
gravatar for Mike Smith
8 months ago by
Mike Smith3.0k
EMBL Heidelberg / de.NBI
Mike Smith3.0k wrote:

In that example ref <- system.file("extdata","reference.fa",package="Rsubread") gives the path to a file.  You can see this simply by printing it out.  For me this is

> ref
[1] "/g/R-devel/library/Rsubread/extdata/reference.fa"

So you need to get a FASTA file for the genome you want to use, and then use the path to that file for ref, so something like ref <- "c:/path/to/hg19.fa" or ref <- "/path/to/hg19.fa" depending on if you're using Windows or not.

You should then be able to use the buildindex() function in the way described.

ADD COMMENTlink written 8 months ago by Mike Smith3.0k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 246 users visited in the last hour