search for: rbern

...SeaW") #,"Month" #Month Taken out so that model runs #create initial values for the coefficients in the model (this is model specific) #construct inits for the three chain inits1 <- list(gamma0 = 0, delta0 = 0,beta1 = 0, alpha1 = 0, alpha2 = 0, index = rbern(J,0.5)+1) inits2 <- list(gamma0 = 0, delta0 = 0,beta1 = 0, alpha1 = 0, alpha2 = 0, index = rbern(J,0.5)+1) inits3 <- list(gamma0 = 0, delta0 = 0,beta1 = 0, alpha1 = 0, alpha2 = 0, index = rbern(J,0.5)+1) inits <- list(inits1,inits2,inits3) ## these data.frames need to be trimmed so...