作者lsshno1 (朝右边锋迈进)
看板Statistics
标题[程式] WinBug and R
时间Thu Nov 12 14:56:37 2009
大家好 我在编写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
1F：→ clickhere:困难? 11/12 17:54
2F：→ clickhere:parameters<-c("alpha","beta","sigma") 11/12 17:57
3F：→ lsshno1:我刚刚试过了还是错误耶 11/12 19:23
4F：→ clickhere:lambda[j] <- exp(alpha+beta * venue [j]+ error[j]) 11/13 09:09
5F：→ clickhere:错误讯息是? 11/13 09:10
6F：→ lsshno1:您好,我已经寄站内信给您了 11/13 17:36
7F：→ lsshno1:另外错误讯息为 invalid or unexpected token scanned 11/13 17:36
8F：→ clickhere:你要先确定你的bug code在winbugs中能跑! 11/14 02:40
9F：→ clickhere:然後在把它转到R中. 11/14 02:41
10F：→ clickhere:debug=TRUE, 会给你很多log file了 11/14 02:42
11F：→ clickhere:你有3个chain,但只给1个inits? 11/14 02:43
12F：→ clickhere:你有1个alpha, 却给J个rnorm? J个error, 却给1个rnorm? 11/14 02:44
13F：→ clickhere:tau的prior给太小了 winbug 不收 11/14 02:45