HyPhy message board | |
http://www.hyphy.org/cgi-bin/hyphy_forums/YaBB.pl
HYPHY Package >> HyPhy bugs >> problems with GARDProcessor.bf http://www.hyphy.org/cgi-bin/hyphy_forums/YaBB.pl?num=1221161656 Message started by luis on Sep 11th, 2008 at 12:34pm |
Title: problems with GARDProcessor.bf Post by luis on Sep 11th, 2008 at 12:34pm
dear sergei
I again, then run GARDProcessor.bf show error: (echo data/demo.seq; echo /home/jlacosta/Download/programas/HYPHY_Source/GARD_demo.out;) | mpirun -np 2 HYPHY TemplateBatchFiles/GARDProcessor.bf /home/jlacosta/Download/programas/HYPHY_Source/TemplateBatchFiles/Please load a nucleotide data file:: /home/jlacosta/Download/programas/HYPHY_Source/TemplateBatchFiles/Please load a GA partition analysis result file:: Loaded 1023 partitions Sequences :4 Sites :273 Variable :44 f(A) = 0.458791 f(C) = 0.141941 f(G) = 0.159341 f(T) = 0.239927 c-AIC = 1180.51 Log Likelihood = -576.553622490902; Shared Parameters: GT=0.0702169=0.0702169 CT=0.637412=0.637412 betaP=0.856559=0.856559 CG=0.133953=0.133953 AT=0.201307=0.201307 AC=0.388101=0.388101 alpha=0.0234925=0.0234925 betaQ=2.97442=2.97442 Tree givenTree=(a:0.0326137,b:0.00689777,(c:0.0266298,og:0.0874902)Node3:0.0488684); Mean divergence : 10.9395% Fitting a single-tree, multiple partition model mpirun noticed that job rank 0 with PID 5382 on node FIERA exited on signal 11 (Segmentation fault). 1 additional process aborted (not shown) ?????? sorry!!! |
Title: Re: problems with GARDProcessor.bf Post by Sergei on Sep 11th, 2008 at 6:17pm
Dear Luis,
I think you may be using an old vertion of GARDProcessor.bf The line [code] Loaded 1023 partitions [/code] tells me that the code is not reading the number of partitions correctly. Try using GARDProcessor.bf from the Standard Analyses menu (under Recombination). I attach a folder containing a simple example. Cheers, Sergei http://www.hyphy.org/cgi-bin/hyphy_forums/YaBB.pl?action=downloadfile;file=GARD_test.zip (4 KB | )
|
Title: Re: problems with GARDProcessor.bf Post by luis on Sep 11th, 2008 at 7:10pm
dear sergei
whole procedure was:: $ ( echo ; echo data/demo.seq; echo ; echo 010010; echo 1; echo \n;) | mpirun -np 8 HYPHYMPI TemplateBatchFiles/GARD.bf show it: GA has considered 355/ 13244 (1230 over all runs) unique models Total run time 0 hrs 1 mins 25 seconds Throughput 14.47 models/second Allocated time remaining 999 hrs 58 mins 35 seconds (approx. 5.20929e+07 more models.) GENERATION 50 with 3 breakpoints (~96% converged) Breakpoints c-AIC Delta c-AIC [BP 1] [BP 2] [BP 3] 0 1177.93 1 1170.30 7.629 114 2 1167.12 3.174 114 206 3 1167.12 0.000 114 206 GA has considered 362/ 13244 (1237 over all runs) unique models Total run time 0 hrs 1 mins 25 seconds Throughput 14.55 models/second Allocated time remaining 999 hrs 58 mins 35 seconds (approx. 5.23894e+07 more models.) GENERATION 51 with 3 breakpoints (~98% converged) Breakpoints c-AIC Delta c-AIC [BP 1] [BP 2] [BP 3] 0 1177.93 1 1170.30 7.629 114 2 1167.12 3.174 114 206 3 1167.12 0.000 114 206 GA has considered 369/ 13244 (1244 over all runs) unique models Total run time 0 hrs 1 mins 26 seconds Throughput 14.47 models/second Allocated time remaining 999 hrs 58 mins 34 seconds (approx. 5.20732e+07 more models.) GENERATION 52 with 3 breakpoints (~100% converged) Breakpoints c-AIC Delta c-AIC [BP 1] [BP 2] [BP 3] 0 1177.93 1 1170.30 7.629 114 2 1167.12 3.174 114 206 3 1167.12 0.000 114 206 GA has considered 376/ 13244 (1251 over all runs) unique models Total run time 0 hrs 1 mins 26 seconds Throughput 14.55 models/second Allocated time remaining 999 hrs 58 mins 34 seconds (approx. 5.23662e+07 more models.) Performing the final optimization... I generate these files: messages.log.mpinode0 messages.log.mpinode1 messages.log.mpinode2 messages.log.mpinode3 messages.log.mpinode4 messages.log.mpinode5 messages.log.mpinode6 messages.log.mpinode7 where is gard.splits ??? I take the last file that generates GARD.bf (echo data/demo.seq; echo /home/jlacosta/Download/programas/HYPHY_Source/messages.log.mpinode7;) | ./HYPHYMP TemplateBatchFiles/RecombinationProcessor.bf show me: TemplateBatchFiles/RecombinationProcessor.bf /home/jlacosta/Download/programas/HYPHY_Source/TemplateBatchFiles/Please load a nucleotide data file:: /home/jlacosta/Download/programas/HYPHY_Source/TemplateBatchFiles/Please load a GA partition analysis result file:: Loaded 2970.5 partitions Sequences :4 Sites :273 Variable :44 f(A) = 0.458791 f(C) = 0.141941 f(G) = 0.159341 f(T) = 0.239927 c-AIC = 1180.51 Log Likelihood = -576.553622490902; Shared Parameters: GT=0.0702169=0.0702169 CT=0.637412=0.637412 betaP=0.856559=0.856559 CG=0.133953=0.133953 AT=0.201307=0.201307 AC=0.388101=0.388101 alpha=0.0234925=0.0234925 betaQ=2.97442=2.97442 Tree givenTree=(a:0.0326137,b:0.00689777,(c:0.0266298,og:0.0874902)Node3:0.0488684); Mean divergence : 10.9395% Fitting a single-tree, multiple partition model my question is how to generate file GARD.split??? Sorry for the inconvenience and thank you in advance. |
Title: Re: problems with GARDProcessor.bf Post by Sergei on Sep 11th, 2008 at 8:40pm
Dear Luis,
GARD.bf prompts you (the last echo statement in your example) for the location of the results file. That file will be written by GARD along with several other files in the same path. On of them will end with _splits - that's the one you want to feed to the processor. E.g. if you request that the results be written to /home/user/results/main then the splits file will be in /home/user/results/main_splits Simply edit your echo statements to include a non-trivial path and all should go well. Cheers, Sergei PS message.log files only contain diagnostic messages about the run |
Title: Re: problems with GARDProcessor.bf Post by luis on Sep 11th, 2008 at 8:56pm
dear sergei!!
everything is perfect! thank you very much!!! |
HyPhy message board » Powered by YaBB 2.5.2! YaBB Forum Software © 2000-2024. All Rights Reserved. |