hzar (version 0.2-5)

hzar.make.obsDataGroup: Collect optimizer output from multiple models for analysis.

Description

This method collects optimizer output based on the same hzar.obsData object. It automatically creates hzar.dataGroup objects as needed, as well as combining any objects which use the same model as determined by hzar.sameModel.

Usage

hzar.make.obsDataGroup(dataGroups, obsDataGroup = NULL)

Arguments

dataGroups

A list of hzar.dataGroup or hzar.fitRequest objects to include.

obsDataGroup

The hzar.obsDataGroup to which dataGroups will be added. If NULL, an empty hzar.obsDataGroup will be added.

Value

A hzar.obsDataGroup object.

data.groups

A list of hzar.dataGroup objects, each with a unique meta models. See hzar.sameModel.

obsData

The hzar.obsData object shared bydata.groups

See Also

hzar.sameModel hzar.obsData hzar.dataGroup

Examples

data(manakinMolecular);
mknAdaA <-
  hzar.doMolecularData1DPops(manakinMolecular$distance,
                             manakinMolecular$ada.A,
                             manakinMolecular$ada.nSamples);
hzar.plot.obsData(mknAdaA);
mknAdaAmodel <-
  hzar.makeCline1DFreq(mknAdaA, scaling="fixed",tails="none");
mknAdaAmodel <-
  hzar.model.addBoxReq(mknAdaAmodel,-30,600);
mknAdaAmodelFitR <-
   hzar.first.fitRequest.old.ML(model=mknAdaAmodel ,
                                mknAdaA,
                                verbose=FALSE);
mknAdaAmodelFitR$mcmcParam$chainLength <- 5e3;
mknAdaAmodelFitR$mcmcParam$burnin <- 5e2;
mknAdaAmodelFit <- hzar.doFit(mknAdaAmodelFitR)
mknAdaAmodelData <-
  hzar.dataGroup.add(mknAdaAmodelFit);
## Not run: 
# mknAdaAmodelFitL <-
#     hzar.chain.doSeq(hzar.next.fitRequest(mknAdaAmodelFit), count=3);
# ## End(Not run)
mknAdaAmodelNull <- hzar.dataGroup.null(mknAdaA);
mknAdaAdGs <- list(clineModel = mknAdaAmodelData,
                   nullModel  = mknAdaAmodelNull);
mknAdaAoDG <- hzar.make.obsDataGroup(mknAdaAmodelFit);
## Not run: mknAdaAoDG <- hzar.make.obsDataGroup(mknAdaAmodelFitL,mknAdaAoDG);
mknAdaAoDG <- hzar.make.obsDataGroup(mknAdaAmodelNull,mknAdaAoDG);

mknAdaAoDG <- hzar.copyModelLabels(mknAdaAdGs,mknAdaAoDG);
hzar.plot.cline(mknAdaAoDG);
print(hzar.AICc.hzar.obsDataGroup(mknAdaAoDG));