Simulatore di catene di Markov

Share:

Description

Questa funzione simula una catena di Markov a stati finiti.

Usage

1
Markov(x0, n, x, P)

Arguments

x0

stato iniziale

n

lunghezza della traiettoria

x

insieme degli stati

P

matrice di probabilita' di transizione

Details

La funzione Markov2 e' basata sulla funzione sample.

Value

Una lista contente la traiettoria della catena di Markov:

X

valori assunti dalla catena di Markov

t

tempi

See Also

gen.vc,trajectory,lewis.

Examples

1
2
3
4
5
6
7
8
9
x <- c("P","S","N")
P <- matrix(c(0.5,0.5,0.25,0.25,0,0.25,0.25,0.5,0.5),3,3)
Markov("S",15,x,P)  -> traj
traj
plot(traj$t,unclass(factor(traj$X)),type="s",axes=FALSE,
 xlab="t",ylab="Che tempo fa'")
axis(1)
axis(2,c(1,2,3),levels(factor(traj$X)))
box()