### Data

dat <- read.table("../AUD_JPY.dat")
lab <- 'AUD_JPY'
xlabT <- 'XLAB'
ylabT <- 'YLAB'


### Estimate

# Approx. percentage
dat$V1 <- dat$V1
dat$V2 <- dat$V2

ldat <- log(dat)

Min <- min(ldat$V1)
Max <- max(ldat$V1)

X <- ldat$V1[ldat$V1 >= Min & ldat$V1 <= Max]
Y <- ldat$V2[ldat$V1 >= Min & ldat$V1 <= Max]

res <- lm(Y~X)

# Error propagation
A <- summary(res)$coefficients[1,1]
dA <- summary(res)$coefficients[1,2]
B <- summary(res)$coefficients[2,1]
dB <- summary(res)$coefficients[2,2]
C <- exp(-A/B)
C.err <- sqrt((-C/B*dA)^2+(A/B^2*C*dB)^2)

slEmpRSx <- function(t) (t/C)^B


### Plot
start <- 0
end <- 68
pchi <- c(rep(NA,start), 16, rep(NA,end));

plot(exp(X),exp(Y), xlim=c(minX, maxX), log='xy', ylim=c(minY, maxY), cex=1.5, lwd=2.5,  col='yellow3', ylab = ylabT, xlab = xlabT, main='TITLE', pch=pchi, type='l', lty=3, bg='white', yaxt='n', xaxt='n', cex.main=2, cex.lab=1.5, mgp=c(3.2,0,0))

par(new=TRUE)

plot(exp(X),exp(Y), xlim=c(minX, maxX), log='xy', ylim=c(minY, maxY), cex=1.5, lwd=1.5,  col='yellow3', ylab = ylabT, xlab = xlabT, main='TITLE', pch=pchi, type='p', lty=0, bg='white', yaxt='n', xaxt='n', cex.main=2, cex.lab=1.5, mgp=c(3.2,0,0))

cexaxis <- 1

minv  <- round((minY),digits=2)
sequ <- 10^((ceiling(log10(minY)):ceiling(log10(maxY))))
axis(2, at =  format(sequ, digits=2,sci = T), cex.axis=cexaxis)

axis.at <- 10 ^ c(floor(log10(minY)):ceiling(log10(maxY)))
sq <- 1:10 * rep(axis.at[-1] / 10, each = 10)
axis(2, at = sq[-(1:3)] , tcl = -0.3, labels = FALSE, cex.axis=cexaxis) 

sequ <- 10^((floor(log10(minX)):ceiling(log10(maxX))))
axis(1, at =  format(sequ, digits=2,sci = TRUE), cex.axis=cexaxis)

axis.at <- 10 ^ c(floor(log10(minX)):ceiling(log10(maxX)))
axis(1, at = 1:10 * rep(axis.at[-1] / 10, each = 10),tcl = -0.3, labels = FALSE, cex.axis=cexaxis) 

title(main='LAB')