d1 <- read.table("results_rmsip")
d2 <- read.table("hei_results")

pdf("rmsip.pdf", onefile=TRUE, width=6, height=3)
par(mar=c(3,2,1,1))
#postscript("methods_compare_overlap.ps", paper='special', width=6, height=6, horizontal=FALSE)
par(mfrow=c(1,1))

plot(d1[,1]/1000, d1[,2], type='o', xlab="Window size (ns)", ylab="RMSIP", main="", ylim=c(min(d2[,2]),1), xlim=c(0,200), cex=.5)
#lines(d2[,1]/1000, d2[,2], type='o', col=2, cex=.5)


#legend(12,0.7, c("C-alpha NMA", "All atom NMA", "PCA"), col=c(1,2,3), lty=c(1,1,1))



dev.off()
graphics.off()