Displaying 1 result from an estimated 1 matches for "codachain1".
2012 Jan 02
0
Reading mcmc/coda into a big.matrix efficiently
...ad.table("Big.CODAindex.txt", col.names =
c("var","start","end"))
n = index[1,3] - index[1,2] + 1
k = dim(index)[1]
X = filebacked.big.matrix( nrow = n, ncol = k, backingfile =
"Big.CODA.backing")
for(i in 1:k) { X[,i] = read.table("Big.CODAchain1.txt", skip =
(i-1)*n, nrows = n)[,2]
print(i)
print(Sys.time())
}
Also, here are the first few rows of the index and chain files, so you
can see the formatting. The index file tells you each variable's name
and the range or rows in the chain file containing the variab...