
## Simulated t data
y<-read.table('data-generated-student-coda.txt')

ymat<-array(y[, 2], dim=c(1000,5, 500))

## write data into txt files

for (i in 1001:2000){
 t(ymat[i-1000,,])->temp
 filename<-paste('./Tdata/TdataN500I',i,'.txt', sep='')
 write.table(temp, filename, row.names=F, col.names=F)
}

source("http://nd.psychstat.org/_media/classdata/rtobugs.r")

for (i in 1001:2000){
 t(ymat[i-1000,,])->temp
 filename<-paste('./Tdata/TdataN500I',i,'bugs.txt', sep='')
 rtobugs(list(N=500, y=temp), filename)
}