Làm thế nào để kiểm tra nếu trạng thái trước đó của Wikipedia có ảnh hưởng đến Bang tiếp theo


10

Hãy tưởng tượng một tình huống: Chúng ta có hồ sơ lịch sử (20 năm) của ba mỏ. Có phải sự hiện diện của bạc làm tăng xác suất tìm thấy vàng trong năm tới? Làm thế nào để kiểm tra câu hỏi như vậy?


nhập mô tả hình ảnh ở đây

Dưới đây là dữ liệu ví dụ:

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)


1
Bạn có thể quan tâm đến việc tính toán ma trận chuyển tiếp .
Andy W

Xin chào @AndyW! Cảm ơn bạn đã bình luận. Tôi quen thuộc với gói ma trận chuyển tiếp: makkovchain - markovchainFit (). Tôi có thể sử dụng các giá trị xác suất từ ​​ma trận chuyển tiếp làm giá trị p không? Có cách nào để kiểm tra giả thuyết: "Tồn tại một mối quan hệ" bạc-vàng "." (giá trị p = xx)?
Ladislav Naďo

1
Xác suất chuyển đổi @LadislavNado không thể được hiểu là giá trị p (chúng không cho bạn biết bất cứ điều gì về việc từ chối bất kỳ H0 nào), xem stats.stackexchange.com/questions/31/ tựa để tìm hiểu thêm về giá trị p.
Tim

1
Tôi thấy một vấn đề với cách bạn đã trích xuất dữ liệu của mình. Hãy xem xét kịch bản "bạc: không" và "vàng: có" của bạn, bạn cũng nên tính số lần chạy "vàng" liên tiếp của mình vì điều đó đáp ứng các tiêu chí logic.

1
Với một ô được sửa từ 1 đến 14, mô hình thay đổi thành: Hệ số: Ước tính Std. Lỗi z giá trị Pr (> | z |) (Chặn) -1.2528 0.8018 -1.562 0.118 as.factor (c (0, 1)) 1 0.3655 0.8624 0.424 0.672

Câu trả lời:


4

Cố gắng tốt nhất của tôi: ... sử dụng ma trận chuyển tiếp được đề xuất bởi @AndyW có lẽ không phải là giải pháp tôi đang tìm kiếm (dựa trên nhận xét của @ Tim). Vì vậy, tôi đã thử một cách tiếp cận khác. Tôi tìm thấy liên kết này liên quan đến cách thực hiện hồi quy logistic trong đó biến trả lời y và biến dự đoán x đều là nhị phân .

Theo ví dụ, tôi nên tạo bảng 2 × 2 dựa trên dữ liệu của mình:

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

Cách tôi trích xuất các giá trị: nhập mô tả hình ảnh ở đây

Và xây dựng một mô hình:

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

Nó có phải là một giải pháp tốt? Giá trị p (0,673) có nghĩa là sự hiện diện của bạc không làm tăng xác suất tìm thấy vàng?


Làm thế nào bạn tạo ra những biểu đồ tốt đẹp? Tikz?
Shadowtalker

Xin chào @ssdecontrol! Biểu đồ được làm bằng tay trong Inkscape.
Ladislav Naďo

Vâng, đó là một cách giải thích hợp lý. Ngoài ra, nếu bạn chỉ nhìn vào các hàng của bảng 2x2, ở hàng trên cùng (bạc: có), bạn có 9 trường hợp, 2 trong số đó có vàng, do đó, xác suất vàng của năm tới là 2/9 = 0,222. Ở hàng dưới cùng (bạc: không) bạn có 48 trường hợp, 14 trong số đó có vàng vào năm tới, do đó không có xác suất bạc nào là vàng là 14 / (14 + 34) = 0,222. Với tất cả những điều đó, có vẻ như bạc làm tổn thương cơ hội tìm vàng của bạn, mặc dù từ giá trị p của bạn không "có ý nghĩa thống kê".
Gregor Thomas

Ngoài ra, hãy chú ý đến mã hóa của bạn, bạn bắt đầu với yes = c(2, 14), no = c(7, 34), điều đó có nghĩa là bạn đặt Bạc: có trước. Vì vậy, khi bạn làm as.factor(c(0, 1))0 tương ứng với bạc: có, đó là mức tham chiếu của bạn và do đó, việc đánh chặn của bạn. Giá trị p 0,67 tương ứng với vết sưng dương nhỏ mà bạn có xác suất tìm thấy vàng chuyển từ bạc: có sang bạc: không.
Gregor Thomas

Một nhận xét cuối cùng: bạn đang sử dụng ma trận chuyển tiếp. Ma trận 2, 7, 14, 34 của bạn là ma trận chuyển tiếp.
Gregor Thomas
Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.