posteriors: Joint And Marginal Model Posterior Probabilities

Description Usage Arguments Details Value Author(s) See Also Examples

Description

Extract the joint and marginal model posterior probabilities from an object of class mdmixmod.

Usage

1
2
## S3 method for class 'mdmixmod'
posteriors(x, ...)

Arguments

x

an object of class mdmixmod.

...

currently unused.

Details

posteriors returns the joint posterior probability for x and the marginal probabilities for each of the data sources used in the fitting of x.

Value

A list. The first element has the name JOINT, indicating the posterior probability for the joint model; the other elements have the names of the data sources used in the fitting of x.

Author(s)

Daniel Dvorkin

See Also

mixmod, mdmixmod, marginals

Examples

 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
## Not run: 
data(CiData)
joint <- mdmixmod(CiData, c(2,3,2), family=c("pvii", "norm", "norm"),
    topology="chained")
pj = posteriors(joint)
lapply(pj, function(p) p[1:3,])
# $JOINT
#           [,1]      [,2]
# [1,] 0.0500650 0.9499350
# [2,] 0.1592164 0.8407836
# [3,] 0.5812072 0.4187928
# 
# $binding
#           [,1]      [,2]
# [1,] 0.1924187 0.8075813
# [2,] 0.2221210 0.7778790
# [3,] 0.8393147 0.1606853
# 
# $expression
#           [,1]         [,2]       [,3]
# 1 0.0003780798 9.543479e-01 0.04527400
# 2 0.3904756217 2.536705e-18 0.60952438
# 3 0.0008795689 9.806697e-01 0.01845077
# 
# $conservation
#           [,1]       [,2]
# [1,] 0.5916767 0.40832326
# [2,] 0.8832577 0.11674231
# [3,] 0.9465522 0.05344779

## End(Not run)

lcmix documentation built on May 31, 2017, 5 a.m.