<?xml version="1.0" standalone="yes"?>
<!-- Generated by genBEASTuniforPriorText.sage --> 
<beast>

  <taxa id="taxa">
<taxon id="tip1"/>
<taxon id="tip2"/>
<taxon id="tip3"/>
<taxon id="tip4"/>
<taxon id="tip5"/>
<taxon id="tip6"/>
  </taxa>

  <newick id="tree">
    (((tip1:1,tip2:1):2,((tip3:1,tip4:1):1,tip5:2):1):1,tip6:4)
  </newick>

    <treeModel id="treeModel">
      <tree idref="tree"/>
      <rootHeight>
	<parameter id="treeModel.rootHeight"/>
      </rootHeight>
      <nodeHeights internalNodes="true">
	<parameter id="treeModel.internalNodeHeights"/>
      </nodeHeights>
      <nodeHeights internalNodes="true" rootNode="true">
	<parameter id="treeModel.allInternalNodeHeights"/>
      </nodeHeights>
    </treeModel>

    <dummyModel id="dummy">
       <parameter idref="treeModel.rootHeight"/>
       <parameter idref="treeModel.internalNodeHeights"/>
       <parameter idref="treeModel.allInternalNodeHeights"/>
    </dummyModel>

    <tmrcaStatistic id="n1">
      <treeModel idref="treeModel"/>
      <mrca>
	<taxa>
<taxon idref="tip1"/>
<taxon idref="tip2"/>
<taxon idref="tip3"/>
<taxon idref="tip4"/>
<taxon idref="tip5"/>
	</taxa>
      </mrca>
    </tmrcaStatistic>

    <tmrcaStatistic id="n2">
      <treeModel idref="treeModel"/>
      <mrca>
	<taxa>
<taxon idref="tip1"/>
<taxon idref="tip2"/>
	</taxa>
      </mrca>
    </tmrcaStatistic>

    <tmrcaStatistic id="n5">
      <treeModel idref="treeModel"/>
      <mrca>
	<taxa>
<taxon idref="tip3"/>
<taxon idref="tip4"/>
<taxon idref="tip5"/>
	</taxa>
      </mrca>
    </tmrcaStatistic>

    <tmrcaStatistic id="n6">
      <treeModel idref="treeModel"/>
      <mrca>
	<taxa>
<taxon idref="tip3"/>
<taxon idref="tip4"/>
	</taxa>
      </mrca>
    </tmrcaStatistic>

    <operators id="operators">
      <treeUniform count="3" weight="1">
	<treeModel idref="treeModel"/>
      </treeUniform>
    </operators>

    <mcmc id="mcmc" chainLength="1000000" autoOptimize="true" preBurnin="0">
      <posterior id="posterior">
	<prior id="prior">
	  <uniformPrior lower="0" upper="10">
	    <parameter idref="treeModel.rootHeight"/>
	  </uniformPrior>
	</prior>
	<dummyModel idref="dummy"/>
      </posterior>

      <operators idref="operators"/>

      <log id="fileLog" logEvery="10" fileName="test3.log">
<tmrcaStatistic idref="n1"/>
<tmrcaStatistic idref="n2"/>
<tmrcaStatistic idref="n5"/>
<tmrcaStatistic idref="n6"/>
      </log>
      
      <logTree id="treeFileLog" logEvery="10" nexusFormat="true" fileName="test3.trees">
	<treeModel idref="treeModel"/>
      </logTree>
    </mcmc>

    <traceAnalysis fileName="test3.log" burnIn="0">

      <expectation name="n1" value="3.2"/> <!-- 4*x1^3 --> 
      <expectation name="n2" value="1.6"/> <!-- -(4*x2^3 - 4)/3 --> 
      <expectation name="n5" value="2.13333"/> <!-- 4*x5 - 4*x5^3 --> 
      <expectation name="n6" value="1.06667"/> <!-- (4*x6^3 - 12*x6 + 8)/3 --> 
    </traceAnalysis>

</beast>

