Zufälliges Gehen mit Schwung

18

Betrachten Sie einen ganzzahligen Zufallsrundgang ab 0 unter den folgenden Bedingungen:

  • Der erste Schritt ist plus oder minus 1 mit gleicher Wahrscheinlichkeit.

  • Jeder zukünftige Schritt ist: 60% wahrscheinlich in die gleiche Richtung wie der vorherige Schritt, 40% wahrscheinlich in die entgegengesetzte Richtung

Welche Verteilung bringt das?

Ich weiß, dass ein zufälliger Gang ohne Impuls eine Normalverteilung ergibt. Ändert der Impuls nur die Varianz oder die Art der Verteilung vollständig?

Ich suche nach einer generischen Antwort, also meine ich mit 60% und 40% mehr wirklich p und 1-p

barrycarter
quelle
Tatsächlich, @Dilip, brauchen Sie eine Markov-Kette mit Zuständen, die durch geordnete Paare und ( i , i - 1 ) , i Z , indiziert sind . Die Übergänge sind ( i , i + 1 ) ( i + 1 , i + 1 ) und ( i , i - 1 ) ( i - 1 , i(i,i+1)(i,i1)iZ(i,i+1)(i+1,i+1) mit der Wahrscheinlichkeit p und ( i , i + 1 ) ( i + 1 , i ) und ( i , i - 1 ) ( i - 1 , i - 2 ) mit der Wahrscheinlichkeit 1 - p . (i,i1)(i1,i)p(i,i+1)(i+1,i)(i,i1)(i1,i2)1p
whuber
Beachten Sie, dass die Schrittgrößen eine Markov-Kette auf und Sie sie zufällig (?!) bei einer stationären Verteilung gestartet haben. {1,+1}
Kardinal
Sind Sie wollen eine Begrenzung (marginal) Verteilung für , wo die X n{ - 1 , + 1 } sind die Schritte des Weges? Sn=ich=1nXnXn{-1,+1}
Kardinal
Ein anderer Ansatz könnte darin bestehen, sich alternierende Summen geometrischer Zufallsvariablen anzusehen und dann eine Martingaltheorie anzuwenden. Das Problem ist, dass Sie eine Art Stoppzeit definieren müssen, die schwierig sein kann.
Shabbychef

Antworten:

8

Um sofort zum Schluss zu kommen, ändert der "Impuls" nicht die Tatsache, dass die Normalverteilung eine asymptotische Annäherung an die Verteilung des Zufallslaufs ist, sondern die Varianz ändert sich von 4np(1p) zu . Dies kann in diesem speziellen Fall durch relativ elementare Überlegungen abgeleitet werden. Es ist nicht schwierig, die folgenden Argumente auf eine CLT für Markov-Ketten im endlichen Zustandsraum zu verallgemeinern, aber das größte Problem ist tatsächlich die Berechnung der Varianz. Für das jeweilige Problem kann esnp/(1p)berechnet werden, und hoffentlich können die folgenden Argumente den Leser davon überzeugen, dass es sich um die richtige Varianz handelt.

Unter Verwendung der Einsicht, die Cardinal in einem Kommentar liefert, wird der Zufallslauf als wobei X k{ - 1 , 1 } und die X k eine Markov-Kette mit Übergangswahrscheinlichkeitsmatrix bilden ( p 1 - p 1 - p p ) . Aus asymptotischen Überlegungen spielt die Anfangsverteilung von X 1 keine Rolle, wenn n ∞ ist

Sn=k=1nXk
Xk{1,1}Xk
(p1p1pp).
nX1 für das folgende Argument, und nehme auch an, dass 0 < p < 1X1=10<p<1 . Eine raffinierte Technik besteht darin, die Markov-Kette in unabhängige Zyklen zu zerlegen. Sei das erste Mal nach dem Zeitpunkt 1, dass die Markov-Kette zu 1 zurückkehrt. Das heißt, wenn X 2 = 1, dann σ 1 = 2 , und wenn X 2 = X 3 = - 1 und X 4 = 1, dann σ 1 = 4σ1X2=1σ1=2X2=X3=1X4=1σ1=4 . Im Allgemeinen sei die i -te Rückkehrzeit zu 1 und sei τ i = σ i - σ i - 1σiiτi=σiσi1 die Zwischenrückkehrzeiten (mit ). Mit diesen Definitionen haben wirσ0=1
  • Mit ist S σ n = X 1 + n i = 1 U i .Ui=k=σi1+1σiXk
    Sσn=X1+i=1nUi.
  • Da den Wert annimmt , - 1 für k = σ i - 1 + 1 , ... , σ i - 1 und X & sgr; i = 1 gilt , daß U i = 2 - τ i .Xk1k=σi1+1,,σi1Xσi=1
    Ui=2τi.
  • Die Inter-Return-Zeiten für eine Markov-Kette sind iid (formal aufgrund der starken Markov-Eigenschaft) und in diesem Fall mit dem Mittelwert E ( τ i ) = 2 und der Varianz V ( τ i ) = 2 pτiE(τi)=2 . Nachfolgend wird angegeben, wie der Mittelwert und die Varianz berechnet werden.V(τi)=2p1p
  • Die gewöhnliche CLT für iid Variablen ergibt , dass
    SσnasympN(0,2np1p).
  • Das letzte , was zu beachten, die einen kleinen Sprung des Glaubens erfordert, weil ich die Details auslassen, ist , dass , die Ausbeuten , dass S n asymp ~ N ( 0 , n pσn=1+i=1nτi2n
    SnasympN(0,np1p).

τ1P(τ1=1)=pm2P(τ1=m)=(1p)2pm2X1pZ=1(τ1=1)1+X(1Z)τ1

NRH
quelle
+1 schön. Ich hätte nur assymptotische Distribution für geschrieben1/nSn
2

ρρ=2p1

True standard error of x¯p1psn,
nx¯nsn1x¯2nx¯np/(1p)

edit : ich hatte die falsche autokorrelation (bzwphätte anders ausgelegt werden sollen); ist jetzt konsequent (hoffe ich!)

shabbychef
quelle
Interessant. Ich bin mir nicht sicher, ob das etwas sehr Sinnvolles für das ergibtp=0Subcase; Dies könnte jedoch auf mit diesem Fall verbundene Pathologien zurückzuführen sein.
Kardinal
@ Kardinal guter Fang, sollte die Autokorrelation sein ρ=2p-1, nicht 1-2p. Korrekturen ...
Shabbychef