R/mini_rescale_start.R

In brentscott93/biophysr: Analysis tools for muscle biophysics

#' Find better event start
#'
#' @param raw_data
#' @param run_mean
#'
#' @return
#' @export
#'
#' @examples
mini_rescale_start <- function(raw_data, run_mean){
  #Determine if run_mean is in an event or baseline noise by using >10 as event
  on_off <- ifelse(run_mean > 10, 2, 1)

  rle_object<- as_tibble(do.call("cbind", rle(on_off)))

  #### FIND EVENT DURATIONS ####
  #If the rle_object's last row is in state 1, get rid of that last row. This needs to end in state 2 to capture the end of the last event
  mini_rle_object <- if(tail(rle_object, 1)$values == 1){
    slice(rle_object, -length(rle_object$values))
  } else {
    rle_object
  }

  split_data <- mini_rle_object %>%
    dplyr::mutate(cumsum = cumsum(lengths)) %>%
    dplyr::group_by(values) %>%
    split(mini_rle_object$values)

  #data is recmombined in a state_1 column and a state_2
  #the values in these columns represent the last data point in either state 1 or state 2
  #So the range of values between the end of state 1 (or start of state 2) and the end of state 2 is the event duration
  regroup_data <- bind_cols(state_1_end = split_data[[1]]$cumsum, state_2_end = split_data[[2]]$cumsum) %>%
    mutate(filtered_event_duration_dp = state_2_end - state_1_end)

  #filter out state 2s that are less than 10 ms (50 data points)
  events <- regroup_data %>%
    filter(filtered_event_duration_dp > 50)


}