In tdanker/ephys2: Read, analyze and plot HEKA Patchmaster files

rm(CURSORS)
library(ephys2)

library(tidyverse)
hergDRC<-get_treeinfo(examplefile("hergDRC"))
set_cursor('hERG-drc','peak',curMax_,c(2.3,2.4))
set_cursor('hERG-drc','peak2',curMax_,c(3.3,3.4))
set_cursor('hERG-drc','base',curMean_,c(0.18,0.22))
set_resultmethod("hERG-drc","corr", function(x) {x$corr<-x$peak-x$base;x})
theme_set(theme_classic())
s <-  getSeries(hergDRC, 1,1,2)
s$results() %>% names

s$results()
# columns not here but in s$sweeps:
# filename

s$sweeps %>% names

# those columns of $resutls that should befiltered are also in $sweeps
all(s$sweeps$trace == s$results()$trace)
all(s$sweeps$Trace_Time == s$results()$Trace_Time)
all(s$sweeps$tracename == s$results()$tracename)
all(s$sweeps$Serieslabel == s$results()$serlabel)

# not all are consistent:
#all(s$sweeps$filename == s$results()$###)
all(s$sweeps$exp == s$results()$exp_)
#all(s$sweeps$ser == s$results()$ser)  #ser ==> ser+serlabel, s.u.

#all(s$sweeps$Stimulus == s$results()$###)
#all(s$sweeps$ExperimentLabel == s$results()$explabel) NOT true

# utility to select sweeps. things like auto-pick "labelled" or "lpends" should go here!
select_sweeps<-function(s,selvec){
  s2<-proto(s)# copy of s. since it is a proto, it would be cahnged otherwise!
  s2$sweeps$x<-s2$sweeps$x[,selvec]
s2$sweeps$y<-s2$sweeps$y[,selvec]
s2$sweeps$Trace_Time <-s2$sweeps$Trace_Time[selvec]
s2
}

s %>% select_sweeps(c(1,3,13)) ->s2 ; s2$results()#swp and Sweep are recalulated which is not good (now 1.2.3)