InitialLawSM: Estimation of the initial law (semi-Markov model)

In SMM: Simulation and Estimation of Multi-State Discrete-Time Semi-Markov and Markov Models

Description

For one sequence, estimation of the initial law by computating the limit law of the semi-Markovian chain.

For several sequences, estimation of the initial law by computating the first state frequencies.

Usage

InitialLawSM(E, seq, q)

Arguments

E	Vector of state space
seq	List of sequence(s)
q	Array of size SxSxKmax with the semi-Markov kernel (S = length(E))

Value

init	Vector of the initial distribution

Author(s)

Vlad Stefan Barbu, [e-mail]
Caroline Berard, [e-mail]
Dominique Cellier, [e-mail]
Mathilde Sautreuil, [e-mail]
Nicolas Vergne, [e-mail]

See Also

StationaryLaw, estimSM, simulSM, estimMk, simulMk

Examples

seq = list(c("a","c","c","g","t","a","a","a","a","g","c","t","t","t","g"))
res = estimSM(seq = seq, E = c("a","c","g","t"), distr = "NP")
q = res$q
p = res$Ptrans

InitialLawSM(E = c("a","c","g","t"), seq = seq, q = q)

SMM documentation built on Jan. 31, 2020, 5:07 p.m.