est = new MCMC(new MetropolisHastings(like,like,prop));
est.setChainLength(10); //Number of samples to keep in memory...
est.setMaxIterations(300);
String[] parameters = like.parameterSet();
MCMCBlock varcompBlock = new MCMCBlock();
MCMCBlock otherBlock = new MCMCBlock();
varcompBlock.setName("Latent pH/Variance components");
otherBlock.setName("Model components");
varcompBlock.setCycles(25);
otherBlock.setCycles(25);
for (int p=0;p<parameters.length;p++)
{
if (parameters[p].contains("sigma"))
{
varcompBlock.addParameter(parameters[p]);
}
else if (parameters[p].startsWith("pH"))
{
varcompBlock.addParameter(parameters[p]);
}
else if (parameters[p].startsWith("linearpH"))
{
varcompBlock.addParameter(parameters[p]);
}
else if (parameters[p].startsWith("quadraticpH"))
{
varcompBlock.addParameter(parameters[p]);
}
else if (parameters[p].startsWith("cubicpH"))
{
varcompBlock.addParameter(parameters[p]);
}
else
{
otherBlock.addParameter(parameters[p]);
}
}
est.addBlock(varcompBlock);
est.addBlock(otherBlock);