SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Bioinformatics (http://seqanswers.com/forums/forumdisplay.php?f=18)
-   -   SplicingCompass: error in processing CoverageBed files (http://seqanswers.com/forums/showthread.php?t=75579)

big_picture 04-21-2017 12:52 AM

SplicingCompass: error in processing CoverageBed files
 
1 Attachment(s)
Hey,

I want to do differential splicing analysis using SplicingCompass (http://www.leibniz-hki.de/files/cont...ssTutorial.pdf).

Necessary files:
- bam files (ideally from TopHat mapping)

- counted read numbers using CoverageBed, aligned to all exons in union transcript (script is provided by the authors, annotation file downloaded from UCSC database; hg38, annotation file: genes and gene predictions, NCBI_Refseq, GTF; id to symbol mapping file: genes and gene predictions, NCBI_Refseq, selected fields from primary and related tables --> name and name2)

- junction files from mapping.

I used mapping with STAR, so I converted the SJ_out files to junction.bed files using this script on github:
https://gist.github.com/fabiolib/ffb...0074aeffc54901
to convert them into the right format for SplicingCompass.

I had no problems until line 38, where I got the following error message:

Number of input files ok?
[1] TRUE
> countTable=new("CountTable")
> countTable=setExperimentInfo(countTable,expInf)
> countTable=constructCountTable(countTable,nCores=1,printDotPerGene=TRUE)
Processing CoverageBed output files:
|=============== | 12%
Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec, :
scan() expected 'a real', got 'D00614:183:CABNAANXX:3:1107:7896:66038'
>


It is unable to process the gff files, I couldn't identify where I made a mistake. Could you help me with it, please?

See the examples attached.

Thank you for the help!


All times are GMT -8. The time now is 11:18 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.