vv<-read.table("c:/vincentg/rfile.txt",header=TRUE); iter<-vv$iter; for (n in 0:49) { x<-vv$x[iter == n]; y<-vv$y[iter == n]; z<-vv$new[iter == n]; plot(x,y,xlim=c(0,1),ylim=c(0,1),pch=20,cex=0.75,col=1+z,xlab="",ylab="",axes=FALSE); dev.copy(png,filename=paste("c:/vincentg/test2_",n,".png",sep="")); dev.off (); }