/**
* Le paramètre modelResult ne doit pas être null, il permet d'obtenir les valeurs de MonitorBurnin, NBurnin, NIter.
*/
internal Mcmc(ModelResult modelResult, params string[] chainNamePrefixes) : this(modelResult)
{
foreach (string prefix in chainNamePrefixes)
{
ChainPairCreator.ChainPair pair = this.cpc.GetPair(prefix);
ChainDictionary.Add(pair.Burnin.Label, pair.Burnin);
ChainDictionary.Add(pair.Sample.Label, pair.Sample);
}
}
private Mcmc(ModelResult modelResult)
{
this.burninSize = modelResult.ParentModel.MonitorBurnin ? modelResult.ParentModel.NBurnin : 0;
this.sampleSize = modelResult.ParentModel.NIter;
this.cpc = new ChainPairCreator(burninSize, sampleSize);
this.ChainDictionary = new Dictionary <string, Chain>();
this.quantileDef = new Statistics.Quantile();
this.QuantileHeader = "," + string.Join(",", quantileDef.Probabilities.Select(p => string.Format("{0:P1}", p)));
//chainNo = chainNumber;
}