R/ElectroProject.R

In Tobias-Ruff/relectro_bak: Analysis of electrophysiological data

#' A S4 class to represent an electrophysiological project containing several recording sessions.
#'
#' An ElectroProject object contains a group of RecSession objects. It can be used to run functions on several recording sessions.
#'
#' It assumes that your data are organized in a specific way on your hard drive.
#' There should be a top directory containing one directory for each subject of the experiment.
#' Within each subject directory, there should be a directory for each recording session with this subject.
#' By convention, the names of the session directories contain at least one hyphen.
#' There should be no other directories with a hyphen in their names.
#'
#' This object can be used to get an overview of the progress during the data acquisition period.
#'
#' Two useful methods are getSessionList and runOnSessionList
#'
#' @slot directory Top directory of the project
#' @slot resultsDirectory Directory where results will be saved by default
#' @slot nSessions Number of recording sessions in the project
#' @slot sessionNameList Name of the recording sessions
#' @slot sessionPathList Path to all the recording sessions
#' @slot sessionList List of RecSession objects.
ElectroProject <- setClass(
  "ElectroProject", ## name of the class
  slots=c(directory="character",
          resultsDirectory="character",
          nSessions="numeric",
          sessionNameList="character",
          sessionPathList="character",
          sessionList="list"),
  prototype = list(directory=""))


#' Create a list of RecSession objects for the ElectroProject
#'
#' Will list directories in the project directories.
#' Only directories with an hyphen in their names are considered recSession directories.
#' So only recording session directories should have a hyphen in their names
#'
#' @param ep ElectroProject object
#' @param loadSessions logical, whether the RecSession object should be created
#' @return ElectroProject object with a list of session names and directories
#'
#' @docType methods
#' @rdname setSessionList-methods
setGeneric(name="setSessionList",
           def=function(ep,loadSessions=TRUE)
           {standardGeneric("setSessionList")}
)
#' @rdname setSessionList-methods
#' @aliases setSessionList,ANY,ANY-method
setMethod(f="setSessionList",
          signature="ElectroProject",
          definition=function(ep,loadSessions=TRUE)
          {
            if(ep@directory=="")
              stop("ep@directory not set")

            ## remove any "/" at the end of the directory name
            if(substr(ep@directory,nchar(ep@directory),nchar(ep@directory))=="/")
              ep@directory<-substr(ep@directory,1,nchar(ep@directory)-1)

            ep@resultsDirectory=paste(ep@directory,"results",sep="/")

            if(!file.exists(ep@resultsDirectory)){
              dir.create(path=ep@resultsDirectory)
            }

            ## list all directories in the project path
            ## this is really slow if over nfs.
            dirs<-list.dirs(path=ep@directory)
            lastLevelNames<-unlist(lapply(strsplit(dirs,"/"),function(x){x[length(x)]}))

            ## only keep the directories with 2 hyphens in the name of the last level
            dirs<-dirs[lengths(gregexpr("-", lastLevelNames))==2]
            ep@nSessions<-length(dirs)
            ep@sessionPathList<-dirs
            ## make sure all directories have the same depth
            x<-sapply(ep@sessionPathList,function(x){length(unlist(strsplit(x,split="/")))})
            xMode<-modeRelectro(as.numeric(x))
            if(any(x!=xMode)){
              print(paste("The depth of some directories differs"))
              print(paste("Most directories have a depth of",xMode))
              print(paste("For example,",head(ep@sessionPathList[which(x==xMode)],n=1)))
              print("The following directories will not be included")
              print(ep@sessionPathList[which(x!=xMode)])
              ep@sessionPathList<-ep@sessionPathList[which(x==xMode)]
              ep@nSessions<-length(ep@sessionPathList)
            }
            dirDepth<-xMode
            ep@sessionNameList<-
              unlist(strsplit(ep@sessionPathList,split="/"))[seq(from=dirDepth,to=dirDepth*ep@nSessions,by=dirDepth)]

            if(loadSessions==TRUE){
              ep<-loadSessionsInList(ep)
              }
          return(ep)
          })


#' Create a List of RecSessions objects from the SessionNameList
#'
#' The RecSession objects will be initialized.
#'
#' @param ep ElectroProject object
#' @return ElectroProject object with RecSession loaded
#'
#' @docType methods
#' @rdname loadSessionsInList-methods
setGeneric(name="loadSessionsInList",
           def=function(ep)
           {standardGeneric("loadSessionsInList")}
)
#' @rdname loadSessionsInList-methods
#' @aliases loadSessionsInList,ANY,ANY-method
setMethod(f="loadSessionsInList",
          signature="ElectroProject",
          definition=function(ep)
          {
            if(ep@directory=="")
              stop("ep@directory not set")

            if(length(ep@sessionNameList)==0)
              stop("ep@sessionNameList has length of 0")
            if(length(ep@sessionPathList)==0)
              stop("ep@sessionPathList has length of 0")

            for (i in 1:length(ep@sessionNameList))
            {ep@sessionList[[i]]<-new("RecSession",session=ep@sessionNameList[i],path=ep@sessionPathList[i])}
            ep@sessionList<-lapply(ep@sessionList,loadRecSession)
            return(ep)
          })



#' Return a RecSession object from the list in an ElectroProject objects
#'
#'
#' @param ep ElectroProject object
#' @param name Name of the RecSession, e.g. jp19841-01072015-0108
#' @return RecSession
#'
#' @docType methods
#' @rdname getRecSession-methods
setGeneric(name="getRecSession",
           def=function(ep,name)
           {standardGeneric("getRecSession")}
)
#' @rdname getRecSession-methods
#' @aliases getRecSession,ANY,ANY-method
setMethod(f="getRecSession",
          signature="ElectroProject",
          definition=function(ep,name="")
          {
            if(ep@directory=="")
              stop("ep@directory not set")
            if(!name%in%ep@sessionNameList){
              stop(paste(name,"is not in the session list"))
            }
            return(ep@sessionList[sapply(ep@sessionList,function(x,name){x@session==name},name)][[1]])
          })


#' Return a list of RecSession objects from a character vectors containing session names
#'
#' @param ep ElectroProject object
#' @param sessionNames Character vector containing session names
#' @return list of RecSession objects
#'
#' @docType methods
#' @rdname getSessionListFromSessionNames-methods
setGeneric(name="getSessionListFromSessionNames",
           def=function(ep,sessionNames)
           {standardGeneric("getSessionListFromSessionNames")}
)
#' @rdname getSessionListFromSessionNames-methods
#' @aliases getSessionListFromSessionNames,ANY,ANY-method
setMethod(f="getSessionListFromSessionNames",
          signature="ElectroProject",
          definition=function(ep,sessionNames)
          {
            if(ep@directory=="")
              stop("ep@directory not set")
            if(class(sessionNames)!="character")
              stop("sessionNames should be a character vector")
            return(ep@sessionList[which(ep@sessionNameList %in% sessionNames)])
          })


#' Return a list of clustered RecSession objects
#'
#'
#' @param ep ElectroProject object
#' @return list of clustered RecSession objects
#'
#' @docType methods
#' @rdname getClusteredSessionList-methods
setGeneric(name="getClusteredSessionList",
           def=function(ep)
           {standardGeneric("getClusteredSessionList")}
)
#' @rdname getClusteredSessionList-methods
#' @aliases getClusteredSessionList,ANY,ANY-method
setMethod(f="getClusteredSessionList",
          signature="ElectroProject",
          definition=function(ep)
          {
            if(ep@directory=="")
              stop("ep@directory not set")
            return(ep@sessionList[sapply(ep@sessionList,getIsClustered)])
          })


#' Return a list of RecSession objects that have some common properties
#'
#'
#' @param ep ElectroProject object
#' @param clustered logical indicating whether the session should be clustered or not
#' @param region Set to a given brain region to select only sessions with tetrodes in this brain region
#' @param environment Set to a given environment code to select only sessions during which this environment was presented
#' @param stim Set to a given stimulation code to select only sessions during which this stimulatoin was presented
#' @param fileExtension Keep sessions that have a session file ending with the value of fileExtension (e.g. kld-19021016-0101.fileExtension)
#' @return list of RecSession objects
#'
#' @docType methods
#' @rdname getSessionList-methods
setGeneric(name="getSessionList",
           def=function(ep,clustered="",region="",environment="",stim="",fileExtension="")
           {standardGeneric("getSessionList")}
)
#' @rdname getSessionList-methods
#' @aliases getSessionList,ANY,ANY-method
setMethod(f="getSessionList",
          signature="ElectroProject",
          definition=function(ep,clustered="",region="",environment="",stim="",fileExtension="")
          {
            if(ep@directory=="")
              stop("ep@directory not set")
            if(length(ep@sessionList)==0)
              stop("sp@sessionList has a length of 0")
            myList<-ep@sessionList
            if(clustered==T)
              myList<-myList[sapply(myList,getIsClustered)]
            if(region!=""){
              myList<-myList[sapply(myList,containsElectrodeLocation,location=region)]
            }
            if(environment!=""){
              myList<-myList[sapply(myList,containsEnvironment,environment=environment)]
            }
            if(stim!=""){
              myList<-myList[sapply(myList,containsStimulation,stimulation=stim)]
            }
            if(fileExtension!=""){
              myList<-myList[sapply(myList,fileExists,extension=fileExtension)]
            }
            return(myList)
          })



#' Running a function on a set of recording sessions
#'
#' This applies a function to a list of RecSession objects.
#' If save is set to TRUE, the results returned by the function will be saved.
#' The function applied to single recording sessions should return the results in a list.
#' The results are saved in the resultsDirectory of the ElectroProject object.
#' The names of the files saved will be the name of the elements in the list returned by the function
#' If the element is an array or matrix, they will be bound with abind.
#' If not a array or matrix, the data from each recording session will be bound togheter using rbind.
#' You can use snow::parLapply instead of lapply by setting parallel to TRUE and passing a valid cluster to the function
#' If save is set to FALSE, the data returned by the lapply function will not be saved but retured.
#' If overwrite is set to TRUE, previous data will be overwrite when saving the data.
#' If overwrite is set to FALSE, the new data will be appended to the old one.
#'
#' @param ep ElectroProject object
#' @param sessionList List of RecSession objects on which the function will be applied
#' @param fnct A function to run on each RecSession
#' @param save Whether you want to save the data returned by the function
#' @param overwrite Logical indicating if the data returned by the function will overwrite the old one when saving into a file
#' @param parallel Whether you want to run the function in parallel
#' @param cluster A cluster generated from the makeCluster function of the snow package
#' @param ... optional arguments to fnct
#' @return list of list returned by the lapply function
#'
#' @docType methods
#' @rdname runOnSessionList-methods
setGeneric(name="runOnSessionList",
            def=function(ep,sessionList,fnct=function(x){NA},save=T,overwrite=T,parallel=F,cluster="",...)
            {standardGeneric("runOnSessionList")}
          )

#' @rdname runOnSessionList-methods
#' @aliases runOnSessionList,ANY,ANY-method
setMethod(f="runOnSessionList",
         signature="ElectroProject",
         definition=function(ep,sessionList,fnct=function(x){NA},save=T,overwrite=T,parallel=F,cluster="",...)
         {
           if(class(ep)!="ElectroProject")
             stop("runOnSessionList, ep should be a ElectroProject object")
           if(!is.list(sessionList))
             stop("runOnSessionList, sessionList has to be a list")
           if(length(sessionList)==0)
             stop("runOnSessionList, sessionList has size 0")
           if(class(sessionList[[1]])!="RecSession")
             stop("runOnSessionList, sessionList should contain RecSession objects")
           if(!is.function(fnct))
             stop("runOnSessionList, fnct needs to be a function")
           if(parallel==T)
             if(class(cluster)[2]!="cluster")
               stop("runOnSessionList, give a valid snow cluster if you want to run the function on several threads")
           if(save==T){
             if(!file.exists(ep@resultsDirectory))
             {
               stop(paste("runOnSessionList,",ep@resultsDirectory,"does not exist"))
             }
           }

           if(parallel==T){
             list.res<-snow::parLapply(cluster,sessionList,fnct,...)
           } else {
             list.res<-lapply(sessionList,fnct,...)
           }

           if(save==T){
             ## check that list.res is a list of list
             if(!is.list(list.res))
               stop("runOnSessionList, list.res is not a list")
             if(!is.list(list.res[[1]]))
               stop(paste("runOnSessionList, first item of the list list.res is not a list, fnct should return a list"))

             ## list of objects to merge and save
             objectNames<-names(list.res[[1]])
             if(any(objectNames==""))
               stop(paste("runOnSessionList, not all elements of the list returned by the function have a name"))

             for(n in objectNames){
               obj<-list.res[[1]]
               if(overwrite==T){
                 if(class(obj[[n]])=="array"| class(obj[[n]])=="matrix"){
                  assign(n,do.call(abind::abind,sapply(list.res,function(x){x[n]})))  ## bind along the last dimension
                 }else{
                  assign(n,do.call("rbind", sapply(list.res,function(x){x[n]})))
                 }
               }else{## concatonate to existing data
                 if(class(obj[[n]])=="array"|class(obj[[n]])=="matrix"){
                   assign(paste(n,"new",sep="."),do.call(abind::abind, sapply(list.res,function(x){x[n]}))) ## bind along the last dimension
                   load(file=paste(ep@resultsDirectory,n,sep="/"))
                   assign(n,abind::abind(get(n),get(paste(n,"new",sep="."))))
                 }else{
                   assign(paste(n,"new",sep="."),do.call("rbind", sapply(list.res,function(x){x[n]})))
                   load(file=paste(ep@resultsDirectory,n,sep="/"))
                   assign(n,rbind(get(n),get(paste(n,"new",sep="."))))
                 }
               }
               print(paste("saving",paste(ep@resultsDirectory,n,sep="/")))
               save(list=n,file=paste(ep@resultsDirectory,n,sep="/"))
             }
           }else{
             if(!is.null(list.res[[1]]))
              return(list.res)
           }
         })

### show ###
setMethod("show", "ElectroProject",
          function(object){
            print(paste("directory:",object@directory))
            print(paste("result directory:",object@resultsDirectory))
            print(paste("nSessions:",object@nSessions))

            if(length(object@sessionNameList)!=0){
              print(paste("sessionNameList:"))
              print(object@sessionNameList)
              if(length(object@sessionList)!=0){
                m<-matrix(c(sapply(object@sessionList,getIsClustered),sapply(object@sessionList,getIsEarlyProcessed)),ncol=2)
                print(paste("Clustered sessions:",sum(m[,1])))
                print(object@sessionNameList[m[,1]])
                print(paste("Not clustered, but early processed sessions:",length(which(m[,1]==F&m[,2]==T))))
                print(object@sessionNameList[which(m[,1]==F&m[,2]==T)])
                print(paste("Not early processed sessions:", length(which(m[,1]==F&m[,2]==F))))
                print(object@sessionNameList[which(m[,1]==F&m[,2]==F)])
              }
            }
          })


#' Sort a list of RecSession objects according to the recording date of the sessions
#'
#' Use the date in the name of the RecSession object
#'
#' @param rsl List containing recording sessions
#' @return List containing the recording sessions ordered according to their recording date.
sortRecSessionListChronologically<-function(rsl){
  if(class(rsl)!="list")
    stop("sortRecSessionListChronologically: rsl is not a list")
  if(length(rsl)==0)
    stop("sortRecSessionListChronologically: rsl length is 0")
  if(class(rsl[[1]])!="RecSession")
    stop("sortRecSessionListChronologically: rsl[[1]] is not a RecSession object")
  o<-order(sapply(rsl,recordingDate))
  return(rsl[o])
}


#' Get a character vector with the session names of a list of RecSession object
#'
#' @param rss List containing recording sessions
#' @return Charcter vector of the session names
sessionNamesFromSessionList<-function(rss){
  x<-function(rs){rs@session}
  unlist(lapply(rss,x))
}



#' Create a copy of the data from an experiment
#'
#' There is an argument to specify the type of backup you want.
#' This can be use to backup different files depending on
#' what needs to be copied.
#'
#' List of files exported if copyType is set to medium
#' - .res, .clu, .par, .desen, .px_per_cm, .resofs, .desel
#' .res_samples_per_whl_sample .sampling_rate_dat, .whl .whd
#' Also includes the tetrode specific .clu. and .res. files
#'
#'
#' @param ep ElectroProject object
#' @param sessionList List of RecSession objects that will be included in the copy
#' @param copyType Type of copy. Only "medium" is supported at the moment.
#' @param extensions Character vector with the list of file types you want to export
#' @param destination Path of the directory in which the copy is made
#' @param compress Logical indicating whether to call tar to make a tar.gzip from the copy
#'
#' @docType methods
#' @rdname copyDatabase-methods
setGeneric(name="copyDatabase",
           def=function(ep,sessionList,copyType="medium",extensions,destination,compress=FALSE)
           {standardGeneric("copyDatabase")}
)

#' @rdname copyDatabase-methods
#' @aliases copyDatabase,ANY,ANY-method
setMethod(f="copyDatabase",
          signature="ElectroProject",
          definition=function(ep,sessionList,copyType="medium",extensions,destination,compress=FALSE)
          {
            if(!dir.exists(destination)){
              if(dir.create(path=destination,recursive=T)==FALSE)
                stop()
            }
            if(ep@directory=="")
              stop("copyExperiment: ep@directory is not set")
            if(class(sessionList)!="list")
              stop("copyExperiment: sessionList should be a list")
            if(class(sessionList[[1]])!="RecSession")
              stop("copyExperiment: sessionList[[1]] is not a RecSession object")
            if(copyType!="medium")
              stop("copyExperiment: copyType is not set to medium")
            if(copyType=="medium"){
              #each session will take care of copying itself.
              sessionSpecificExtensions=c("res","clu","par","desen","px_per_cm","resofs","desel",
                                          "res_samples_per_whl_sample","sampling_rate_dat","whl","whd")
              tetrodeSpecificExtensions=c("res","clu")
            }
            print(paste("destination:",destination))
            print(paste("session specific extensions:"))
            print(sessionSpecificExtensions)
            print(paste("tetrode specific extensions:"))
            print(tetrodeSpecificExtensions)
            print(paste("additional session specific extensions:"))
            print(extensions)
            ## combine all the session specific extensions
            sessionSpecificExtensions<-c(sessionSpecificExtensions,extensions)
            ## each RecSession object will copy itself ##
            res<-sapply(sessionList,copyRecSessionFiles,destination,sessionSpecificExtensions,tetrodeSpecificExtensions)

            if(compress==TRUE)
            {
              comp=paste(destination,"tar.gzip",sep=".")
              print(paste("creating",comp, "with the tar funciton"))
              tar(tarfile=comp,
                  files=destination,
                  compression="gzip")
            }
      })