Überprüfen Sie die memorylose Eigenschaft einer Markov-Kette

17

Ich vermute, dass eine Reihe von beobachteten Sequenzen eine Markov-Kette sind ...

X=(ACDDBACBAACADABCADABE)

Wie kann ich jedoch überprüfen, ob sie tatsächlich die memorylose Eigenschaft von respektieren

P(Xi=xi|Xj=xj)?

Oder zumindest beweisen, dass sie Markov in der Natur sind? Beachten Sie, dass dies empirisch beobachtete Sequenzen sind. Irgendwelche Gedanken?

BEARBEITEN

Um nur hinzuzufügen, das Ziel ist es, einen vorhergesagten Satz von Sequenzen mit den beobachteten zu vergleichen. Wir würden uns daher über Kommentare freuen, wie diese am besten verglichen werden können.

Übergangsmatrix erster Ordnung

Mij=xijmxik
wobei m = A..E ist

M=(0.18340.30770.07690.14790.28400.46970.11360.00760.25000.15910.18270.24040.22120.19230.16350.23780.18180.06290.33570.18180.24580.17880.11730.17880.2793)

Eigenwerte von

E=(1.0000000000.2283000000.1344000000.11360.0430i000000.1136+0.0430i)

Eigenvektoren von M

V=(0.44720.58520.42190.23430.0421i0.2343+0.0421i0.44720.78380.42110.44790.2723i0.4479+0.2723i0.44720.20060.37250.63230.63230.44720.00100.70890.21230.0908i0.2123+0.0908i0.44720.05400.05890.2546+0.3881i0.25460.3881i)
HCAI
quelle
Die Spalten enthalten die Reihen und die Zeilen die Elemente der Sequenzen? Wie viele Zeilen und Spalten wurden beobachtet?
mpiktas
2
Mögliches Duplikat: stats.stackexchange.com/questions/29490/…
mpiktas
@mpiktas Die Zeilen repräsentieren die unabhängig beobachteten Folgen von Übergängen durch die Zustände AD. Es gibt ungefähr 400 Sequenzen ... Beachten Sie, dass die beobachteten Sequenzen nicht alle gleich lang sind. Tatsächlich wird die obige Matrix in vielen Fällen durch Nullen ergänzt. Danke übrigens für den Link. Es scheint, dass es auf diesem Gebiet noch viel Raum für Arbeit gibt. Hast du noch weitere Gedanken? Grüße,
HCAI
1
Die lineare Regression war ein Beispiel, um den Punkt meiner Argumentation zu stärken. Dh, dass Sie die Markov-Eigenschaft möglicherweise nicht direkt testen müssen, müssen Sie nur ein Modem anpassen, das die Markov-Eigenschaft übernimmt, und dann die Gültigkeit des Modells überprüfen.
mpiktas
1
Ich erinnere mich vage, dass ich irgendwo einen Hypothesentest für H0 = {Markov} vs H1 = {Markov order 2} gesehen habe. Das könnte helfen.
Stéphane Laurent

Antworten:

5

Ich frage mich, ob das Folgende einen gültigen Pearson Test für die folgenden Proportionen ergeben würde.χ2

  1. Schätzen Sie die Ein-Schritt-Übergangswahrscheinlichkeiten ab - das haben Sie getan.
  2. Erhalten , um die Zwei-Schritt - Modellwahrscheinlichkeiten: p U , V = P r o b [ X i + 2 = U | X i = V ] = Σ W { A , B , C , D } P R o b [ X i + 2 = U | X i + 1 = W ] P r
    p^U,V=Prob[Xi+2=U|Xi=V]=W{A,B,C,D}Prob[Xi+2=U|Xi+1=W]Prob[Xi+1=W|Xi=V]
  3. Erhalten Sie die zweistufigen empirischen Wahrscheinlichkeiten
    p~U,V=i#Xi=V,Xi+2=Ui#Xi=V
  4. TV=#{Xi=V}U(p^U,Vp~U,V)2p^U,V,T=TA+TB+TC+TD

TUχ32Tχ122p^p¯

StasK
quelle
Müssen die Wahrscheinlichkeiten nicht eine Normalverteilung mit Mittelwert 0 und Varianz = 1 haben, damit dies gilt? Es würde mich sehr interessieren, was hier jemand denkt.
HCAI
Das ist, was die Ausdrücke in der Summe sein sollen, asymptotisch mit großen Zählimpulsen.
StasK
6

Die Markov-Eigenschaft ist möglicherweise nur schwer direkt zu testen. Es kann jedoch ausreichend sein, ein Modell anzupassen, das die Markov-Eigenschaft annimmt, und dann zu testen, ob das Modell gültig ist. Es kann sich herausstellen, dass das angepasste Modell eine gute Annäherung ist, die für Sie in der Praxis nützlich ist, und Sie müssen sich keine Sorgen machen, ob die Markov-Eigenschaft tatsächlich gilt oder nicht.

Die Parallele kann zur linearen Regression gezogen werden. In der Regel wird nicht geprüft, ob die Linearität gilt, sondern ob das lineare Modell eine nützliche Annäherung darstellt.

mpiktas
quelle
Dies scheint die beste Option in der Realität zu sein, nur kann ich ein lineares Modell nicht mit tatsächlichen experimentellen Daten vergleichen. Oder hattest du etwas anderes im Sinn?
HCAI
6

Um den Vorschlag der vorherigen Antwort zu konkretisieren, möchten Sie zunächst die Markov-Wahrscheinlichkeiten abschätzen - vorausgesetzt, es handelt sich um Markov. Siehe die Antwort hier Schätzen der Markov-Kettenwahrscheinlichkeiten

MM2M2

MM2

Eine andere Möglichkeit wäre zu sehen, ob die Proportionen des Grundzustands: die in A verbrachte Proportionszeit, die in B verbrachte Zeit, mit dem Eigenvektor des Einheitseigenwerts von M übereinstimmen Staat sollte zu dieser Grenze neigen.

Placidia
quelle
MM2
Der letztere Kommentar ist auch sehr interessant, obwohl ich nicht die Zeit habe, die in jedem Zustand meiner beobachteten Sequenzen verbracht wird. Ich habe nur die Gesamtzeit für jede Zeile. Dies kann die Anwendbarkeit dieser Methode einschränken. Was sind deine Gedanken?
HCAI
1
M2
RE: Gleichgewicht. Ich ging davon aus, dass die Übergänge zu festgelegten Zeitpunkten stattfinden - sagen wir, Sie wechseln jede Sekunde vom aktuellen zum nächsten Zustand. Sie können die Häufigkeit von A-, B-, C- und D-Zuständen in der Nähe der Enden der Sequenzen oder über Sequenzen hinweg bestimmen, um das Grenzverhalten abzuschätzen.
Placidia,
Wenn Sie in R eigen (M) tun, sollten Sie die Eigenwerte und Eigenvektoren von M erhalten. Ein Eigenwert ist 1. Der entsprechende Eigenvektor sollte proportional zu Ihren stationären Proportionen sein .... wenn Markov.
Placidia,
2

XtP(t)ttXtXt7Xt1

XtXt2Xt1nXtXt2{Xt1=xj}nxjXt>1Xt2

applysweepp(Xt|Xt1=xj,Xt2=xi)ij als spaltenindex im spalier soll unter MP zu ähnlichen verteilungen innerhalb einer spalte führen.

Der Junge. 5 des Buches Die statistische Analyse stochastischer Prozesse in der Zeit von JK Lindsey enthält andere Ideen zur Überprüfung von Annahmen.

enter image description here

[## simulates a MC with transition matrix in 'trans', starting from 'ini'
simMC <- function(trans, ini = 1, N) {
  X <- rep(NA, N)
  Pcum <- t(apply(trans, 1, cumsum))
  X[1] <- ini 
  for (t in 2:N) {
    U <- runif(1)
    X[t] <- findInterval(U, Pcum[X[t-1], ]) + 1
  }
  X
}
set.seed(1234)
## transition matrix
P <- matrix(c(0.1, 0.1, 0.1, 0.7,
              0.1, 0.1, 0.6, 0.2,
              0.1, 0.3, 0.2, 0.4,
              0.2, 0.2, 0.3, 0.3),
            nrow = 4, ncol = 4, byrow = TRUE)
N <- 2000
X <- simMC(trans = P, ini = 1, N = N)
## it is better to work with factors
X <- as.factor(X)
levels(X) <- LETTERS[1:4]
## table transitions and normalize each row
Phat <- table(X[1:(N-1)], X[2:N])
Phat <- sweep(x = Phat, MARGIN = 1, STATS = apply(Phat, 1, sum), FUN = "/")
## explicit dimnames
dimnames(Phat) <- lapply(list("X(t-1)=" ,"X(t)="),
                         paste, sep = "", levels(as.factor(X)))
## transition 3-fold contingency array
P3 <- table(X[1:(N-2)], X[2:(N-1)], X[3:N])
dimnames(P3) <- lapply(list("X(t-2)=", "X(t-1)=" ,"X(t)="),
                       paste, sep = "", levels(as.factor(X)))
## apply ONE indendence test 
fisher.test(P3[ , 1, ], simulate.p.value = TRUE)
## plot conditional distr.
library(lattice)
X3 <- data.frame(X = X[3:N], lag1X =  X[2:(N-1)], lag2X = X[1:(N-2)])
histogram( ~ X | lag1X + lag2X, data = X3, col = "SteelBlue3")

]

Yves
quelle
2

Ich denke, Placida und Mpiktas haben beide sehr nachdenkliche und ausgezeichnete Ansätze gegeben.

Ich antworte, weil ich nur hinzufügen möchte, dass man einen Test erstellen könnte, um zu sehen, ob P(Xich=x|Xich-1=y) unterscheidet sich von P(Xich=x|Xich-1=y und Xich-2=z).

Ich würde Werte für auswählen x, y und z für die es eine große Anzahl von Fällen gibt, in denen der Übergang von z zu y zu xtritt ein. Berechnen Sie Stichprobenschätzungen für beide Wahrscheinlichkeiten. Testen Sie dann die Proportionen. Der schwierige Aspekt dabei ist, die Varianzen der beiden Schätzungen unter der Nullhypothese zu ermitteln, die besagen, dass die Proportionen gleich sind und die Kette stationär und Markov ist. Wenn wir in diesem Fall unter der Nullhypothese nur alle 2-Stufen-Übergänge betrachten und sie mit ihren entsprechenden 3-Stufen-Übergängen vergleichen, aber nur Ergebnisse einschließen, bei denen diese gepaarten Ergebnissätze um mindestens 2 Zeitpunkte voneinander getrennt sind, dann ist die Folge der gemeinsamen Ergebnisse wo Erfolg ist definiert alsz zu y zu x Übergang und alle anderen zweistufigen Übergänge zu xals Misserfolge repräsentieren eine Reihe von unabhängigen Bernoulli-Versuchen unter der Nullhypothese. Dasselbe würde für die Definition aller funktioniereny zu x Übergänge als Erfolge und andere einstufige Übergänge zu x als Fehlschläge.

Then the test statistic would be the difference between these estimated proportions. The complication to the standard comparison of the Bernoulli sequences is that they are correlated. But you could do a bootstrap test of binomial proportions in this case.

The other possibility is to construct a two by two table of the two stage and three stage paired outcomes where 0 is failure and 1 is success and the cell frequencies are counts for the pairs (0,0), (0,1), (1,0) and (1,1) where the first component is the two stage outcome and the second is the corresponding three stage outcome. You can then apply McNemar's test to the table.

Michael R. Chernick
quelle
I see what you are referring to here although I'm finding the first paragraph very terse however. For example "Compute sample estimates[...], then test for difference in proportions". What do you mean by sample estimates? Surely there would be no variance in
P(Xi|Xi1=y)
or am I misunderstanding your train of thought?
HCAI
@user1134241 You mentioned "empirically observed", I assumed that you have data from this stochastic sequence. If you want to estimate P(Xi=x|Xi1=y) for each index i-1 where Xi1=y, count the number of times Xi = x and divide it by the number of times Xi1 = y (regardless of what Xi equals). That is an estimate because the observed finite sequence is just a sample of a portion of a sequence of the stochastic process.
Michael R. Chernick
Lassen Sie mich in Ihrem letzten Absatz fragen, was genau einen Erfolg ausmacht. Wenn Sie einen zweistufigen Übergang sagen: Sagen Sie das?ichjich und ein 3-stufiger wäre ichjkich?
HCAI
1

You could bin the data into evenly spaced intervals, then compute the unbiased sample variances of subsets {Xn+1:Xn=x1,Xnk=x2}. By the law of total variance,

Var[E(Xn+1|Xn,Xnk)|Xn]=Var[Xn+1|Xn]E(Var[Xn+1|Xn])

The LHS, if it is almost zero, provides evidence that the transition probabilities do not depend on Xnk, though it is clearly a weaker statement: e.g., let Xn+1N(Xn,Xn1). Taking the expected value of both sides of the above equation, the RHS can be computed from the sample variances (i.e., replacing expected values with averages). If the expected value of the variance is zero then the variance is 0 almost always.

Luke O'Connor
quelle