Dear SUPPA developers,
I’m encountering an error when trying to compute PSI values using the psiPerEvent subcommand. I suspect it’s related to the format of either my .ioe or expression file, but I cannot resolve it.
What I am trying to do
I’m analyzing alternative splicing differences between freshwater and seawater samples using the following workflow:
Generate AS events from GTF:
python /home/software/SUPPA2/SUPPA-2.4/suppa.py generateEvents \
-i ../1.6.merge.agat.gtf \
-o merge.events \
--event-type SE SS MX RI FL \
-f ioe \
-b S
Merge all .ioe files:
awk 'FNR==1 && NR!=1 { next } 1' *.ioe > merge.events.all.ioe
Run PSI estimation:
python /home/software/SUPPA2/SUPPA-2.4/suppa.py psiPerEvent \
-i /bio/sup_drs/6.2.AS/1_step1/merge.events.all.ioe \
-e suppa_tpm.tsv \
-o general_events
Error Message
When running psiPerEvent, I get the following error repeatedly:
ERROR:lib.tools:<line_number>. Unexpected number of fields. 6 expected, 5 given. Skipping line...
Input Files
merge.events.all.ioe (first 3 lines):
seqname gene_id event_id alternative_transcripts total_transcripts
ptg000001l DRS_000562 DRS_000562;A3:ptg000001l:29280-30515:29220-30515:- Lm_00001430 Lm_00001431,Lm_00001430
ptg000001l DRS_000562 DRS_000562;A3:ptg000001l:29280-30515:29118-30515:- Lm_00001430 Lm_00001432,Lm_00001430
suppa_tpm.tsv (first 5 lines):
sample1 sample2 sample3 sample4
Lm_00000001 7.16778 8.68441 7.65165 8.27315
Lm_00000002 0.40268 0.81928 0.40271 0.30084
Lm_00000003 2.33556 2.04821 2.65794 1.27857
Lm_00000004 24.8053 25.3159 24.1631 21.8110
All fields are tab-delimited.
Question
What could be causing this "6 expected, 5 given" error? Is there something wrong with my .ioe format or TPM file?
Thank you for your help!
Dear SUPPA developers,
I’m encountering an error when trying to compute PSI values using the psiPerEvent subcommand. I suspect it’s related to the format of either my .ioe or expression file, but I cannot resolve it.
What I am trying to do
I’m analyzing alternative splicing differences between freshwater and seawater samples using the following workflow:
Generate AS events from GTF:
Merge all .ioe files:
Run PSI estimation:
Error Message
When running psiPerEvent, I get the following error repeatedly:
Input Files
merge.events.all.ioe (first 3 lines):
suppa_tpm.tsv (first 5 lines):
All fields are tab-delimited.
Question
What could be causing this "6 expected, 5 given" error? Is there something wrong with my .ioe format or TPM file?
Thank you for your help!