list(A=c(NA, .25, .5, .75, NA),k=4, bL=c(6), bS = c(1), Pre_e2 = 1, rho=0, Cov= structure(.Data = c(1,NA,NA,1),.Dim=c(2,2)))
list(A=c(NA, .25, .5, .75, NA),k=4, bL=c(0), bS = c(0), Pre_e2 = 1, Inv_cov= structure(.Data = c(1,0,0,1),.Dim=c(2,2)))
modelCheck("model.txt") modelData("data.txt") modelCompile() modelInits("inits.txt") modelGenInits() modelUpdate(10000) samplesSet("para") modelUpdate(50000) samplesCoda("para","coda") modelQuit()
tar czvf Ndata-Tdata-0.1.tgz NdataR