InitialLawSM {SMM}

Estimation of the initial law (semi-Markov model)

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)