作者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