createScaleOperator(ClockType.UCGD_SHAPE, demoTuning, rateWeights);
createScaleOperator(ClockType.LOCAL_CLOCK + ".relativeRates", demoTuning, treeWeights);
createOperator(ClockType.LOCAL_CLOCK + ".changes", OperatorType.BITFLIP, 1, treeWeights);
createDiscreteStatistic("rateChanges", "number of random local clocks"); // POISSON_PRIOR
if (USE_DIRICHLET_PRIOR_FOR_MUS) {
<BUG>createNonNegativeParameterDirichletPrior("allMus", "relative rates amongst partitions parameter", this, PriorScaleType.SUBSTITUTION_PARAMETER_SCALE, 1.0);
</BUG>
createOperator("scaleMus", "allMus",
"Scale partition rates relative to each other", "allMus",
OperatorType.SCALE_INDEPENDENTLY, 0.75, 3.0);

createScaleOperator(ClockType.UCGD_SHAPE, demoTuning, rateWeights);
createScaleOperator(ClockType.LOCAL_CLOCK + ".relativeRates", demoTuning, treeWeights);
createOperator(ClockType.LOCAL_CLOCK + ".changes", OperatorType.BITFLIP, 1, treeWeights);
createDiscreteStatistic("rateChanges", "number of random local clocks"); // POISSON_PRIOR
if (USE_DIRICHLET_PRIOR_FOR_MUS) {
createNonNegativeParameterDirichletPrior("allMus", "relative rates amongst partitions parameter", this, PriorScaleType.SUBSTITUTION_PARAMETER_SCALE, 1.0, true);
createOperator("scaleMus", "allMus",
"Scale partition rates relative to each other", "allMus",
OperatorType.SCALE_INDEPENDENTLY, 0.75, 3.0);
