map.mmoment {mapfit}

Marginal moments of MAP

Description

Compute up to k-th marginal moments for a given MAP

Usage

map.mmoment(k, map, ...)

Arguments

k

An integer for the moments to be computed

map

An instance of MAP

...

Others

Value

A vector of moments

Examples

## create an MAP with specific parameters
(param1 <- map(alpha=c(1,0,0),
               D0=rbind(c(-4,2,0),c(2,-5,1),c(1,0,-4)),
               D1=rbind(c(1,1,0),c(1,0,1),c(2,0,1))))

## create an ER-HMM with specific parameters
(param2 <- erhmm(shape=c(2,3), alpha=c(0.3,0.7),
                 rate=c(1.0,10.0),
                 P=rbind(c(0.3, 0.7), c(0.1, 0.9))))

map.mmoment(k=3, map=param1)
map.mmoment(k=3, map=param2)