MrBayes
MrBayes copied to clipboard
Error in sumt vstat file
The sumt vstat file erroneously outputs the PSRF values for each of the runs instead of presenting them for the ensemble sample just once. This causes a mismatch with the header, which causes problems with reading the file in R etc.
Not sure I see the issue. Currently (mb v pre-3.2.7) I have this (see below). Can you give an example on the header mismatch?
[ID: 3363454712]
Parameter Mean Variance CredInt_Lower CredInt_Upper Median PSRF Nruns
length[1] 2.960091e-02 6.941899e-05 1.551849e-02 4.609698e-02 2.870080e-02 9.993874e-01 2
length[2] 5.683918e-02 1.318669e-04 3.603541e-02 7.957600e-02 5.598001e-02 1.002068e+00 2
length[3] 2.219691e-02 5.806434e-05 7.565149e-03 3.656863e-02 2.147934e-02 1.000370e+00 2
length[4] 5.700314e-02 1.272752e-04 3.619497e-02 7.887455e-02 5.669569e-02 1.003029e+00 2
length[5] 6.226569e-02 1.322869e-04 4.241672e-02 8.682871e-02 6.147461e-02 9.997443e-01 2
length[6] 9.764845e-03 2.070111e-05 2.476778e-03 1.874374e-02 8.935915e-03 9.989104e-01 2
length[7] 2.855160e-02 5.977353e-05 1.370642e-02 4.287688e-02 2.773313e-02 9.989123e-01 2
length[8] 1.327044e-01 3.095640e-04 9.799957e-02 1.650352e-01 1.323964e-01 9.996849e-01 2
length[9] 3.225995e-02 7.523941e-05 1.637728e-02 4.912450e-02 3.138446e-02 9.997713e-01 2
length[10] 1.892585e-02 4.972666e-05 6.993651e-03 3.349350e-02 1.832597e-02 9.989399e-01 2
length[11] 8.867930e-03 2.422600e-05 1.272899e-03 1.820825e-02 8.037385e-03 1.006504e+00 2
length[12] 9.748204e-03 2.794165e-05 1.652509e-03 2.041604e-02 9.141705e-03 9.987066e-01 2
length[13] 9.207161e-03 2.785406e-05 1.501596e-04 1.884100e-02 8.138787e-03 9.987114e-01 2
length[14] 1.007304e-02 2.721849e-05 7.376355e-04 2.001096e-02 9.302668e-03 1.009816e+00 2
length[15] 7.051515e-03 2.186726e-05 3.786873e-04 1.721006e-02 6.208944e-03 9.975316e-01 2
length[16] 8.181499e-03 1.956783e-05 1.267793e-03 1.578982e-02 7.921050e-03 9.950886e-01 2
length[17] 7.365950e-03 1.790256e-05 1.713367e-04 1.454458e-02 6.807642e-03