Displaying 1 result from an estimated 1 matches for "initials1".
Did you mean:
initials
2005 Jul 27
1
error message running R2WinBUGS
...,n[i])
#Hier.prior
p[i, j] ~ dbeta(alpha[i, j], beta[i, j])
alpha[i, j] < - pi[ j]*(1-theta[i])/theta[i]
beta[i, j] < -(1-pi[ j])*(1-theta[i])/theta[i]
}
}
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(
3, 47, 8, 19, 87, 69,
2, 4, 75, 24, 16, 81,
10, 78, 87, 44, 17, 56,
23, 75...