MrBayes icon indicating copy to clipboard operation
MrBayes copied to clipboard

Error in sumt vstat file

Open ronquist opened this issue 7 years ago • 1 comments

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.

ronquist avatar Dec 04 '17 18:12 ronquist

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

nylander avatar Feb 13 '19 13:02 nylander