spec1<-read.table("pah1-ir.txt.out") names(spec1)<-c("1/cm","absorb") spec2<-read.table("pah2-ir.txt.out") names(spec2)<-c("1/cm","absorb") spec3<-read.table("pah3-ir.txt.out") names(spec3)<-c("1/cm","absorb") spec4<-read.table("fa1-ir.txt.out") names(spec4)<-c("1/cm","absorb") spec5<-read.table("fa2-ir.txt.out") names(spec5)<-c("1/cm","absorb") spec6<-read.table("fa3-ir.txt.out") names(spec6)<-c("1/cm","absorb") par(mfcol=c(2,1)) plot(spec1,type="l") lines(spec2, col="red") lines(spec3, col="blue") plot(spec4,type="l") lines(spec5, col="red") lines(spec6, col="blue")