मैं एक मॉडल बनाने की कोशिश कर रहा हूं, जहां प्रतिक्रिया एक अनुपात है (यह वास्तव में एक पार्टी को निर्वाचन क्षेत्रों में मिले वोटों का हिस्सा है)। इसका वितरण सामान्य नहीं है, इसलिए मैंने इसे बीटा वितरण के साथ मॉडल करने का निर्णय लिया। मेरे पास कई भविष्यवक्ता भी हैं।
हालाँकि, मुझे नहीं पता कि इसे BUGS / JAGS / STAN में कैसे लिखा जाता है (JAGS मेरी सबसे अच्छी पसंद होगी, लेकिन यह वास्तव में मायने नहीं रखता है)। मेरी समस्या यह है कि मैं भविष्यवक्ताओं द्वारा मापदंडों का योग बनाता हूं, लेकिन फिर मैं इसके साथ क्या कर सकता हूं?
कोड कुछ इस तरह होगा (JAGS सिंटैक्स में), लेकिन मुझे पता है कि पैरामीटर y_hat
और y
मापदंडों को "लिंक" कैसे किया जाता है ।
for (i in 1:n) {
y[i] ~ dbeta(alpha, beta)
y_hat[i] <- a + b * x[i]
}
( y_hat
सिर्फ मापदंडों और भविष्यवाणियों का क्रॉस-उत्पाद है, इसलिए निर्धारक संबंध हैं। a
और b
वे गुणांक हैं जो मैं अनुमान लगाने की कोशिश करता हूं, x
एक भविष्यवक्ता होने के नाते)।
आपके सुझाव के लिए धन्यवाद!