MISO (Mixture-of-Isoforms) is a probabilistic framework that quantitates the expression level of alternatively spliced genes from RNA-Seq data, and identifies diff erentially regulated isoforms or exons across samples. By modeling the generative process by which reads are produced from isoforms in RNA-Seq, the MISO model uses Ba yesian inference to compute the probability that a read originated from a particular isoform.

/usr/local/MISO

The following variables need to be set first:

csh or tcsh users:

% setenv LD_LIBRARY_PATH /usr/local/gsl-1.14/lib:$LD_LIBRARY_PATH
% setenv PATH /usr/local/Python/2.7.2/bin:$PATH

bash users:

$ export LD_LIBRARY_PATH=/usr/local/gsl-1.14/lib:$LD_LIBRARY_PATH
$ export PATH=/usr/local/Python/2.7.2/bin:$PATH

If you use the program frequently, it may be convenient to add the program directory to your path.

For a single session:

To add this directory to your path at login time automatically, add the appropriate line above to the end of your ~/.cshrc or ~/.bashrc file.

1. Create a script file. The file will contain the lines similar to the lines below. Modify the path of program location before running.

#!/bin/bash
# This file is YourOwnFileName
#
#PBS -N yourownfilename
#PBS -m be
#PBS -k oe

export PATH=/usr/local/Python/2.7.2/bin:$PATH
export LD_LIBRARY_PATH=/usr/local/gsl-1.14/lib:$LD_LIBRARY_PATH
cd /data/user/somewhereWithInputFile
python sam_to_bam.py --convert my_sample.sam my_sample/ --ref mm9_complete.fa.fai
python index_gff.py --index SE.mm9.gff indexed
....
....

2. Submit the script using the 'qsub' command on Biowulf, e.g. Note, user is recommend to run benchmarks to determine what kind of node is suitable for his/her jobs.

Useful commands:

freen: see http://biowulf.nih.gov/user_guide.html#freen

qstat: search for 'qstat' on http://biowulf.nih.gov/user_guide.html for it's usage.

jobload: search for 'jobload' on http://biowulf.nih.gov/user_guide.html for it's usage.

Using the 'swarm' utility, one can submit many jobs to the cluster to run concurrently.

Set up a swarm command file (eg /data/username/cmdfile). Here is a sample file:

There are one flag of swarm that's required '-f' and two other flags of swarm user most possibly needs to specify when submit a swarm job: '-t' and '-g'.

-f: the swarm command file name above (required)
-t: number of processors per node to use for each line of the commands in the swarm file above.(optional)
-g: GB of memory needed for each line of the commands in the swarm file above.(optional)

By default, each line of the commands above will be executed on '1' processor core of a node and uses 1GB of memory. If this is not what you want, you will need to specify '-t' and '-g' flags when you submit the job on biowulf.

Say if each line of the commands above also will need to use 10gb of memory instead of the default 1gb of memory, make sure swarm understands this by including '-g 10' flag:

[user@biowulf]$ swarm -g 10 -f cmdfile

For more information regarding running swarm, see swarm.html

User may need to run jobs interactively sometimes. Such jobs should not be run on the Biowulf login node. Instead allocate an interactive node as described below, and run the interactive job there.

User may add property of node in the qsub command to request specific interactive node. For example, if you need a node with 24gb of memory to run job interactivel y, do this:

http://genes.mit.edu/burgelab/miso/docs/