วิธีการสร้างการทำนายด้วย rjags?


12

ฉันใช้ rjags เพื่อรัน MCMC ในรูปแบบที่ระบุในภาษา JAGS มีวิธีที่ดีในการแยกโมเดลนั้นและทำการทำนายด้วยหรือไม่ (ใช้การแจกแจงหลังของพารามิเตอร์ของฉัน)? ฉันสามารถระบุรุ่นใน R อีกครั้งและเสียบโหมดของพารามิเตอร์โปสเตอร์ของฉัน ฉันแค่สงสัยว่ามีวิธีการทำซ้ำซ้อนน้อยลงหรือไม่

ฉันเชื่อว่าhttp://sourceforge.net/p/mcmc-jags/discussion/610037/thread/0ecab41cกำลังถามคำถามเดียวกัน


ลิงก์เสีย กรุณาอัพเดทได้ไหม
chl

เสร็จสิ้น ดูเหมือนว่าแหล่งข้อมูลที่จัดเรียงใหม่
นักประวัติศาสตร์เชิงปริมาณ

ขอบคุณ! (ฉันไม่ได้พบด้ายตัวเองตั้งแต่ SF เปลี่ยนแปลงทางญาติของพวกเขา.)
CHL

คำตอบ:


8

โดยปกติคุณสามารถคาดการณ์ได้ใน JAGS ด้านล่างเป็นตัวอย่างการถดถอยของ FEV (บางอย่างเกี่ยวกับความจุปอด) เป็นตัวแปรตามและอายุและตัวบ่งชี้การสูบบุหรี่เป็นตัวทำนาย

FEV20s และ FEV20ns เป็นค่า FEV ที่คาดการณ์ไว้สำหรับผู้สูบบุหรี่อายุ 20 ปีและผู้ไม่สูบบุหรี่อายุ 20 ปี

model
{
for(i in 1:n){
    FEV[i] ~ dnorm(mu[i],tau)
    mu[i] <- beta[1] + beta[2]*Age[i] + beta[3]*Smoke[i]  + beta[4]*Age[i]*Smoke[i]
}

#priors
beta[1] ~ dnorm(0,0.001)
beta[2] ~ dnorm(0,0.001)
beta[3] ~ dnorm(0,0.001)
beta[4] ~ dnorm(0,0.001)
tau ~ dgamma(0.001,0.001)
sigma<-1/sqrt(tau) 

## Predict the FEV for a 20 year old smoker and for a 20 year old nonsmoker
mu20s <-  beta[1] + (beta[2]+beta[4])*20 + beta[3]
mu20ns <-  beta[1] + beta[2]*20 
FEV20s ~ dnorm(mu20s,tau)
FEV20ns ~ dnorm(mu20ns,tau)
}

ตัวอย่างจาก: แนวคิดแบบเบย์และการวิเคราะห์ข้อมูล


ขอบคุณสำหรับตัวชี้ - ฉันไม่เคยคิดที่จะส่งข้อมูลทดสอบของฉันไปยัง JAGS แต่ควรทำเช่นนั้น
นักประวัติศาสตร์เชิงปริมาณ

1
มีวิธีในการสร้างการทำนายเหล่านี้โดยไม่ต้องปรับโฉมทั้งโมเดลหรือไม่? หากมีมันจะง่ายพอที่จะทำให้การคาดการณ์ที่เกิดขึ้นขนานกันอย่างหนาแน่นอย่างไรก็ตามถ้าโมเดลทั้งหมดจำเป็นต้องได้รับการปรับโฉมสิ่งนี้ไม่สามารถทำได้
ลิน
โดยการใช้ไซต์ของเรา หมายความว่าคุณได้อ่านและทำความเข้าใจนโยบายคุกกี้และนโยบายความเป็นส่วนตัวของเราแล้ว
Licensed under cc by-sa 3.0 with attribution required.