search

timbitz / Whippet.jl

Lightweight and Fast; RNA-seq quantification at the event-level
MIT License
105 stars 21 forks source link

Error with whippet-delta.jl #80

Closed paulinefx closed 5 years ago

paulinefx commented 5 years ago

Hello,

It is the first time, I am using whippet and julia, so I am discovering a bit and I don't understand this error. For my two samples I created their customed index with their bam file using whippet-index.jl than I have done the quantification withe whippet-quant.jl. Have you ever seen this error ?

Thank you for your help, Pauline

Whippet v0.11 loading and compiling... 30.385936 seconds Sample A: 1441_T_output.psi.gz Sample B: 1441_S_output.psi.gz Now processing files and calculating posterior distributions... ERROR: LoadError: AssertionError: Incorrect events matched!! Stacktrace: [1] #process_psi_files#100(::Int64, ::Int64, ::Float64, ::Int64, ::Function, ::String, ::Array{BufferedStreams.BufferedInputStream,1}, ::Array{BufferedStreams.BufferedInputStream,1}) at /data/home/hfx405/.julia/v0.6/Whippet/src/diff.jl:125 [2] (::Whippet.#kw##process_psi_files)(::Array{Any,1}, ::Whippet.#process_psi_files, ::String, ::Array{BufferedStreams.BufferedInputStream,1}, ::Array{BufferedStreams.BufferedInputStream,1}) at ./:0 [3] macro expansion at /data/home/hfx405/.julia/v0.6/Whippet/src/timer.jl:5 [inlined] [4] main() at /data/home/hfx405/.julia/v0.6/Whippet/bin/whippet-delta.jl:91 [5] include_from_node1(::String) at ./loading.jl:576 [6] include(::String) at ./sysimg.jl:14 [7] process_options(::Base.JLOptions) at ./client.jl:305 [8] _start() at ./client.jl:371 while loading /data/home/hfx405/.julia/v0.6/Whippet/bin/whippet-delta.jl, in expression starting on line 5

timbitz commented 5 years ago

Hi @paulinefx, yes I have seen this assertion error before-- I created it! :) You would get this error if you are using two different indexes for alignment... You can only compare .psi.gz files that were created from the same index. Hope that helps.

paulinefx commented 5 years ago

Yes that helps a lot. I went too fast in the creation of the index and forgot to merge my bam files. Thank you