# Read in callibration data
layout(matrix(c(1,2), 1, 1, byrow = TRUE))

opar<- par(mar=c(5, 3, 4, 3) + 0.1)

Z<- read.table("EXPOR000.TXT", header=TRUE, skip=0, nrows=1, sep=',')
iConc<-7:9
iVconc<-10:12
Conc<-as.real(Z[iConc])/1000
Vconc<-as.real(Z[iVconc])
Lconc<-log10(Conc)
ELconc<-predict(lm(Lconc~Vconc))
plot(Vconc, ELconc, ty='b', main="-pH vs mV CALLIBRATION")
points(Vconc, Lconc,col='red')
points(90, -5.25, col='red')
text(93, -5.25, "Obs[-pH]", col='red', pos=4)
points(90, -5.35, col='black')
text(93, -5.35, "E[-pH]", col='black', pos=4)
points(90, -5.45, col='blue', cex=2)
text(93, -5.45, "New E[-pH]", col='blue', pos=4)

NewVC <- data.frame(Vconc = c(100, 120, 140, 160, 180, 200))
points(t(NewVC), predict(lm(Lconc~Vconc), newdata=NewVC), cex=2, col='blue')
browser()

layout(matrix(c(1:8), 2, 4, byrow = TRUE))

filnam<-"EXPOR000.TXT"
X<-read.table(filnam, header=TRUE, skip=12, nrows=-1, sep=',')
print(names(X))
attach(X)
Geti<-1:100
Xpos<-as.real(ProbeX[Geti])-24990
Ypos<-as.real(ProbeY)
Zpos<-as.real(ProbeZ)
OmV<-as.real(Origin.1.mV[Geti])
XuV<-O.X.1.uV[Geti]
YuV<-O.Y.1.uV
ZuV<-O.Z.1.uV
detach(X)
plot(Xpos, OmV)
plot(Xpos, XuV)
NewVC <- data.frame(Vconc = OmV)
plot(Xpos, (OpH<- predict(lm(Lconc~Vconc), newdata=NewVC)), cex=2, col='blue')
plot(1/Xpos, (OpH<- predict(lm(Lconc~Vconc), newdata=NewVC)), cex=2, col='blue')