Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Splice analysis using external data cannot be merged with local files #452

Open
zshu-alt opened this issue Apr 12, 2023 · 2 comments
Open

Comments

@zshu-alt
Copy link

Hello!
I was able to successfully run the FRASER analysis with 11 local bam files. To increase the sample size, I added the data in GTEx. The operation went smoothly at first. But in running rule AberrantSplicing_pipeline_Counting_03_filter_expression_FraseR_R:
There was an error like this:
Loading provided [counts]
Wed Apr 12 09:57:14 2023: Merging data ...
Wed Apr 12 09:57:17 2023: Calculate the PSI 5 and 3 values ...
Wed Apr 12 09:57:26 2023: Calculate the PSI site values ...
Error: subscript is a logical vector with out-of-bounds TRUE values
In addition: Warning message:
In .merge_two_Seqinfo_objects(x, y) :
The 2 combined objects have no sequence levels in common. (Use
suppressWarnings() to suppress this warning.)
Execution halted
Is it a merger problem? My bam file is compared by STAR software and fa in STAR database. I am sure it is hg19.
I uploaded my relevant files
How can I solve this problem? Increasing the sample size by external counting is important to me
I would really appreciate it if you could answer my questions.

sample_annotation_put.txt
nohup.log
config_put.txt

@vyepez88
Copy link
Collaborator

Hi! Weird, DROP should be able to recognize different chr styles. Can you double check if indeed your local BAM files have a different chr style than the external junction count matrices?

@AtaJadidAhari
Copy link
Collaborator

Hi @zshu-alt, were you able to sort this out?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants