data(hiv)
### Initialize the proposals
mcmc.hiv<-fit.model(hiv[1:640,c(1:4)],hiv[1:640,c(5:8)],B=1000,min.iter=500,batch=1,shift=30,mcmc.obj=NULL,dye.swap=TRUE,nb.col1=2)
weight.plot(mcmc.hiv,hiv[1:640,9:10],array=3)
Run the code above in your browser using DataLab