hzar.make.obsDataGroup: Collect optimizer output from multiple models for analysis.
In hzar: Hybrid Zone Analysis using R

Description Usage Arguments Value Author(s) See Also Examples

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.

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

`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.

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 by`data.groups`

Graham Derryberry asterion@alum.mit.edu

hzar.sameModel hzar.obsData hzar.dataGroup

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));