In RobWHickman/gettyR:

library(tidyverse)
library(patchwork)
library(gettyR)

dir <- "C:/Users/DHill/Desktop/robert_cell_sorting/sessions"
session <- "U-2020-07-30-s2"
save_dir <- file.path(dir, session, "figures")
if(!dir.exists(save_dir)) {
  dir.create(save_dir)
}

if(!dir.exists(file.path(save_dir, "unsorted_responses"))) {
  dir.create(file.path(save_dir, "unsorted_responses"))
}
if(!dir.exists(file.path(save_dir, "sorted_responses"))) {
  dir.create(file.path(save_dir, "sorted_responses"))
}
if(!dir.exists(file.path(save_dir, "regressions"))) {
  dir.create(file.path(save_dir, "regressions"))
}
if(!dir.exists(file.path(save_dir, "cell_properties"))) {
  dir.create(file.path(save_dir, "cell_properties"))
}

spike_data <- gettyR::open_spike_data(dir, session) 
spike_data <- spike_data %>%
  select(trial, spikes, sorted_spikes) %>%
  pmap_df(gettyR::rough_match_spikes) %>%
  right_join(spike_data, by = "trial")
trace_data <- gettyR::open_cell_trace(dir)

p1 <- gettyR::plot_cell_trace(trace_data)
ggsave(file.path(save_dir, "cell_properties", "trace_plot.png"), p1)

p2 <- gettyR::plot_isi_histogram(spike_data$sorted_spikes)
ggsave(file.path(save_dir, "cell_properties", "isi_plot.png"), p2)

p3 <- gettyR::plot_getty_responses(spike_data)
ggsave(file.path(save_dir, "getty_spikes_plot.png"), p3)

RobWHickman/gettyR documentation built on Sept. 16, 2020, 10:54 p.m.