Periksa properti tanpa memori dari rantai Markov

17

Saya menduga bahwa serangkaian urutan yang diamati adalah rantai Markov ...

X=(ACDDBACBAACADABCADABE)

Namun bagaimana saya bisa mengecek bahwa mereka benar-benar menghargai properti tanpa memori

P(Xi=xi|Xj=xj)?

Atau paling tidak membuktikan bahwa mereka adalah Markov di alam? Perhatikan ini adalah urutan yang diamati secara empiris. Adakah pikiran?

EDIT

Hanya untuk menambahkan, tujuannya adalah untuk membandingkan serangkaian urutan yang diprediksi dari yang diamati. Jadi kami sangat menghargai komentar tentang bagaimana cara terbaik untuk membandingkan ini.

Matriks Transisi Pesanan Pertama

Mij=xijmxik
mana m = A..E menyatakan

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)

Nilai eigen dari M

E=(1.0000000000.2283000000.1344000000.11360.0430i000000.1136+0.0430i)

Vektor eigen dari 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
sumber
Kolom berisi seri, dan baris elemen dari rangkaian? Berapa jumlah baris dan kolom yang diamati?
mpiktas
2
Kemungkinan rangkap: stats.stackexchange.com/questions/29490/…
mpiktas
@mpiktas Baris mewakili urutan transisi yang diamati secara independen melalui status AD. Ada sekitar 400 urutan ... Ingatlah bahwa urutan yang diamati tidak semuanya sama panjang. Bahkan matriks di atas dalam banyak kasus ditambah dengan nol. Omong-omong, terima kasih atas tautannya. Tampaknya masih ada ruang yang cukup untuk bekerja di bidang ini. Apakah Anda memiliki pemikiran lebih lanjut? Salam,
HCAI
1
Regresi linier adalah contoh untuk memperkuat titik argumen saya. Yaitu bahwa Anda mungkin tidak perlu menguji properti Markov secara langsung, Anda hanya perlu memasukkan beberapa modem yang mengasumsikan properti Markov dan kemudian memeriksa validitas model.
mpiktas
1
Samar-samar saya ingat saya telah melihat suatu tempat tes hipotesis untuk H0 = {Markov} vs H1 = {Markov order 2}. Ini bisa membantu.
Stéphane Laurent

Jawaban:

5

Saya ingin tahu apakah yang berikut akan memberikan tes Pearson valid untuk proporsi sebagai berikut.χ2

  1. Perkirakan probabilitas transisi satu langkah - Anda telah melakukannya.
  2. Mendapatkan model probabilitas dua
    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. Dapatkan probabilitas empiris dua langkah
    p~U,V=i#Xi=V,Xi+2=Ui#Xi=V
  4. Form Pearson uji statistik
    TV=#{Xi=V}U(p^U,Vp~U,V)2p^U,V,T=TA+TB+TC+TD

Hal ini menggoda bagi saya untuk berpikir bahwa setiap , sehingga total T ~ χ 2 12 . Namun, saya tidak sepenuhnya yakin akan hal itu, dan akan menghargai pemikiran Anda tentang ini. Saya tidak juga tidak co sertain tentang apakah salah satu kebutuhan untuk menjadi paranoid tentang kemerdekaan, dan ingin membagi sampel di bagian untuk memperkirakan p dan ˉ p .TUχ32Tχ122p^p¯

Tugas
sumber
Tidakkah probabilitas harus memiliki distribusi normal dengan rata-rata 0 dan varians = 1 untuk dapat bertahan? Saya akan sangat tertarik untuk mengetahui apa yang dipikirkan orang di sini.
HCAI
Seharusnya istilah-istilah dalam penjumlahan itu asimtotik dengan jumlah besar.
Tugas
6

Properti Markov mungkin sulit untuk diuji secara langsung. Tetapi mungkin cukup untuk mencocokkan model yang mengasumsikan properti Markov dan kemudian menguji apakah model itu berlaku. Mungkin ternyata model yang cocok adalah pendekatan yang baik yang berguna untuk Anda dalam praktik, dan Anda tidak perlu khawatir apakah properti Markov benar-benar dimiliki atau tidak.

Paralel dapat ditarik ke regresi linier. Praktik yang biasa bukan untuk menguji apakah linearitas berlaku, tetapi apakah model linier adalah pendekatan yang berguna.

mpiktas
sumber
Ini sepertinya pilihan terbaik dalam kenyataan, hanya saja saya tidak bisa benar-benar membandingkan model linear dengan data eksperimental yang sebenarnya. Atau apakah Anda memikirkan hal lain?
HCAI
6

Untuk mengkonkretkan saran dari jawaban sebelumnya, Anda pertama-tama ingin memperkirakan probabilitas Markov - dengan asumsi itu Markov. Lihat balasannya di sini Memperkirakan Kemungkinan Rantai Markov

Anda harus mendapatkan 4 x 4 matriks berdasarkan proporsi transisi dari negara A ke A, A ke B, dll Panggil matriks ini . M 2 kemudian harus menjadi matriks transisi dua langkah: A ke A dalam 2 langkah, dan seterusnya. Anda kemudian dapat menguji apakah matriks transisi 2 langkah yang diamati Anda mirip dengan M 2 .MM2M2

Karena Anda memiliki banyak data untuk jumlah status, Anda dapat memperkirakan dari satu setengah data dan menguji M 2 menggunakan setengah lainnya - Anda menguji frekuensi yang diamati terhadap probabilitas teoritis multinomial. Itu seharusnya memberi Anda gambaran tentang seberapa jauh Anda.MM2

Kemungkinan lain adalah untuk melihat apakah proporsi kondisi dasar: proporsi waktu yang dihabiskan dalam A, waktu yang dihabiskan dalam B, cocok dengan vektor eigen dari nilai eigen unit M. Jika seri Anda telah mencapai semacam kondisi mapan, proporsi waktu di masing-masing negara harus cenderung ke batas itu.

Placidia
sumber
Ada sedikit untuk diterima di sana. Saya telah menghitung matriks Transisi , tapi saya tidak yakin bagaimana Anda akan menghitung M 2 secara empiris. Bisakah Anda mengklarifikasi hal itu? Salam,MM2
HCAI
Juga, komentar terakhir ini sangat menarik, walaupun saya tidak punya waktu dihabiskan di setiap negara bagian dari urutan yang saya amati. Saya hanya punya waktu total untuk setiap baris. Sehingga dapat membatasi penerapan metode itu. Apa yang kamu pikirkan?
HCAI
1
Lakukan dengan cara yang sama Anda lakukan M, hanya alih-alih melihat transisi tetangga terdekat, (katakanlah, urutan AB), lihat pasangan yang terpisah 2. Jadi, jika subjek beranjak ACB, itu diperhitungkan dalam transisi AB Anda. Begitu juga ABB. Buat matriks di mana item dalam baris i, kolom j berisi transisi i ke j. Kemudian bagi dengan jumlah kolom. Anda ingin jumlah kolom menjadi 1. Di bawah properti Markov, matriks ini harus dekat dengan M2
Placidia
RE: keseimbangan. Saya berasumsi bahwa transisi terjadi pada saat yang ditetapkan - katakan setiap detik, Anda transisi dari keadaan saat ini ke keadaan berikutnya. Anda bisa mengambil frekuensi status A, B, C, dan D di dekat ujung urutan, atau melintasi urutan untuk memperkirakan perilaku batas.
Placidia
Dalam R, jika Anda melakukan eigen (M), Anda harus mendapatkan nilai eigen dan vektor eigen dari M. Satu nilai eigen akan menjadi 1. Vektor eigen yang sesuai harus proporsional dengan proporsi kondisi stabil Anda .... jika Markov.
Placidia
2

Di luar Markov Property (MP), properti lebih lanjut Waktu Homogenitas (TH): dapat Markov tetapi dengan matriks transisi P ( t ) tergantung pada waktu t . Misalnya, mungkin tergantung pada hari kerja di t jika pengamatan harian, dan kemudian ketergantungan X t pada X t - 7 syarat X t - 1 dapat didiagnosis jika TH adalah terlalu diasumsikan.XtP(t)ttXtXt7Xt1

Dengan asumsi TH memegang, cek mungkin bagi MP adalah pengujian yang independen dari X t - 2 tergantung pada X t - 1 , sebagai Michael Chernick dan StasK disarankan. Ini dapat dilakukan dengan menggunakan tes untuk tabel kontingensi. Kita dapat membangun n tabel kontingensi X t dan X t - 2 bersyarat pada { X t - 1 = x j } untuk n kemungkinan nilai x jXtXt2Xt1nXtXt2{Xt1=xj}nxj, dan menguji independensi. Ini juga dapat dilakukan dengan menggunakan dengan > 1 sebagai pengganti X t - 2 .Xt>1Xt2

Dalam R, kontingensi tabel atau array mudah diproduksi berkat faktor fasilitas dan fungsi apply, sweep. Gagasan di atas juga dapat dieksploitasi secara grafis. Paket ggplot2 atau kisi dengan mudah menyediakan plot bersyarat untuk membandingkan distribusi kondisional . Misalnya pengaturan saya sebagai indeks baris dan jp(Xt|Xt1=xj,Xt2=xi)ij sebagai indeks kolom dalam terali harus di bawah MP mengarah ke distribusi serupa dalam kolom.

Chap. 5 buku Analisis statistik proses stokastik dalam waktu oleh JK Lindsey berisi ide-ide lain untuk memeriksa asumsi.

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
sumber
2

Saya pikir placida dan mpikta telah memberikan pendekatan yang sangat bijaksana dan sangat baik.

P(Xi=x|Xi1=y)P(Xi=x|Xi1=y and Xi2=z)

xyzzyxzyx transisi dan semua transisi dua tahap ke xkarena kegagalan mewakili serangkaian uji coba Bernoulli independen berdasarkan hipotesis nol. Hal yang sama akan berfungsi untuk mendefinisikan semuay untuk x transisi sebagai keberhasilan dan transisi satu tahap ke x sebagai kegagalan.

Maka statistik uji akan menjadi perbedaan antara estimasi proporsi ini. Komplikasi untuk perbandingan standar dari urutan Bernoulli adalah bahwa mereka berkorelasi. Tetapi Anda bisa melakukan tes bootstrap proporsi binomial dalam kasus ini.

Kemungkinan lainnya adalah membangun dua dua tabel dari dua tahap dan tiga tahap hasil berpasangan di mana 0 adalah kegagalan dan 1 sukses dan frekuensi sel dihitung untuk pasangan (0,0), (0,1), (1,0) dan (1,1)di mana komponen pertama adalah hasil dua tahap dan yang kedua adalah hasil tiga tahap yang sesuai. Anda kemudian dapat menerapkan tes McNemar ke tabel.

Michael R. Chernick
sumber
Saya melihat apa yang Anda maksudkan di sini meskipun saya menemukan paragraf pertama sangat singkat. Misalnya "Hitung taksiran sampel [...], lalu uji perbedaan proporsi". Apa yang Anda maksud dengan perkiraan sampel? Tentunya tidak akan ada perbedaan dalam
P(Xsaya|Xsaya-1=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
In your last paragraph, let me ask what constitute a success and exactly? In the case where you say a two-step transition: are you saying iji and a 3-step would be ijki?
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
sumber