for (SS in c(100,300,500)){ for (TT in 4:5){ for (i in 1:100){ file.name<-paste('data/data',i,'.txt', sep='') y<-read.table(file.name) y<-y[1:SS, 1:TT] file.name<-paste('data.txt', sep='') dset<-cbind(rep(1:TT, each=SS),rep(1:SS,TT),c(y)) write.table(dset,'ndata.txt',row.names=F,col.names=F) system(paste('sas tsas.sas', sep='')) out<-paste('output/norm-SS-',SS,'TT', TT,'-', i, '.txt', sep='') temp<-try(file.copy('histnorm.txt', out, overwrite=T)) out<-paste('output/t-SS-',SS,'TT', TT,'-', i, '.txt', sep='') temp<-try(file.copy('histt.txt', out, overwrite=T)) } } }