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