From: qi zhang <shellyzhang77_at_gmail.com>

Date: Thu 28 Jul 2005 - 01:21:09 EST

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Thu Jul 28 01:31:45 2005

Date: Thu 28 Jul 2005 - 01:21:09 EST

for(j in 1 : T ) {

pi[ j ] ~dbeta(0.1,0.1)I(0.1,0.9)

*}
*

for (i in 1:N) {

theta[i] ~dbeta(2,10)

*}
**}
*

And my R code is as followings:

initials1<-list(theta=c(0.2,0.01), pi=rbeta(50, 0.1, 0.1)*0.8+0.1

, p=matrix(rbeta(100, 2, 10)*0.8+0.1,nr=2,nc=50,byrow=TRUE))

inits<-list(initials1 initials1)

data<-list(N=2,T=50 ,n=c(100,150),x = structure(.Data = c(

structure<-bugs(data,inits,model.file="S:/Common/zhangqi/ebayes/for R/model.txt",

parameters=c("p","pi","theta"),n.chains=2,n.iter=2000,n.burnin=500,

debug=TRUE,DIC=FALSE,bugs.directory="c:/Program Files/WinBUGS14/",

check(S:/Common/zhangqi/ebayes/for R/model.txt)

model is syntactically correct

data(C:/Program Files/R/rw2011/data.txt)

inits(1,C:/Program Files/R/rw2011/inits1.txt)

value of binomial x[1,28] must be between zero and order of x[1,28]

I really have no clue about what I can do to fix it.

Could any of your please take a look at my problem, I truely appreciate any suggestions.

Qi Zhang

[[alternative HTML version deleted]]

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Thu Jul 28 01:31:45 2005

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:34:05 EST
*