how to create index
1
0
Entering edit mode
@bioinformatics-10931
Last seen 2.7 years ago
United States

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")
buildindex(basename="reference_index",reference=ref)​

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

 

Rsubread • 1.1k views
ADD COMMENT
0
Entering edit mode
Mike Smith ★ 6.5k
@mike-smith
Last seen 8 hours ago
EMBL Heidelberg

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 COMMENT

Login before adding your answer.

Traffic: 621 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6