[程式] WinBug and R
大家好 我在編寫Winbug的程式時 出現了一些困難
請大家指點我一下 感激不盡
R 的code :
J<-38
y<-goal
venue<-data[,2]
data.m<-list("J","y","venue")
inits<-function(){
list(alpha=rnorm(J,0,100),beta=rnorm(1,0,100),error=rnorm(1,0,100))}
parameters<-c("Alpha","Beta","Sigma")
goal.sim <-
bugs(data.m,inits,parameters,"poissonmodel.txt",n.chain=3,n.iter=1000,debug=TRUE)
sim<-goal.sim$sims.list
Winbug poissmodel的code:
model {
for ( j in 1:J) {
y[j] ~ dpois (lambda[j])
error[j] ~ dnorm (0,tau)
log (lambda[j] ) <- ( alpha+beta * venue [j]+ error[j] )
}
alpha ~ dnorm (0.59168 , 0.05689)
beta ~ dnorm (0.23537 , 0.07612)
tau ~ dgamma (0.001,0.001)
sigma <- 1 / sqrt(tau)
}
註:我的模型為 log (lamdba) = alpha + beta * venue + error
感謝大家!!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.129.77.1
→
11/12 17:54, , 1F
11/12 17:54, 1F
→
11/12 17:57, , 2F
11/12 17:57, 2F
→
11/12 19:23, , 3F
11/12 19:23, 3F
→
11/13 09:09, , 4F
11/13 09:09, 4F
→
11/13 09:10, , 5F
11/13 09:10, 5F
→
11/13 17:36, , 6F
11/13 17:36, 6F
→
11/13 17:36, , 7F
11/13 17:36, 7F
→
11/14 02:40, , 8F
11/14 02:40, 8F
→
11/14 02:41, , 9F
11/14 02:41, 9F
→
11/14 02:42, , 10F
11/14 02:42, 10F
→
11/14 02:43, , 11F
11/14 02:43, 11F
→
11/14 02:44, , 12F
11/14 02:44, 12F
→
11/14 02:45, , 13F
11/14 02:45, 13F
→
01/02 15:00,
5年前
, 14F
01/02 15:00, 14F