gating process

gating_atf<-read.csv("19417001.atf", skip=2, sep='\t')
gating_atf_short = gating_atf[seq(1, nrow(gating_atf), 10), ]
names(gating_atf_short)<-c("Time (ms)", "Vj (mV)", "current (pA)")
gating_atf_short$gj<--gating_atf_short$current/gating_atf_short$Vj
#gating_atf_short$gj[gating_atf_short$gj > 20] = 20
Vj<-append(gating_atf_short$Vj, -gating_atf_short$Vj)
gj<-append(gating_atf_short$gj, gating_atf_short$gj)
plot(Vj[c(500:1400,2400:3300)],gj[c(500:1400,2400:3300)], xlab="Vj, mV", ylab="gj, nS")