-
Notifications
You must be signed in to change notification settings - Fork 8
Description
Hi
I ran two chains in PhyloBayes-MPI with the following commands:
mpirun -n 24 pb_mpi -d alig.phy -catfix C20 -lg -t intree.tree chain_0
mpirun -n 24 pb_mpi -d alig.phy -catfix C20 -lg -t intree.tree chain_1
I stopped the chains after 10000 generations.
tail chain_0.trace
10650 353.805 93 -1195936.045 0.2836247507 0.146153911 20 2.427761501 20
10651 339.699 92 -1195329.104 0.2837520955 0.1446938163 20 2.428734088 20
10652 353.853 92 -1195153.129 0.2843045745 0.1452920026 20 2.427999614 20
tail chain_1.trace
10668 362.94 93 -1195727.005 0.2851022857 0.1412945893 20 2.428278018 20
10669 358.809 93 -1195032.63 0.285462945 0.1458959495 20 2.427681612 20
10670 359.08 93 -1195585.968 0.2865275251 0.1460400457 20 2.427547846 20
I stopped the runs softly with:
echo 0 > chain_0.run
echo 0 > chain_1.run
Then, I ran:
tracecomp chain_0 chain_1
setting upper limit to : 10539
chain_0.trace burnin : 2107 sample size : 8432
chain_1.trace burnin : 2107 sample size : 8432
name effsize rel_diff
loglik 3 0.00649395
length 3 1.67442
alpha 3 1.81706
Nmode 4217 1.90942
statent 3 1.87607
statalpha 3 1.89889
And I obtained those bad numbers
Then I ran:
bpcomp chain_0 chain_1
chain_0.treelist error
(
and I obtained an error with the tree file.
What would be happening??
Could you help me please?
Thanks,
lisy