User Tools

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision
Previous revision
anthill_sri011 [2019/07/06 13:00] tmatejukanthill_sri011 [2023/08/01 01:08] (current) – external edit 127.0.0.1
Line 6: Line 6:
  
 ==== prepare software ==== ==== prepare software ====
-Login to anthill23 and download rMATS along with some test dataset, and gft files:+Login to anthill23 and download rMATS along with some test datasets, and gft files((downloading this files will take about 3 minutes)):
  
   @anthill23:~$   @anthill23:~$
Line 20: Line 20:
   wget ftp://ftp.ensembl.org/pub/grch37/current/gtf/homo_sapiens/Homo_sapiens.GRCh37.87.gtf.gz   wget ftp://ftp.ensembl.org/pub/grch37/current/gtf/homo_sapiens/Homo_sapiens.GRCh37.87.gtf.gz
   gunzip Homo_sapiens.GRCh37.87.gtf.gz   gunzip Homo_sapiens.GRCh37.87.gtf.gz
-  wget ftp://ftp.ensembl.org/pub/release-91/gtf/homo_sapiens/Homo_sapiens.GRCh38.91.gtf.gz 
-  gunzip Homo_sapiens.GRCh38.91.gtf.gz 
  
 ==== prepare gsl lib ==== ==== prepare gsl lib ====
-rMATS uses gsl library that need to be prepared separately((takes about five minutes)). This library have to be compiled on computing node, so it is done in interactive mode ( +rMATS uses gsl library that need to be prepared separately((takes about five minutes)). This library have to be compiled on computing node, so it is done in interactive mode ([[howtouseslurm_000|SLURM, batch and interactive mode]]) :
-[[howtouseslurm_000|SLURM, batch and interactive mode]]) :+
  
   @anthill23:~$   @anthill23:~$
      
-  srun -N1 -n4 --pty bash -l+  srun -J gsl_compile -N1 -n4 --pty bash -l
      
   cd ~/anthill23_rmats/   cd ~/anthill23_rmats/
Line 41: Line 38:
   cd ~/anthill23_rmats/gsl-2.4/lib   cd ~/anthill23_rmats/gsl-2.4/lib
   ln -s libgsl.so libgsl.so.0   ln -s libgsl.so libgsl.so.0
-   +   
 +  exit
  
 +==== simple test ====
 +If above steps ( prepare software, and prepare gsl lib ) finish successfully this simple test should finish with success. Do this test in  interactive mode ([[howtouseslurm_000|SLURM, batch and interactive mode]]) :
  
 +  @anthill23:~$
 +  
 +  srun -J rMATS_test --pty bash -l
 +  
 +  cd /home/users/${USER}/anthill23_rmats
 +  export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:/home/users/${USER}/anthill23_rmats/gsl-2.4/lib
 +  python rMATS.4.0.2/rMATS-turbo-Linux-UCS4/rmats.py --b1 testData/b1.txt --b2 testData/b2.txt --gtf gtf/Homo_sapiens.GRCh37.87.gtf --od bam_test -t paired --readLength 50 --cstat 0.0001 --libType fr-unstranded
 +  
 +  ...
 +  #read output and exit interactive mode when finished
 +  exit
  
 ==== sbatch example ==== ==== sbatch example ====
-...+Bellow example was taken form [[http://rnaseq-mats.sourceforge.net/user_guide.htm|rMATS User Guide]]. This description assumes that above steps ( prepare software, and prepare gsl lib ) were done beforehand if not, remember to modify them accordingly.  
 + 
 +File ''rmats_test.batch'' 
  
 <code> <code>
-... +#!/bin/bash -l 
-</code>+#SBATCH --partition=test 
 +#SBATCH --ntasks=4 
 +#SBATCH --mem 4G 
 +#SBATCH --job-name rMATS_test
  
 +cd /home/users/${USER}/anthill23_rmats
 +export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:/home/users/${USER}/anthill23_rmats/gsl-2.4/lib
  
-==== performance tests ==== +python rMATS.4.0.2/rMATS-turbo-Linux-UCS4/rmats.py \ 
-...+ --b1 b1.txt \ 
 + --b2 b2.txt \ 
 + --gtf gtf/Homo_sapiens.GRCh37.87.gtf \ 
 + --od bam_test \ 
 + -t paired \ 
 + --readLength 50 \ 
 + --cstat 0.0001 \ 
 + --libType fr-unstranded
  
-<code> 
-... 
 </code> </code>
 +
 +==== performance tests ====
 +None. At this time.
 +
  
  
  

This website uses cookies. By using the website, you agree with storing cookies on your computer. Also, you acknowledge that you have read and understand our Privacy Policy. If you do not agree, please leave the website.

More information