Mam pełny zestaw sekwencji (dokładnie 432 obserwacje) z 4 stanów : np
EDYCJA : Sekwencje obserwacji mają nierówne długości! Czy to coś zmienia?
Czy istnieje sposób obliczania macierzy przejścia
Mam pełny zestaw sekwencji (dokładnie 432 obserwacje) z 4 stanów : np
EDYCJA : Sekwencje obserwacji mają nierówne długości! Czy to coś zmienia?
Czy istnieje sposób obliczania macierzy przejścia
Odpowiedzi:
Sprawdź komentarze powyżej. Oto szybkie wdrożenie w R.
x <- c(1,2,1,1,3,4,4,1,2,4,1,4,3,4,4,4,3,1,3,2,3,3,3,4,2,2,3)
p <- matrix(nrow = 4, ncol = 4, 0)
for (t in 1:(length(x) - 1)) p[x[t], x[t + 1]] <- p[x[t], x[t + 1]] + 1
for (i in 1:4) p[i, ] <- p[i, ] / sum(p[i, ])
Wyniki:
> p
[,1] [,2] [,3] [,4]
[1,] 0.1666667 0.3333333 0.3333333 0.1666667
[2,] 0.2000000 0.2000000 0.4000000 0.2000000
[3,] 0.1428571 0.1428571 0.2857143 0.4285714
[4,] 0.2500000 0.1250000 0.2500000 0.3750000
(Prawdopodobnie głupia) implementacja w MATLAB (z której nigdy nie korzystałem, więc nie wiem, czy to zadziała. Właśnie przejrzałem „zadeklaruj macierz wektorową MATLAB”, aby uzyskać składnię):
x = [ 1, 2, 1, 1, 3, 4, 4, 1, 2, 4, 1, 4, 3, 4, 4, 4, 3, 1, 3, 2, 3, 3, 3, 4, 2, 2, 3 ]
n = length(x) - 1
p = zeros(4,4)
for t = 1:n
p(x(t), x(t + 1)) = p(x(t), x(t + 1)) + 1
end
for i = 1:4
p(i, :) = p(i, :) / sum(p(i, :))
end
for
pętli. Trzecia linia oryginalnego kodu liczy, ile razy przechodzi ze stanu x i do stanu x j ? Jeśli mógłbyś to powiedzieć słowami, bardzo bym to docenił. Pozdrawiam
Oto moja implementacja w języku R.
x <- c(1,2,1,1,3,4,4,1,2,4,1,4,3,4,4,4,3,1,3,2,3,3,3,4,2,2,3)
xChar<-as.character(x)
library(markovchain)
mcX<-markovchainFit(xChar)$estimate
mcX
Oto sposób, aby to zrobić w Matlabie:
x = [1,2,1,1,3,4,4,1,2,4,1,4,3,4,4,4,3,1,3,2,3,3,3,4,2,2,3];
counts_mat = full(sparse(x(1:end-1),x(2:end),1));
trans_mat = bsxfun(@rdivide,counts_mat,sum(counts_mat,2))
Acknowledgement owed to SomptingGuy: http://www.eng-tips.com/viewthread.cfm?qid=236532