Hello,

I'm trying to implement a PAML M8a model via Hyphy using the NielsenYang.bf batch file. However, I keep getting the message "Abort Trap: 6" when running. The errors.log and messages.log remain empty, as well.

I have tried to use both the direct command line interface as well as an inputRedirect, and both yield the error. I'm attaching my file for inputRedirect here. I've tried it with a variety of alignments/trees, all of which are fasta and newick, respectively, and the abort trap keeps occurring.


Any thoughts would be much appreciated.

Thanks,
Stephanie Spielman

Hi Stephanie,

I haven't seen this behavior on OS X before: could you provide more details for me?

1). Which version of HyPhy are you running (command line, GUI, version number)?
2). What is the OS X version?
3). How much of the analysis finishes before you see the Abort Trap raised?

Sergei

Hi Stephanie,

Your HyPhy version seems a bit out of date -- please use the current 2.2 version (Multimedia File Viewing and Clickable Links are available for Registered Members only!!  You need to Login). There have been occasional crashes on OS 10.7 and 10.8 due to some multithreading issues. Try including the CPU=1 argument on the command line: this will force HyPhy to run in single threaded mode, and if the issue no longer occurs, then the problem is due to multithreading (this issue is not limited to HyPhy, btw, but many packages that use OpenMP).

Beta & omega should indeed have one category greater than 1 (which can actually lead to bad statistical performance, see Multimedia File Viewing and Clickable Links are available for Registered Members only!!  You need to Login. If you see more than one omega > 1 (presumably the others are close to 1?), then this is a numerical issue with HyPhy, but also indicative of the fact that beta & omega is a bad fit to your data (i.e. the beta distribution gets 'bunched' on the boundary). Are the p and q parameters high or low?

Unless you have compelling reasons to use Beta & omega, I would suggest you use FUBAR instead: it's much faster and more reliable: Multimedia File Viewing and Clickable Links are available for Registered Members only!!  You need to Login

Sergei

Hi Stephanie,

I agree with your sentiment on PAML's models implementations; I think what you see is a numerical artifact due to a 'weird' beta distribution shape. Can you send me the model fit that HyPhy writes out as a part of the analysis, so that i trouble-shoot the problem?

Sergei