postscript("myFigure.ps", height=6, width=8, horizontal=T, onefile=F)
par(mar=c(4,4.5,3,1))
par(mfrow=c(1,1))


### Find global min/max values
dat0 <- read.table("../AUD_JPY.dat")
dat0$V1 <- dat0$V1
dat0$V2 <- dat0$V2
# ...
# All available currency pairs
# ...
dat13 <- read.table("../USD_JPY.dat")
dat13$V1 <- dat13$V1
dat13$V2 <- dat13$V2

minX <- min(dat0$V1, dat1$V1, dat2$V1, dat3$V1, dat4$V1, dat5$V1, dat6$V1, dat7$V1, dat8$V1, dat9$V1, dat10$V1, dat11$V1, dat12$V1, dat13$V1)
maxX <- 15*max(dat0$V1, dat1$V1, dat2$V1, dat3$V1, dat4$V1, dat5$V1, dat6$V1, dat7$V1, dat8$V1, dat9$V1, dat10$V1, dat11$V1, dat12$V1, dat13$V1)
minY <- min(dat0$V2, dat1$V2, dat2$V2, dat3$V2, dat4$V2, dat5$V2, dat6$V2, dat7$V2, dat8$V2, dat9$V2, dat10$V2, dat11$V2, dat12$V2, dat13$V2)
maxY <- max(dat0$V2, dat1$V2, dat2$V2, dat3$V2, dat4$V2, dat5$V2, dat6$V2, dat7$V2, dat8$V2, dat9$V2, dat10$V2, dat11$V2, dat12$V2, dat13$V2)


### Compute scaling law and generate plot
source("AUD_JPY.R")
par(new=TRUE)
# ...
# All available currency pairs
# ...
par(new=TRUE)
source("USD_JPY.R")
par(new=TRUE)
legend(maxX, maxY, legend=(c('AUD_JPY','AUD_USD','CHF_JPY','EUR_AUD','EUR_CHF','EUR_GBP','EUR_JPY','EUR_USD','GBP_CHF','GBP_JPY','GBP_USD','GRW_GRW','USD_CHF','USD_JPY')), col=c('yellow3','gray45','pink1','blue','green2','cyan2','orange','black','magenta','seagreen','gray65','red','brown4','salmon'), pch=c(16,21,17,24,18,23,15,22,12,25,10,8,9,11), lty=rep(0,14), cex=1.2, xjust=1, yjust=1)


dev.off()