G-PhoCS on XSEDE
1.3.2
A Generalized Phylogenetic Coalescent Sampler
Gronau I, Hubisz MJ, Gulko B, Danko CG, Siepel A.
Gronau I, Hubisz MJ, Gulko B, Danko CG, Siepel A. Bayesian inference of ancient human demography from individual genome sequences. Nature Genetics 43 1031–1034. 2011
Phylogeny / Alignment
gphocs_xsede_132
infile
Input sequence file
sequences.txt
gphocs_invocation
perl
"convertFile.sh ctrlfile.ctrl;g-phocs_1.3.2_expanse ctrlfile.ctrl"
1
gphocs_invocation3
perl
" -n $numthreads -v"
80
gphocs_scheduler
scheduler.conf
perl
"threads_per_process=$numthreads\\n" .
"mem=" . (int(($numthreads)*(248/128))) . "G\\n" .
"node_exclusive=0\\n" .
"nodes=1\\n"
0
all_output
All output
*
runtime
1
scheduler.conf
Maximum Hours to Run (click here for help setting this correctly)
perl
"runhours=$value\\n"
0.25
Maximum Hours to Run must be less than 168
perl
$runtime > 168.0
Maximum Hours to Run must be greater than 0.1
perl
$runtime < 0.1
The job will run on $numthreads processors as configured. If it runs for the entire configured time, it will consume $numthreads x $runtime cpu hours
perl
$runtime ne 0
Estimate the maximum time your job will need to run. We recommend testimg initially with a < 0.5hr test run because Jobs set for 0.5 h or less depedendably run immediately in the "debug" queue.
Once you are sure the configuration is correct, you then increase the time. The reason is that jobs > 0.5 h are submitted to the "normal" queue, where jobs configured for 1 or a few hours times may
run sooner than jobs configured for the full 168 hours.
primary_ctrlfile
Select CTRL file (it must name the input file sequences.txt)
ctrlfile.ctrl
Plese make sure your ctrl file specifies that the seq-file name is sequences.txt
perl
defined $primary_ctrlfile
secondary_ctrlfile
Select secondary CTRL file (optional)
perl
defined $value ? "ctrlfile2.ctrl" : ""
3
ctrlfile2.ctrl
numthreads
How many threads?
Sorry, threads must be no less than 6
perl
$numthreads < 6
Sorry, threads must be no more than 24
perl
$numthreads > 24
Please enter a thread value
perl
!defined $numthreads