Dev/plot_comparison_multi.R

In pnojai/rwalk: Random Walk Model of Neurotransmitter Release and Uptake

# Random Walk.
# Build CV model and plot it.
# Plot data from file.
# Report fit in terms of r-squared.

library(tidyverse)

fil <- "input/181015_10mg-kgAMPH_50mM Nimo_2_outlier_scrub.csv"
sample_rate <- 100

dat <- read_experiment_csv(fil, sr = sample_rate)
dat_list <- split_stims(dat)

vmax <- 0
km <- .4
pulses <- 30
pulse_freq <- 50
release <- 3.05
bin_size <- 2.0
electrode_distance <- 50
dead_space_distance <- 4
diffusion_coefficient <- 2.7 * 10^-6
convert_current = TRUE
calibration_current = 7500.0
calibration_concentration = 5.0

lapply(dat_list, function(x) {
        compare_pulse(x, vmax = vmax, km = km, pulses = pulses, pulse_freq = pulse_freq,
                      release = release, bin_size = bin_size,
                      electrode_distance = electrode_distance, dead_space_distance = dead_space_distance,
                      diffusion_coefficient = diffusion_coefficient, convert_current = convert_current,
                      calibration_current = calibration_current, calibration_concentration = calibration_concentration)
})

# compare_pulse(dat, vmax = vmax, km = km, pulses = pulses, pulse_freq = pulse_freq,
#               release = release, bin_size = bin_size,
#               electrode_distance = electrode_distance, dead_space_distance = dead_space_distance,
#               diffusion_coefficient = diffusion_coefficient, convert_current = convert_current,
#               calibration_current = calibration_current, calibration_concentration = calibration_concentration)