# Lobster_Lesion_CaH_080613.R
filnam="SD04SmExuLesionGrid.csv"
ZZ<- read.csv(filnam, skip=25)
attach(ZZ)
print(names(ZZ))
Min_From<-ZZ[,1]
X_microns<-ZZ[,2]
Y_microns<-ZZ[,3]
Z_microns<-ZZ[,4]

quartz(width=18, height=6)
plot(c(min(Min_From),max(Min_From)), c(min(Y_microns),max(Y_microns)), 
     typ='n', 
     main=paste(filnam,"Position, X=red, Y=green, Z=blue"),
     xlab="time (min)",
     ylab="position um")
points(Min_From, Z_microns, col='blue', pch=20)
points(Min_From, Y_microns, col='green', pch=20)
points(Min_From, X_microns, col='red', pch=20)

browser()
quartz(width=18, height=6)
plot(X_microns, Y_microns, main=paste(filnam,"Y_microns vs X_microns")) 
browser()
quartz(width=18, height=6)

## AvggradHuV,AvggradCauV
batch1<- 1:825 
plot(c(min(Min_From), max(Min_From)), c(-150,400), typ='n', 
      main=paste(filnam,"uV vs Time, H(black), Ca(red), batch-1")) 
points(Min_From[batch1], AvggradHuV[batch1], col='black', typ='l')
points(Min_From[batch1], AvggradCauV[batch1], col='red', typ='l')

points(Min_From[batch1], gradH_3_uV[batch1], col='black', typ='l', lwd=3)
points(Min_From[batch1], gradCa_3_uV[batch1], col='red', typ='l', lwd=3)
for (i in seq(1,825,25)) {lines(c(Min_From[i],Min_From[i]), c(-150,400), col='green')}

browser()
## AvggradHuV,AvggradCauV
quartz(width=8, height=6)
batch1<- 1:825 
plot(c(min(X_microns[batch1]), max(X_microns[batch1])), c(-150,400), typ='n', 
      main=paste(filnam,"uV vs X_microns, H(black), Ca(red), batch-1")) 
points(X_microns[batch1], AvggradHuV[batch1], col='black', typ='l')
points(X_microns[batch1], AvggradCauV[batch1], col='red', typ='l')

browser()

detach(ZZ)