Nói rằng tôi có một thí nghiệm trong đó tôi kiểm tra thời gian phản ứng của một số đối tượng trong đó mỗi đối tượng thực hiện nhiều thử nghiệm thời gian phản ứng. Trong khung Bayes, thời gian phản ứng ( ) có thể được mô hình hóa bằng mô hình phân cấp với phân phối trước cả ở cấp độ chủ đề và cho toàn bộ nhóm đối tượng. Một sơ đồ của mô hình, kiểu Kruschke , có thể là:
... và mã BUGS / JAGS tương ứng sẽ là:
for(i in 1:length(y)) {
y[i] ~ dnorm(mu[subj[i]], tau[subj[i]])
}
for(j in 1:nbr_of_subjects)
mu[subj[i]] ~ dnorm(M_mu, P_mu)
tau[subj[i]] ~ dgamma(S_tau, R_tau)
}
M_mu ~ dnorm(M_M, P_M)
P_mu ~ dgamma(S_P, R_P)
S_tau <- pow(m , 2) / pow(sd, 2)
R_tau <- m / pow(sd, 2)
m ~ dgamma(S_m, R_m)
sd ~ dgamma(S_sd, R_sd)
Nếu tôi muốn so sánh thời gian phản ứng của hai đối tượng thì tôi sẽ so sánh phân phối tương ứng của họ . Nếu các thử nghiệm thời gian phản ứng được chia thành bốn khối, tôi cũng có thể mô hình hóa bằng cách thêm một cấp độ khối bổ sung với các mục sư giữa cấp độ đối tượng và cấp độ thử nghiệm trong sơ đồ (vì có thể trường hợp thời gian phản ứng của các đối tượng khác nhau một chút đối với một số lý do).
Câu hỏi của tôi là bây giờ, nếu tôi muốn so sánh hai đối tượng thì tôi nên so sánh phân phối nào? Tôi có thể so sánh phân phối của phương tiện ở cấp độ chủ đề (hiện đã xác định một phần mức trung bình trước mức trung bình ở cấp độ khối) nhưng tôi cũng có thể so sánh phân phối phương tiện ở cấp độ khối tương ứng với trong mô hình cũ . Theo một cách nào đó có vẻ hợp lý hơn khi so sánh các đối tượng ở cấp độ chủ đề, nhưng nó có tạo ra sự khác biệt nào không? Và nếu có rất ít khối, giả sử hai khối, việc phân phối phương tiện ở cấp độ đối tượng sẽ rất "rộng" phải không?