Mb-sump-output

From Earlham CS Department
Revision as of 12:39, 14 April 2008 by Charliep (talk | contribs)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search
MrBayes > sump filename=Testudinoidea_all-1.nex.mb nruns=1

   Setting sump filename to Testudinoidea_all-1.nex.mb
   Setting sump output file to Testudinoidea_all-1.nex.mb.stat
   Setting sump nruns to 1
   Summarizing parameters in file Testudinoidea_all-1.nex.mb.p
   Writing output to screen but not to file ('Printtofile = No')
   UNIX line termination
   Longest line length = 124
   Found 5001 parameter lines in file "Testudinoidea_all-1.nex.mb"
   All 5001 lines will be summarized (starting at line 3)
   (Only the last set of lines will be read, in case multiple
   parameter blocks are present in the same file.)
   5001 rows and 14 columns in each row
   Successfully read 5001 lines from last parameter block

   Below is a rough plot of the generation (x-axis) versus the log   
   probability of observing the data (y-axis). You can use this      
   graph to determine what the burn in for your analysis should be.  
   When the log probability starts to plateau you may be at station- 
   arity. Sample trees and parameters after the log probability      
   plateaus. Of course, this is not a guarantee that you are at sta- 
   analysis should be. When the log probability starts to plateau    
   tionarity. When possible, run multiple analyses starting from dif-
   ferent random trees; if the inferences you make for independent   
   analyses are the same, this is reasonable evidence that the chains
   have converged. You can use MrBayes to run several independent    
   analyses simultaneously. During such a run, MrBayes will monitor  
   the convergence of topologies. After the run has been completed,  
   the 'sumt' and 'sump' functions will provide additional conver-   
   gence diagnostics for all the parameters in your model. Remember  
   that the burn in is the number of samples to discard. There are   
   a total of ngen / samplefreq samples taken during a MCMC analysis.

   +------------------------------------------------------------+ -50896.09
   | ***********************************************************|
   |                                                            |
   |                                                            |
   |                                                            |
   |                                                            |
   |                                                            |
   |                                                            |
   |                                                            |
   |                                                            |
   |                                                            |
   |                                                            |
   |                                                            |
   |                                                            |
   |                                                            |
   |*                                                           |
   +------+-----+-----+-----+-----+-----+-----+-----+-----+-----+ -52237.00
   ^                                                            ^
   1                                                            5000000


   Estimated marginal likelihoods for run sampled in file "Testudinoidea_all-1.nex.mb.p":
      (Use the harmonic mean for Bayes factor comparisons of models)

   Arithmetic mean   Harmonic mean
   --------------------------------
      -50878.70 *       -87193.62 *
   --------------------------------
   * These estimates may be unreliable because 
     some extreme values were excluded


   Model parameter summaries for run sampled in file "Testudinoidea_all-1.nex.mb":
      (Based on a total of 5001 samples out of a total of 5001 samples from this analysis)

                                              95% Cred. Interval
                                            ----------------------
   Parameter       Mean        Variance       Lower         Upper         Median
   ------------------------------------------------------------------------------
   TL            5.363889      0.364596      5.104000      5.576000      5.328000
   r(A<->C)      0.052355      0.000010      0.047490      0.057619      0.052272
   r(A<->G)      0.224344      0.000102      0.206142      0.243204      0.224514
   r(A<->T)      0.053276      0.000011      0.047998      0.059155      0.053220
   r(C<->G)      0.017005      0.000007      0.013740      0.020402      0.016936
   r(C<->T)      0.629575      0.000189      0.609984      0.650019      0.629242
   r(G<->T)      0.023446      0.000011      0.018576      0.028516      0.023368
   pi(A)         0.329686      0.000039      0.317861      0.342043      0.329402
   pi(C)         0.278162      0.000031      0.268185      0.288695      0.278227
   pi(G)         0.171375      0.000031      0.160837      0.181418      0.171369
   pi(T)         0.220777      0.000023      0.211556      0.229830      0.220873
   alpha         0.229991      0.000050      0.218927      0.241060      0.229703
   ------------------------------------------------------------------------------