Welcome, Guest. Please Login
YaBB - Yet another Bulletin Board
 
  HomeHelpSearchLogin  
 
Tips vs Internals in SLAC (Read 1845 times)
noob
YaBB Newbies
*
Offline


Feed your monkey!

Posts: 8
Tips vs Internals in SLAC
Mar 4th, 2010 at 5:43am
 
Hi,

When I run SLAC within QuickSelectionDetection.bf and select "Tips vs Internals" I get the following error:

Invalid Matrix Index [1830][56] in a 61 by 61 matrix.
Call stack
2 : _SITE_OS_COUNTInternals[cd1]
[cd2]=_SITE_OS_COUNTInternals[cd1][cd2]+1
1: ExecuteComands in string "#include \"SGlvL.bf\";"
reading input from Macintosh
HD:Research:Software:HyPhySecond

When I run it over the full tree it seems to work fine. What am I doing wrong?

Also, could you recommend some reasonably appropriate analyses to detect sites under selection on an alignment of drug naive and post drug HIV RT sequences?

Thanks in advance
Back to top
 
 
IP Logged
 
Sergei
YaBB Administrator
*****
Offline


Datamonkeys are forever...

Posts: 1658
UCSD
Gender: male
Re: Tips vs Internals in SLAC
Reply #1 - Mar 4th, 2010 at 7:27am
 
Hi there,

Sorry about the error message; it is a bug in our code that crops up with gappy sites and I will fix it shortly.
I would recommend running your sequences through all the methods (SLAC/FEL/REL) in datamonkey.org and seeing what you find. FEL is probably the best single method for this purpose, if you had to choose one.

Sergei
Back to top
 

Associate Professor
Division of Infectious Diseases
Division of Biomedical Informatics
School of Medicine
University of California San Diego
WWW WWW  
IP Logged
 
noob
YaBB Newbies
*
Offline


Feed your monkey!

Posts: 8
Re: Tips vs Internals in SLAC
Reply #2 - Mar 4th, 2010 at 12:01pm
 
Hi Sergei,

Thanks for the advice! Let me know when you have fixed that bug. I also get some matrix error when running DEPS, but my computer is currently tied up running something else, so I won't be able to replicate it for a while. Is it perhaps something similar?

Regards,
noob
Back to top
 
 
IP Logged