Tôi đang cố gắng xây dựng một mô hình trong đó phản hồi là một tỷ lệ (thực ra đó là tỷ lệ phiếu bầu mà một đảng nhận được trong các khu vực bầu cử). Phân phối của nó không bình thường, vì vậy tôi quyết định mô hình hóa nó với phân phối beta. Tôi cũng có một vài dự đoán.
Tuy nhiên, tôi không biết cách viết nó bằng BUGS / JAGS / STAN (JAGS sẽ là lựa chọn tốt nhất của tôi, nhưng nó không thực sự quan trọng). Vấn đề của tôi là tôi tạo ra một tổng các tham số bằng các yếu tố dự đoán, nhưng sau đó tôi có thể làm gì với nó?
Mã sẽ giống như thế này (theo cú pháp JAGS), nhưng tôi không biết cách "liên kết" các tham số y_hat
và y
.
for (i in 1:n) {
y[i] ~ dbeta(alpha, beta)
y_hat[i] <- a + b * x[i]
}
( y_hat
chỉ là sản phẩm chéo của các tham số và yếu tố dự đoán, do đó có mối quan hệ xác định. a
và b
là các hệ số mà tôi cố gắng ước tính, x
là một yếu tố dự báo).
Cảm ơn lời đề nghị của bạn!