Comments (7)
Hi, in order to run AbSplice for a certain tissue, you need to provide a SpliceMap of that tissue. We currently provide SpliceMaps for tissues available in GTEx (retina is not there).
However, creating a SpliceMap for a new tissue is straightforward, as long as you have RNA-seq data. To create SpliceMaps see this example notebook. As an input you need to provide a count table with split-read counts for each junction and sample (see Readme of splicemap package).
For the split-read counting you should use a tool that can detect de-novo splice sites. We compared FRASER, RegTools and STAR; which all led to similar results (see Supplementary Fig. 2a,c,d of the AbSplice publication). You do not need a huge cohort in order to create a SpliceMap and get robust estimates for Reference PSI values (see Supplementary Fig. 2e,f of the AbSplice publication).
Let me know if you need any further assistance on that.
from absplice.
Thanks for your reply! Let me do it. If I have other problem, I will post here. Thanks again for your help!
from absplice.
I generated splicemap_5 and splicemap_3 in juyputer notebook. I should throw them to some folders of absplice and run it?
from absplice.
yes. What genome version did you use? For hg19 and hg38 it is easy to do. In that case you just save your splicemaps as gzipped files and put them in the folder that is indicated in this config. It is the same folder that already contains the splicemaps that have been downloaded from the example usecase (/data/resources/downloaded_files/splicemap_{genome}/
, where genome is either hg19 or hg38). In case you did not run the example, just create this folder and put the splicemaps there (so create a folder called downloaded_files/splicemap_{genome}/
in here. You need to name your splicemaps in the same way as in the config, so in your case e.g Retina_splicemap_psi3_method=kn_event_filter=median_cutoff.csv.gz
, same for psi5.
Then just add Retina
to the list of tissues in the main config. From there just follow the steps from the Readme to run the workflow on your own vcf files.
from absplice.
I added retina 3 and 5 like you said, but it can't process
error here
Waiting at most 5 seconds for missing files.
MissingOutputException in rule download_splicemaps in file /home/cheneyu/absplice/example/workflow/./download/Snakefile, line 97:
Job 5 completed successfully, but some output files are missing. Missing files after 5 seconds. This might be due to filesystem latency. If that is the case, consider to increase the wait time with --latency-wait:
../data/resources/downloaded_files/splicemap_hg19/Retina_splicemap_psi3_method=kn_event_filter=median_cutoff.csv.gz
../data/resources/downloaded_files/splicemap_hg19/Retina_splicemap_psi5_method=kn_event_filter=median_cutoff.csv.gz
Shutting down, this might take some time.
Exiting because a job execution failed. Look above for error message
Complete log: .snakemake/log/2023-07-07T114724.268107.snakemake.log
from absplice.
How to forbid the download and get the right result?
from absplice.
I resolve those problems, thank you so much!
from absplice.
Related Issues (16)
- Add AbSplice to environment.yaml
- Convert EBM model to ONNX to get rid of interpret dependency HOT 3
- Can't get attribute 'EBMPreprocessor' on <module 'interpret.glassbox.ebm.ebm' HOT 2
- Add SpliceMap info to AbSplice output
- Bugs in Absplice implementation HOT 18
- Absplice-RNA HOT 3
- T2T CHM13v2.0 (hs1) support HOT 1
- can't install, stop at "Solving environment: |" HOT 2
- Container Environment problem HOT 2
- How to get CATs data HOT 7
- Splicemap_tissues HOT 3
- License HOT 1
- Can't find cause of the error HOT 2
- SpliceAI running time HOT 5
- Numpy Error appears in Example-Workflow-Execution HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from absplice.