Customizing the model:(back to readme) The comments included in plot_struct.m should provide sufficient information for easy customization of the script. Below, we will draw attention to several variables that users may wish to modify. plot_raw = 1; plot_meanbar = 1; Turns on plot of raw traces and also bar graphs showing mean value of raw traces (0=off, 1=on). channel_list = {'hill2', ...} exp_channel=channel_list{X}; Cell array "channel_list" contains the names of all the channels implemented in our minimal model. To select a channel to simulate, change X to the corresponding number (X=1-6). Choices are: 'hill2s' - Hill equation determining nonstationary gating of SK2 'vTau2s' - A typical 2-state Markov process. 'sAHP5s' - Implementation of the sAHP channel (eqn 11 in Stanley et al), without the final O state 'sAHP6s' - Full implementation of sAHP (eqn 11) 'SK2_6s' - SK2 channel kinetic scheme for low open probability gating regime (eqn 6) 'SK2h_6s' - SK2 channel kinetic scheme for high open probability gating regime (eqn 7) scalingfactors = [1 2 4 8 16 32 64 128]; Specifies scaling factor by which to scale pulse height and spacing. See "Model details" section above and comments in plot_struct.m. p.tf = 10; This is the stop time (in seconds) of the simulation. Here we use 10 seconds. Default for Figure 9 is 180 seconds (very slow). p.stats_start = 3.1; When calculating the mean current through a channel, the first 3.1 seconds are discarded to allow for initial settling time. This value should be changed according to the settling time of the channel under consideration. p.rate_scale = 1; Scales all rate constants of the channel's kinetic scheme by this value. Setting this to 100 generates the "fast" traces in Figure 9b and 9c. Good luck and have fun!