search for: st90

...#number of days Y <- subsetA$concentr #concentration per day/subject t <- 1:13 #days dataA <- list("n","nt","Y","t") parameters <- c("tau","C0","st90","C0.pop","st90.pop","tau.cpop","tau.stpop") inits <- list(tau=c(1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1,1),C0=5,st90=4,C0.pop=5,st90.pop=4,tau.cpop=0.2,tau.stpop=1) mcmcA <- bugs(dataA,inits,parameters,modelA,debug=T,n.chains=7,bugs.directory=&q...