“Önceki durumun” R'deki “sonraki durum” üzerinde etkisi olup olmadığı nasıl test edilir?

Bir durum düşünün: Üç mayının tarihsel kayıtlarımız (20 yıl) var. Gümüşün varlığı gelecek yıl altın bulma olasılığını artırıyor mu? Böyle bir soru nasıl test edilir?

resim açıklamasını buraya girin

Örnek veriler:

mine_A <- c("silver","rock","gold","gold","gold","gold","gold",
            "rock","rock","rock","rock","silver","rock","rock",
            "rock","rock","rock","silver","rock","rock")
mine_B <- c("rock","rock","rock","rock","silver","rock","rock",
            "silver","gold","gold","gold","gold","gold","rock",
            "silver","rock","rock","rock","rock","rock")
mine_C <- c("rock","rock","silver","rock","rock","rock","rock",
            "rock","silver","rock","rock","rock","rock","silver",
            "gold","gold","gold","gold","gold","gold")
time <- seq(from = 1, to = 20, by = 1)

r time-series hypothesis-testing stochastic-processes

— Ladislav Naďo
kaynak

Geçiş matrislerini hesaplamak ilginizi çekebilir .

— Andy W

Merhaba @AndyW! Yorum için teşekkürler. Geçiş matrisleri paketini tanıyorum: makkovchain - markovchainFit (). Geçiş matrisindeki olasılık değerlerini p-değerleri olarak kullanabilir miyim? Hipotezin nasıl test edileceğinin bir yolu var mı: "" Gümüş-altın "ilişkisi var." (p-değeri = xx)?

— Ladislav Naďo

@LadislavNado geçiş olasılıkları p-değerleri olarak yorumlanamaz (herhangi bir H0'ı reddetme hakkında hiçbir şey söylemezler), p-değerleri hakkında daha fazla bilgi için bkz. Stats.stackexchange.com/questions/31/… .

— Tim

Verilerinizi alma yönteminizle ilgili bir sorun görüyorum. "Gümüş: hayır" & "altın: evet" senaryolarınızı düşünün, bu mantık kriterlerini karşıladığı için ardışık "altın" işlemlerinizi de saymalısınız.

Bir hücre 1'den 14'e düzeltildiğinde, model şu şekilde değişir: Katsayılar: Tahmini Std. Hata z değeri Pr (> | z |) (Kesme noktası) -1.2528 0.8018 -1.562 0.118 as.Faktör (c (0, 1)) 1 0.3655 0.8624 0.424 0.672

En iyi denemem: ... @AndyW tarafından önerilen geçiş matrislerinin kullanımı muhtemelen aradığım çözüm değil (@ Tim'in yorumuna dayanarak). Bu yüzden farklı bir yaklaşım denedim. Yanıt değişkeni y ve bir tahmin değişkeni x'in ikilik olduğu lojistik regresyonun nasıl yapılacağı ile ilgili bu bağlantıyı buldum .

Örneğe göre verilerime dayanarak 2x2 tablo oluşturmalıyım:

               gold (yes)  gold (no)
silver (yes)       2           7
silver (no)       14          34

Nasıl değerleri çıkardı: resim açıklamasını buraya girin

Ve bir model oluşturun:

response <- cbind(yes = c(2, 14), no = c(7, 34))

mine.logistic <- glm(response ~ as.factor(c(0,1)),
                      family = binomial(link=logit))

summary(mine.logistic)
# Coefficients:
#                     Estimate Std. Error z value Pr(>|z|)
# (Intercept)          -1.2528     0.8018  -1.562    0.118
# as.factor(c(0, 1))1   0.3655     0.8624   0.424    0.672

İyi bir çözüm mü? P-değeri (0.673) gümüşün varlığının altın bulma olasılığını artırmayacağı anlamına mı geliyor?

— Ladislav Naďo
kaynak

Bu güzel grafikleri nasıl oluşturdunuz? TikZ?

— shadowtalker

Merhaba @ssdecontrol! Grafikler Inkscape'te elle yapıldı.

— Ladislav Naďo

Evet, bu iyi bir yorum. Ayrıca, sadece 2x2 tablonuzun satırlarına bakarsanız, üst sırada (gümüş: evet), 2'si altın olan 9 dava var, bu nedenle gelecek yıl gümüş altın olasılığı verildiğinde 2/9 = 0.222'dir. En alt sırada (gümüş: hayır), 14'ü gelecek yıl altın olan 48 vakanız var, bu nedenle gümüş verilmemesi altın 14 / (14 + 34) = 0.292'dir. Tüm bunlar göz önüne alındığında , p değerinizden "istatistiksel olarak önemli" olmasa da, gümüşün altın bulma şansınıza zarar verdiği görülüyor .

— Gregor Thomas

Ayrıca kodlamanıza dikkat edin, başlarsınız yes = c(2, 14), no = c(7, 34), yani Gümüş koymanız gerekir: evet önce. Böylece as.factor(c(0, 1))0 yaptığınız zaman gümüşe karşılık gelir: evet, bu sizin referans seviyeniz ve dolayısıyla da kesişmenizdir. 0.67 p değeri, gümüşten gümüşe doğru hareket eden altın bulma olasılığınız olan küçük pozitif yumruya karşılık gelir: evet, gümüş: hayır.

— Gregor Thomas

Son bir açıklama: Eğer edilmektedir geçiş matrisleri kullanarak. 2, 7, 14, 34 matrisiniz bir geçiş matrisidir.

— Gregor Thomas