Rsubreads - featureCount - gtf - error
3
1
Entering edit mode
@samuel-collombet-6574
Last seen 6.9 years ago
France

Hi,

I would like to use Rsubreads to coutn feature on exon (summarizing at the gene ID level), however featureCounts does not accept my gtf:

...
> || Load annotation file ~/genomes/mus_musculus/mm10/annotation/ensembl_mu ... ||
> Failed to open the annotation file ~/genomes/mus_musculus/mm10/annotation/ensembl_mus_musculus.GRCm38.75.gtf, > or its format is incorrect, or it contains no 'exon' features.
> Error in file(file, "rt") : cannot open the connection
> In addition: Warning message:
> In file(file, "rt") :
>   cannot open file './.Rsubread_featureCounts_pid8772': No such file or directory

What surprise me is that when I run it in bash using featureCounts, same annotation file, I have no problem and got my resut.

 

Am I doing somehting wrong?

Samuel

 

rsubread featurecounts • 4.0k views
ADD COMMENT
2
Entering edit mode
Wei Shi ★ 3.6k
@wei-shi-2183
Last seen 2 days ago
Australia/Melbourne/Olivia Newton-John …

Providing the full path to your GTF file should resolve this problem.

Wei

ADD COMMENT
0
Entering edit mode
Johannes Rainer ★ 2.0k
@johannes-rainer-6987
Last seen 7 weeks ago
Italy

I am not sure, but this error message could also mean that you have no write permission in your current directory where you started R (what I mean is the working directory).

Alternatively, what about providing the full path instead of the ~/?

cheers, jo
 

ADD COMMENT
0
Entering edit mode
@samuel-collombet-6574
Last seen 6.9 years ago
France

Thanks Wei and Johanness, full path of gtf was the thing! Should have thought about it...

(sorry for not answering before)

ADD COMMENT

Login before adding your answer.

Traffic: 776 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