This is a test version of Biostars. For the public version, visit https://www.biostars.org.
featureCounts error???
# get gtf
# wget #https://ftp.ensembl.org/pub/release-110/gtf/mus_musculus/Mus_musculus.GRCm39.110.gtf.gz

input.dir <- "/home/laudy/data/featurecounts/"
setwd(input.dir)

featureCounts -p -O -T -a /input.dir/Mus_musculus.GRCm39.110.gtf -o /input.dir/quants.txt /input.dir/PMN_CTRAligned.sortedByCoord.RD.RG.RC.out.bam

please can someone tell me what's wrong? I tried all the options and he give me the same error:

Error: object 'p' not found
or
Error: unexpected symbol in "featureCounts

Thanks

featurecounts
• 1,241 views • link
updated by Ram Ram • written by Samer • Samer

-T specifies the number of threads, so should have an integer after it.

It's also important to make sure you have the latest version. There were some changes in how the program handles paired-end reads, so to be sure you are doing it correctly check whether you version is current and add the --countReadPairs argument.

• 0 views • link
written by rpolicastro rpolicastro

Subread is implemented both as an Unix program or as an R package. OP seems to be mixing the two versions, and I suspect they have typed a Unix command at the R prompt instead of at the Unix prompt. Typing featureCounts -p etc at the R prompt with Rsubread installed will lead to the object p not found error message.

• 0 views • link
written by Gordon Smyth Gordon Smyth

0 answers

No answers yet.

Log in to answer this question.

More posts like this