rm(list=ls()) quartz(width=18, height=6) layout(matrix(c(1,1,1,1,1,2,2,2,2,2), 2, 5, byrow = TRUE)) filnam<- "jgk080724_031.csv" kline<-290 id<-1.1 XX<-read.csv(filnam, skip=145, nrows=67) # 1533 attach(XX) print(names(XX)) ni<-length(X_microns) minHr<-min(Hrs_From) maxHr<-max(Hrs_From) minX<-min(X_microns) maxX<-max(X_microns) minYloc<-min(Y_microns) maxYloc<-max(Y_microns) minZloc<-min(Z_microns) maxZloc<-max(Z_microns) leng<- length(X_microns) medCauV<-medHuV<- rep(0,leng) for (i in 1:leng) {medHuV[i]<- median(c(gradH_8_uV[i],gradH_9_uV[i],gradH_10_uV[i])) medCauV[i]<- median(c(gradCa_8_uV[i],gradCa_9_uV[i],gradCa_10_uV[i]))} minY<- -500; maxY<- 500 plot(c(minHr, maxHr),c(minY, maxY),typ='n', main=paste(filnam,"H/Ca uV blue/purple median red/green #10", id)) # lines(Hrs_From,AvggradHuV, typ='b', col='blue', cex=2, lwd=2) # lines(Hrs_From, AvggradCauV, typ='b', col='purple', cex=2, lwd=2) lines(Hrs_From,medHuV, typ='b', col='blue', cex=2, lwd=2) lines(Hrs_From, medCauV, typ='b', col='purple', cex=2, lwd=2) lines(Hrs_From,gradH_10_uV, typ='b', col='red', cex=2, lwd=2) lines(Hrs_From, gradCa_10_uV, typ='b', col='green', cex=2, lwd=2) browser() minY<- 87; maxY<- 97 plot(c(minHr, maxHr),c(minY, maxY),typ='n', main=paste(filnam,"H/Ca mV blue/purple 0/Zavg red/green 0", id)) lines(Hrs_From,(NagelSamp0H_10_mV+NagelSampZH_10_mV)/2, typ='b', col='red', cex=2, lwd=2) lines(Hrs_From,100*DA0_V_, typ='l', col='black', cex=2, lwd=2, lty=2) par(new=TRUE, ann=FALSE) minY<- -5; maxY<- 18 plot(c(minHr, maxHr),c(minY, maxY),typ='n', yaxt='n') lines(Hrs_From,(NagelRead0Ca_10_mV+NagelReadZCa_10_mV)/2, typ='b', col='green', cex=2, lwd=2) lines(Hrs_From,100*DA1_V_, typ='l', col='black', cex=2, lwd=2, lty=3) par(ann=TRUE) browser() quartz(width=8, height=8) layout(matrix(1, 1, 1, byrow = TRUE)) minY<- 87; maxY<- 97 plot(c(minZloc, maxZloc),c(minY, maxY),typ='n', main=paste(filnam,"lesion", id)) lines(Z_microns,(NagelSamp0H_10_mV+NagelSampZH_10_mV)/2, typ='b', col='red', cex=2, lwd=2) lines(Z_microns,100*DA0_V_, typ='l', col='black', cex=2, lwd=2, lty=2) par(new=TRUE, ann=FALSE) minY<- -5; maxY<- 18 plot(c(minZloc, maxZloc),c(minY, maxY),typ='n', yaxt='n') lines(Z_microns, (NagelRead0Ca_10_mV+NagelReadZCa_10_mV)/2, typ='b', col='green', cex=2, lwd=2) lines(Z_microns,100*DA1_V_, typ='l', col='black', cex=2, lwd=2, lty=3) browser() par(ann=TRUE) plot(c(minX, maxX),c(minZloc, maxZloc),typ='n', main=paste(filnam,"XY location")) lines(X_microns,Z_microns, typ='b', col='blue', cex=2, lwd=2) browser() minY<- -500; maxY<- 500 plot(c(minZloc, maxZloc),c(minY, maxY),typ='n', main=paste(filnam,"uV vs Zloc")) lines(Z_microns[X_microns < -200],medHuV[X_microns < -200], typ='p', col='blue', cex=2, lwd=2) lines(Z_microns[X_microns < -200],medCauV[X_microns < -200], typ='p', col='purple', cex=2, lwd=2) browser() minY<- -500; maxY<- 500 plot(c(minHr, maxHr),c(minY, maxY),typ='n', main=paste(filnam,"uV vs Zloc")) lines(Hrs_From[Z_microns == 0.1],medHuV[Z_microns == 0.1], typ='b', col='blue', cex=2, lwd=4) lines(Hrs_From[Z_microns == 0.1],medCauV[Z_microns == 0.1], typ='b', col='purple', cex=2, lwd=4) lines(Hrs_From[Z_microns == 200.1],medHuV[Z_microns == 200.1], typ='b', col='blue', cex=2, lwd=3) lines(Hrs_From[Z_microns == 200.1],medCauV[Z_microns == 200.1], typ='b', col='purple', cex=2, lwd=3) lines(Hrs_From[Z_microns == 400.1],medHuV[Z_microns == 400.1], typ='b', col='blue', cex=2, lwd=2) lines(Hrs_From[Z_microns == 400.1],medCauV[Z_microns == 400.1], typ='b', col='purple', cex=2, lwd=2) browser() minY<- 87; maxY<- 97 HmV<- (NagelSamp0H_10_mV+NagelSampZH_10_mV)/2 CamV<- (NagelRead0Ca_10_mV+NagelReadZCa_10_mV)/2 plot(c(minHr, maxHr),c(minY, maxY),typ='n', main=paste(filnam,"mV vs Time at pH 7.88")) lines(Hrs_From[Z_microns == 0.1],HmV[Z_microns == 0.1], typ='b', col='red', cex=2, lwd=4, lty=2) lines(Hrs_From[Z_microns == 200.1],HmV[Z_microns == 200.1], typ='b', col='green', cex=2, lwd=3, lty=3) lines(Hrs_From[Z_microns == 400.1],HmV[Z_microns == 400.1], typ='b', col='blue', cex=2, lwd=2, lty=4) lines(Hrs_From[Z_microns == 600.1],HmV[Z_microns == 600.1], typ='b', col='orange', cex=2, lwd=2, lty=5) lines(Hrs_From[Z_microns == 800.1],HmV[Z_microns == 800.1], typ='b', col='purple', cex=2, lwd=2, lty=6) par(new=TRUE, ann=FALSE) minY<- -5; maxY<- 18 plot(c(minHr, maxHr),c(minY, maxY),typ='n', yaxt='n') lines(Hrs_From[Z_microns == 0.1],CamV[Z_microns == 0.1], typ='b', col='black', cex=2, lwd=4, lty=2) lines(Hrs_From[Z_microns == 200.1],CamV[Z_microns == 200.1], typ='b', col='brown', cex=2, lwd=3, lty=3) lines(Hrs_From[Z_microns == 400.1],CamV[Z_microns == 400.1], typ='b', col='grey', cex=2, lwd=2, lty=4) lines(Hrs_From[Z_microns == 600.1],CamV[Z_microns == 600.1], typ='b', col='purple', cex=2, lwd=1, lty=5) lines(Hrs_From[Z_microns == 800.1],CamV[Z_microns == 800.1], typ='b', col='purple', cex=2, lwd=1, lty=6) browser() detach(XX) rm(list=ls()) quartz(width=18, height=6) layout(matrix(c(1,1,1,1,1,2,2,2,2,2), 2, 5, byrow = TRUE)) filnam<- "jgk080724_031.csv" kline<-290 id<-1.1 XX<-read.csv(filnam, skip=306, nrows=-1) # 1533 attach(XX) print(names(XX)) ni<-length(X_microns) minHr<-min(Hrs_From) maxHr<-max(Hrs_From) minX<-min(X_microns) maxX<-max(X_microns) minYloc<-min(Y_microns) maxYloc<-max(Y_microns) minZloc<-min(Z_microns) maxZloc<-max(Z_microns) leng<- length(X_microns) medCauV<-medHuV<- rep(0,leng) for (i in 1:leng) {medHuV[i]<- median(c(gradH_8_uV[i],gradH_9_uV[i],gradH_10_uV[i])) medCauV[i]<- median(c(gradCa_8_uV[i],gradCa_9_uV[i],gradCa_10_uV[i]))} minY<- -500; maxY<- 500 plot(c(minHr, maxHr),c(minY, maxY),typ='n', main=paste(filnam,"H/Ca uV blue/purple median red/green #10", id)) # lines(Hrs_From,AvggradHuV, typ='b', col='blue', cex=2, lwd=2) # lines(Hrs_From, AvggradCauV, typ='b', col='purple', cex=2, lwd=2) lines(Hrs_From,medHuV, typ='b', col='blue', cex=2, lwd=2) lines(Hrs_From, medCauV, typ='b', col='purple', cex=2, lwd=2) lines(Hrs_From,gradH_10_uV, typ='b', col='red', cex=2, lwd=2) lines(Hrs_From, gradCa_10_uV, typ='b', col='green', cex=2, lwd=2) browser() minY<- 95.5; maxY<- 97.5 plot(c(minHr, maxHr),c(minY, maxY),typ='n', main=paste(filnam,"H/Ca mV blue/purple 0/Zavg red/green 0", id)) lines(Hrs_From,(NagelSamp0H_10_mV+NagelSampZH_10_mV)/2, typ='b', col='red', cex=2, lwd=2) lines(Hrs_From,100*DA0_V_, typ='l', col='black', cex=2, lwd=2, lty=2) par(new=TRUE, ann=FALSE) minY<- -5; maxY<- 18 plot(c(minHr, maxHr),c(minY, maxY),typ='n', yaxt='n') lines(Hrs_From,(NagelRead0Ca_10_mV+NagelReadZCa_10_mV)/2, typ='b', col='green', cex=2, lwd=2) lines(Hrs_From,100*DA1_V_, typ='l', col='black', cex=2, lwd=2, lty=3) par(ann=TRUE) browser() quartz(width=8, height=8) layout(matrix(1, 1, 1, byrow = TRUE)) minY<- 95.5; maxY<- 97.5 plot(c(minZloc, maxZloc),c(minY, maxY),typ='n', main=paste(filnam,"lesion", id)) lines(Z_microns,(NagelSamp0H_10_mV+NagelSampZH_10_mV)/2, typ='b', col='red', cex=2, lwd=2) lines(Z_microns,100*DA0_V_, typ='l', col='black', cex=2, lwd=2, lty=2) par(new=TRUE, ann=FALSE) minY<- -5; maxY<- 18 plot(c(minZloc, maxZloc),c(minY, maxY),typ='n', yaxt='n') lines(Z_microns, (NagelRead0Ca_10_mV+NagelReadZCa_10_mV)/2, typ='b', col='green', cex=2, lwd=2) lines(Z_microns,100*DA1_V_, typ='l', col='black', cex=2, lwd=2, lty=3) browser() par(ann=TRUE) plot(c(minX, maxX),c(minZloc, maxZloc),typ='n', main=paste(filnam,"XY location")) lines(X_microns,Z_microns, typ='b', col='blue', cex=2, lwd=2) browser() minY<- -500; maxY<- 500 plot(c(minZloc, maxZloc),c(minY, maxY),typ='n', main=paste(filnam,"uV vs Zloc")) lines(Z_microns[X_microns < -200],medHuV[X_microns < -200], typ='p', col='blue', cex=2, lwd=2) lines(Z_microns[X_microns < -200],medCauV[X_microns < -200], typ='p', col='purple', cex=2, lwd=2) browser() minY<- -500; maxY<- 500 plot(c(minHr, maxHr),c(minY, maxY),typ='n', main=paste(filnam,"uV vs Zloc")) lines(Hrs_From[Z_microns == 0],medHuV[Z_microns == 0], typ='b', col='blue', cex=2, lwd=4) lines(Hrs_From[Z_microns == 0],medCauV[Z_microns == 0], typ='b', col='purple', cex=2, lwd=4) lines(Hrs_From[Z_microns == 200],medHuV[Z_microns == 200], typ='b', col='blue', cex=2, lwd=3) lines(Hrs_From[Z_microns == 200],medCauV[Z_microns == 200], typ='b', col='purple', cex=2, lwd=3) lines(Hrs_From[Z_microns == 400],medHuV[Z_microns == 400], typ='b', col='blue', cex=2, lwd=2) lines(Hrs_From[Z_microns == 400],medCauV[Z_microns == 400], typ='b', col='purple', cex=2, lwd=2) browser() minY<- 95.5; maxY<- 97.5 HmV<- (NagelSamp0H_10_mV+NagelSampZH_10_mV)/2 CamV<- (NagelRead0Ca_10_mV+NagelReadZCa_10_mV)/2 plot(c(minHr, maxHr),c(minY, maxY),typ='n', main=paste(filnam,"mV vs Time at pH 7.36")) lines(Hrs_From[Z_microns == 0],HmV[Z_microns == 0], typ='b', col='red', cex=2, lwd=4, lty=2) lines(Hrs_From[Z_microns == 200],HmV[Z_microns == 200], typ='b', col='green', cex=2, lwd=3, lty=3) lines(Hrs_From[Z_microns == 400],HmV[Z_microns == 400], typ='b', col='blue', cex=2, lwd=2, lty=4) lines(Hrs_From[Z_microns == 600],HmV[Z_microns == 600], typ='b', col='orange', cex=2, lwd=2, lty=5) lines(Hrs_From[Z_microns == 800],HmV[Z_microns == 800], typ='b', col='purple', cex=2, lwd=2, lty=6) par(new=TRUE, ann=FALSE) minY<- -5; maxY<- 18 plot(c(minHr, maxHr),c(minY, maxY),typ='n', yaxt='n') lines(Hrs_From[Z_microns == 0],CamV[Z_microns == 0], typ='b', col='black', cex=2, lwd=4, lty=2) lines(Hrs_From[Z_microns == 200],CamV[Z_microns == 200], typ='b', col='brown', cex=2, lwd=3, lty=3) lines(Hrs_From[Z_microns == 400],CamV[Z_microns == 400], typ='b', col='grey', cex=2, lwd=2, lty=4) lines(Hrs_From[Z_microns == 600],CamV[Z_microns == 600], typ='b', col='purple', cex=2, lwd=1, lty=5) lines(Hrs_From[Z_microns == 800],CamV[Z_microns == 800], typ='b', col='purple', cex=2, lwd=1, lty=6) browser() detach(XX) rm(list=ls()) quartz(width=18, height=6) layout(matrix(c(1,1,1,1,1,2,2,2,2,2), 2, 5, byrow = TRUE)) filnam<- "jgk080724_031.csv" kline<-290 id<-1.1 XX<-read.csv(filnam, skip=226, nrows=66) # 1533 attach(XX) print(names(XX)) ni<-length(X_microns) minHr<-min(Hrs_From) maxHr<-max(Hrs_From) minX<-min(X_microns) maxX<-max(X_microns) minYloc<-min(Y_microns) maxYloc<-max(Y_microns) minZloc<-min(Z_microns) maxZloc<-max(Z_microns) leng<- length(X_microns) medCauV<-medHuV<- rep(0,leng) for (i in 1:leng) {medHuV[i]<- median(c(gradH_8_uV[i],gradH_9_uV[i],gradH_10_uV[i])) medCauV[i]<- median(c(gradCa_8_uV[i],gradCa_9_uV[i],gradCa_10_uV[i]))} minY<- -500; maxY<- 500 plot(c(minHr, maxHr),c(minY, maxY),typ='n', main=paste(filnam,"H/Ca uV blue/purple median red/green #10", id)) # lines(Hrs_From,AvggradHuV, typ='b', col='blue', cex=2, lwd=2) # lines(Hrs_From, AvggradCauV, typ='b', col='purple', cex=2, lwd=2) lines(Hrs_From,medHuV, typ='b', col='blue', cex=2, lwd=2) lines(Hrs_From, medCauV, typ='b', col='purple', cex=2, lwd=2) lines(Hrs_From,gradH_10_uV, typ='b', col='red', cex=2, lwd=2) lines(Hrs_From, gradCa_10_uV, typ='b', col='green', cex=2, lwd=2) browser() minY<- 122; maxY<- 131 plot(c(minHr, maxHr),c(minY, maxY),typ='n', main=paste(filnam,"H/Ca mV blue/purple 0/Zavg red/green 0", id)) lines(Hrs_From,(NagelSamp0H_10_mV+NagelSampZH_10_mV)/2, typ='b', col='red', cex=2, lwd=2) lines(Hrs_From,100*DA0_V_, typ='l', col='black', cex=2, lwd=2, lty=2) par(new=TRUE, ann=FALSE) minY<- -5; maxY<- 18 plot(c(minHr, maxHr),c(minY, maxY),typ='n', yaxt='n') lines(Hrs_From,(NagelRead0Ca_10_mV+NagelReadZCa_10_mV)/2, typ='b', col='green', cex=2, lwd=2) lines(Hrs_From,100*DA1_V_, typ='l', col='black', cex=2, lwd=2, lty=3) par(ann=TRUE) browser() quartz(width=8, height=8) layout(matrix(1, 1, 1, byrow = TRUE)) minY<- 122; maxY<- 131 plot(c(minZloc, maxZloc),c(minY, maxY),typ='n', main=paste(filnam,"lesion", id)) lines(Z_microns,(NagelSamp0H_10_mV+NagelSampZH_10_mV)/2, typ='b', col='red', cex=2, lwd=2) lines(Z_microns,100*DA0_V_, typ='l', col='black', cex=2, lwd=2, lty=2) par(new=TRUE, ann=FALSE) minY<- -5; maxY<- 18 plot(c(minZloc, maxZloc),c(minY, maxY),typ='n', yaxt='n') lines(Z_microns, (NagelRead0Ca_10_mV+NagelReadZCa_10_mV)/2, typ='b', col='green', cex=2, lwd=2) lines(Z_microns,100*DA1_V_, typ='l', col='black', cex=2, lwd=2, lty=3) browser() par(ann=TRUE) plot(c(minX, maxX),c(minZloc, maxZloc),typ='n', main=paste(filnam,"XY location")) lines(X_microns,Z_microns, typ='b', col='blue', cex=2, lwd=2) browser() minY<- -500; maxY<- 500 plot(c(minZloc, maxZloc),c(minY, maxY),typ='n', main=paste(filnam,"uV vs Zloc")) lines(Z_microns[X_microns < -200],medHuV[X_microns < -200], typ='p', col='blue', cex=2, lwd=2) lines(Z_microns[X_microns < -200],medCauV[X_microns < -200], typ='p', col='purple', cex=2, lwd=2) browser() minY<- -500; maxY<- 500 plot(c(minHr, maxHr),c(minY, maxY),typ='n', main=paste(filnam,"uV vs Zloc")) lines(Hrs_From[Z_microns == 0],medHuV[Z_microns == 0], typ='b', col='blue', cex=2, lwd=4) lines(Hrs_From[Z_microns == 0],medCauV[Z_microns == 0], typ='b', col='purple', cex=2, lwd=4) lines(Hrs_From[Z_microns == 200],medHuV[Z_microns == 200], typ='b', col='blue', cex=2, lwd=3) lines(Hrs_From[Z_microns == 200],medCauV[Z_microns == 200], typ='b', col='purple', cex=2, lwd=3) lines(Hrs_From[Z_microns == 400],medHuV[Z_microns == 400], typ='b', col='blue', cex=2, lwd=2) lines(Hrs_From[Z_microns == 400],medCauV[Z_microns == 400], typ='b', col='purple', cex=2, lwd=2) browser() minY<- 122; maxY<- 131 HmV<- (NagelSamp0H_10_mV+NagelSampZH_10_mV)/2 CamV<- (NagelRead0Ca_10_mV+NagelReadZCa_10_mV)/2 plot(c(minHr, maxHr),c(minY, maxY),typ='n', main=paste(filnam,"mV vs Time at pH 6.72")) lines(Hrs_From[Z_microns == 0],HmV[Z_microns == 0], typ='b', col='red', cex=2, lwd=4, lty=2) lines(Hrs_From[Z_microns == 200],HmV[Z_microns == 200], typ='b', col='green', cex=2, lwd=3, lty=3) lines(Hrs_From[Z_microns == 400],HmV[Z_microns == 400], typ='b', col='blue', cex=2, lwd=2, lty=4) lines(Hrs_From[Z_microns == 600],HmV[Z_microns == 600], typ='b', col='orange', cex=2, lwd=2, lty=5) lines(Hrs_From[Z_microns == 800],HmV[Z_microns == 800], typ='b', col='purple', cex=2, lwd=2, lty=6) par(new=TRUE, ann=FALSE) minY<- -5; maxY<- 18 plot(c(minHr, maxHr),c(minY, maxY),typ='n', yaxt='n') lines(Hrs_From[Z_microns == 0],CamV[Z_microns == 0], typ='b', col='black', cex=2, lwd=4, lty=2) lines(Hrs_From[Z_microns == 200],CamV[Z_microns == 200], typ='b', col='brown', cex=2, lwd=3, lty=3) lines(Hrs_From[Z_microns == 400],CamV[Z_microns == 400], typ='b', col='grey', cex=2, lwd=2, lty=4) lines(Hrs_From[Z_microns == 600],CamV[Z_microns == 600], typ='b', col='purple', cex=2, lwd=1, lty=5) lines(Hrs_From[Z_microns == 800],CamV[Z_microns == 800], typ='b', col='purple', cex=2, lwd=1, lty=6) browser() detach(XX)