Description Usage Arguments Value Author(s) See Also Examples
View source: R/40-hzarPostProcessing.R
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 |
obsDataGroup |
The |
A hzar.obsDataGroup object.
data.groups |
A list of hzar.dataGroup objects, each with a unique
meta models. See |
obsData |
The |
Graham Derryberry asterion@alum.mit.edu
hzar.sameModel
hzar.obsData
hzar.dataGroup
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 | 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));
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.