require(session)
rainboc<- c("red", "yellow", "orange", "green", "blue", "purple", "pink", "brown", "coral")
junk<-read.table("Sym_v109mf-alig.ang", header = TRUE, row.names = 1)
attach(junk)
Nvar<- length(NamVar<-names(junk))
range<- "[1:109]"
range1<- "[1:35]"
range2<- "[36:109]"

N<-length(junk[,1])
All<-matrix(0, nrow = 0.5*Nvar*N, ncol = 2)
for (i in 1:(Nvar/2)) {
    j<- 1+ 2*(i-1)
    All[(i-1)*N+(1:N),]<- as.matrix(junk[,j:(j+1)])
                        }

     op <- par(mfrow = c(1, 2), # 2 x 3 pictures on one plot
               pty = "m")       # square plotting region,

for (i in seq(1, Nvar, 2)) {

plot(All) 
point<- paste("points(", NamVar[i],range2,",", NamVar[i+1], range2,", col='blue')", sep='')
texteval(point)
point<- paste("points(", NamVar[i],range1,",", NamVar[i+1], range1,", col='red')", sep='')
texteval(point)
                       
line<- paste("plot(", NamVar[i],range1,",", NamVar[i+1], range1,", col='red')", sep='')
texteval(line) 
point<- paste("points(", NamVar[i],range2,",", NamVar[i+1], range2,", col='blue')", sep='')
texteval(point)
browser()                   }