Displaying 1 result from an estimated 1 matches for "aasho1".
Did you mean:
aasho
2008 Dec 15
3
R2winbugs : vectorization
...y=y,z=z,n=n,C=C)
inits1<-function(){
list(a1=runif(1),a2=runif(1),a3=runif(1),a4=runif(1),sigma.y=runif(1))}
params1<-c("a1","a2","a3","a4","sigma.y")
## now the bugs call for model 1
proc.time()
aasho.1<-bugs(data1,inits1,params1,"aasho1.bug",n.iter=10000,debug=FALSE)
proc.time()
print(aasho.1,digits=4)
now we try a vector approach
k<-4 # niv
data2<-list(x=x,y=y,z=z,n=n,C=C,k=k)
inits2<-function(){
list(a=runif(k),sigma.y=runif(1))}
params2<-c("a","sigma.y")
## now the bugs call for model 2...