map.jmoment {mapfit}

Joint moments of MAP

Description

Compute joint moments for a given MAP

Usage

map.jmoment(lag, map, ...)

Arguments

lag

An integer for lag

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.jmoment(lag=1, map=param1)
map.jmoment(lag=1, map=param2)