
proc.coda<-function(stem,N, I){
  res<-NULL
  for (i in 1001:2000){
    codafile<-paste('../', stem,I,"/",stem,"CODAN",N,"I",i,"dic.txt",sep='')
    file.res<-paste(stem,'-', I,'-',N,'DIC.txt', sep='')
    if (file.exists(codafile)){
     DIC<-suppressWarnings(readLines(codafile))
     DIC<-DIC[substr(DIC,1,1)=="y"]
     
     cat(c(DIC, i), file=file.res, append=T)
     cat('\n', file=file.res, append=T)
   }  
  }

} 


proc.coda('Tdata',100,'Linear')
proc.coda('Tdata',200,'Linear')
proc.coda('Tdata',300,'Linear')
proc.coda('Tdata',400,'Linear')
proc.coda('Tdata',500,'Linear')

proc.coda('Tdata',100,'Latent')
proc.coda('Tdata',200,'Latent')
proc.coda('Tdata',300,'Latent')
proc.coda('Tdata',400,'Latent')
proc.coda('Tdata',500,'Latent')

proc.coda('Tdata',100,'Quad')
proc.coda('Tdata',200,'Quad')
proc.coda('Tdata',300,'Quad')
proc.coda('Tdata',400,'Quad')
proc.coda('Tdata',500,'Quad')

proc.coda('Tdata',100,'LinearA')
proc.coda('Tdata',200,'LinearA')
proc.coda('Tdata',300,'LinearA')
proc.coda('Tdata',400,'LinearA')
proc.coda('Tdata',500,'LinearA')

proc.coda('Tdata',100,'LatentA')
proc.coda('Tdata',200,'LatentA')
proc.coda('Tdata',300,'LatentA')
proc.coda('Tdata',400,'LatentA')
proc.coda('Tdata',500,'LatentA')

proc.coda('Tdata',100,'QuadA')
proc.coda('Tdata',200,'QuadA')
proc.coda('Tdata',300,'QuadA')
proc.coda('Tdata',400,'QuadA')
proc.coda('Tdata',500,'QuadA')