Simulatore di catene di Markov
Questa funzione simula una catena di Markov a stati finiti.
Markov(x0, n, x, P)
x0 |
stato iniziale |
n |
lunghezza della traiettoria |
x |
insieme degli stati |
P |
matrice di probabilita' di transizione |
La funzione Markov2 e' basata sulla funzione sample.
Una lista contente la traiettoria della catena di Markov:
X |
valori assunti dalla catena di Markov |
t |
tempi |
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()Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.