HyPhy message board | |
http://www.hyphy.org/cgi-bin/hyphy_forums/YaBB.pl
HYPHY Package >> HyPhy bugs >> Error: perBranchTree is undefined in call to ReplicateConstraint http://www.hyphy.org/cgi-bin/hyphy_forums/YaBB.pl?num=1378685638 Message started by stephaneg on Sep 8th, 2013 at 5:13pm |
Title: Error: perBranchTree is undefined in call to ReplicateConstraint Post by stephaneg on Sep 8th, 2013 at 5:13pm
Hi,
When using the MEME model on my data, I get the following error message: Error: perBranchTree is undefined in call to ReplicateConstraint. Function call stack 1 : Replicate the following constraint "this1.?.alpha:=_sFactor*this2.?.synRate__" using these variables:perBranchTree;codonTree; ------- 2 : result_String["BRANCHES"]=obtainBranchWiseEBEstimates(sFactor,nsFactor1,nsFactor2,mixingP,filterString) ------- 3 : Branch under condition 'SendJobMEME()' to <END EXECUTION> else <END EXECUTION> Am I doing something wrong? Regards, |
Title: Re: Error: perBranchTree is undefined in call to ReplicateConstraint Post by Sergei on Sep 9th, 2013 at 1:10pm
Hi there,
This looks like a bug in the MEME script: you are not using MPI, are you (e.g. running HYPHYMPI on a computer cluster)? Sergei |
Title: Re: Error: perBranchTree is undefined in call to ReplicateConstraint Post by stephaneg on Sep 10th, 2013 at 2:29pm
Dear Sergei,
Thanks for your reply. I'm not using MPI. Just MEME_mf.bf with one data partition (I realised aftwerwards that QuickSelectionDetection.bf, which works well, was more appropriate for the analysis I am interested in). On a slightly different topic, I understand that the p-values returned by MEME correspond to the null hypothesis of no positive selection at each codon site. How can I use this information to test for the null hypothesis of no positive selection anywhere in the alignment? Do you recommend applying a Bonferonni correction to the site-specific p-values and reject the null hypothesis of no positive selection anywhere in the alignment if at least one of them is below a nominal level? Or, shall I reject the null hypothesis if at least one non-corrected (i.e., an 'original' site-specific) p-value is below the nominal level? Regards, |
Title: Re: Error: perBranchTree is undefined in call to ReplicateConstraint Post by stephaneg on Sep 27th, 2013 at 12:01am Dear Sergei, Thanks again for your explanations. Just to make sure I understand properly: the type-I error frequencies reported in your Murrell et al. (2012) paper (SI page 2, Null simulations) correspond to the proportion of neutrally evolving sites (not alignments) at which positive selection was detected. Is this correct? Regards, -Stephane- |
Title: Re: Error: perBranchTree is undefined in call to ReplicateConstraint Post by JoseR on Mar 27th, 2014 at 11:08am
Hi,
I am using HYPHYMP to run MEME and I get a similar error message, see below. I used as input a file generated by Gard that works well on the web site. I have compiled HYPHY from the current tar ball, and was able to run Gard and FEL without problems. Any help will be appreciated. Regards, Jose Ribeiro Follows MEME output till crash: [MEME PHASE 0] Retuning branch lengths and nucleotide rates under the codon model... Improved Log(L) BY 14.81274319508702 points [MEME PHASE 1]. Fitting a codon model site-by-site | Codon | Beta- | Weight- | Beta+ | Weight+ | alpha | LRT | p-value |Full Log(L)| Est. remaining time (secs) Error: perBranchTree is undefined in call to ReplicateConstraint. Function call stack 1 : Replicate the following constraint "this1.?.alpha:=_sFactor*this2.?.synRate__" using these variables:perBranchTree;codonTree; ------- 2 : result_String["BRANCHES"]=obtainBranchWiseEBEstimates(sFactor,nsFactor1,nsFactor2,mixingP,filterString) ------- 3 : Branch under condition 'SendJobMEME()' to <END EXECUTION> else <END EXECUTION> |
Title: Re: Error: perBranchTree is undefined in call to ReplicateConstraint Post by JoseR on Mar 30th, 2014 at 9:26am
Hi Sergei,
I actually made it to work now! Just for the record, this is an example of the command line I am using, it may be of help to some other reader. (echo 1;echo 1;echo 1;echo 1;echo /home/hp/gard/AGAP028193-PA.gard_finalout;echo /home/hp/meme/AGAP028193-PA.gard_finalout-meme;echo 0.05;echo /home/hp/meme/AGAP028193-PA.gard_finalout-meme.csv;) | /usr/local/bin/mpirun -np 16 HYPHYMPI /usr/local/lib/hyphy/TemplateBatchFiles/MEME_mf.bf > /home/hp/meme/AGAP028193-PA.gard_finalout-MEME-out Many thanks for making all this available to the community, Jose Ribeiro |
HyPhy message board » Powered by YaBB 2.5.2! YaBB Forum Software © 2000-2024. All Rights Reserved. |