marginals: Extract Marginal Models

In lcmix: Layered and chained mixture models

Description

Extract the marginal models from an object of class mdmixmod.

Usage

## S3 method for class 'mdmixmod'
marginals(x, ...)

Arguments

x	an object of class mdmixmod.
...	currently unused.

Details

marginals fits and returns the marginal models from x, that is, objects of class mixmod for each of the data sources used in the fitting of x, with the corresponding distribution.

Value

A list.

Author(s)

Daniel Dvorkin

See Also

mixmod, mdmixmod, posteriors

Examples

## Not run: 
data(CiData)
joint <- mdmixmod(CiData, c(2,3,2), family=c("pvii", "norm", "norm"),
    topology="chained")
marginals(joint)
# $binding
# PVII mixture model ('pvii')
# Data 'binding' of length 10244 fitted to 2 components
# Model statistics:
#      iter      llik      qval       bic    iclbic 
#    353.00 -14409.37 -20330.50 -28883.38 -40725.65 
# 
# $expression
# Normal mixture model ('mvnorm')
# Data 'expression' of size 10244-by-4 fitted to 3 components
# Model statistics:
#       iter       llik       qval        bic     iclbic 
#      42.00  -47499.54  -50052.71  -95405.40 -100511.73 
# 
# $conservation
# PVII mixture model ('pvii')
# Data 'conservation' of length 10244 fitted to 2 components
# Model statistics:
#       iter      llik      qval       bic    iclbic 
#     897.00 -14259.84 -17797.79 -28584.32 -35660.23

## End(Not run)

lcmix documentation built on May 2, 2019, 6:49 p.m.