SEQanswers

Go Back   SEQanswers > General



Similar Threads
Thread Thread Starter Forum Replies Last Post
ShoRAH pedrolance Bioinformatics 6 02-11-2013 08:09 AM
New ShoRAH installation problems... naragam Bioinformatics 0 01-14-2013 06:17 AM
PubMed: ShoRAH: estimating the genetic diversity of a mixed sample from next-generati Newsbot! Literature Watch 0 09-20-2011 03:00 AM

Reply
 
Thread Tools
Old 02-06-2013, 04:12 AM   #1
ckseq
Junior Member
 
Location: South Africa

Join Date: Feb 2013
Posts: 4
Question ShoRAH

Hi

I am new to the forum. I'm a postgrad student in South Africa, and I'm currently working on viral genomics. I'm trying to reconstruct quasispecies using ShoRAH from illumina short read data.

I'm looking for information about ShoRAH from users but can't seem to find anything anywhere.

We have an analysis running on short read data (illumina), it's been going for two days now. We just chose random numbers for the parameters and fear that my analysis might take weeks and weeks! Can anyone help me by telling me what the following parameters mean:

files
-i <input data file>
parameters
-j <sampling iterations>
-a <alpha>
-K <startvalue for number of clusters> not compat. with -k
-k <avg. number of reads in each startcluster> not compat. with -K
-t <history time>
-R <randomseed>
-----------------------------------------------------
-h this help!

The analysis is currently running on -j = 100, -t = 1000 and -K = 10. We left out the rest. The input data includes 90 000 reads mapped to a 11kb reference genome.

What is the signifficance/function of -a (alpha), -K (start value for number of clusters), -k (number of reads per start cluster), -t (history time) and -R (randomseed)? If I knew what these do it might be easier to try and optimise!!

Note: complete bioinformatics novice but trying
ckseq is offline   Reply With Quote
Reply

Tags
shorah

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 03:37 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO