sta
Compute the spike-triggered average (STA) of the stimulus
Defaults
For ndata files of type "PEP", the second sync pulse is used as the start time if 'stimchan' is NULL (unset).
Options:
- simple - ??
- plot - STA plot
- pca - Principle component analysis
- ars - Adaptive rate scaling
Parameter definitions:
- sampling - (data file value).
- sta_start - (-300) Start time for averaging window.
- sta_period - (300) Duration of averaging window.
- bin_info_flag - (0) 1-estimate mutual info and bits/spk from STA, for binary stimuli only.
- fft_flag - (0) Determines stimulus interpolation
- 0 - simple boxcar, piecewise flat
- 1 - fft interplation
- 2 - impulse
- 3 - piecewise straight (triangular)
- norm_flag - (1) Normalize the stimulus.
- negate_stim - (0) multiply stimulus by -1.0 before processing.
- invert_flag - (0) invert stimulus.
- sr_resp_flag - (0) scale and invert STA to show impulse response in spikes/sec/TFunit
- sr_resp_write_conv - (0) write the nd file with an added channel "conv" that contains the resulting convolution of the stimulus with the STA.
- deconv_flag - (0)
- pos_flag - (0) change direction to position.
- metric - (0) custum metrics for 1,2, ..., sometimes writes to metric_file
- metric_file - ("zzz.zzz") Write stats to this file.
- peak_stat - (0) if 1, write basic STA peak stats to outfile.
- varcond_type - ((char *)NULL) Types of condition STA analyses
- abs_val_avg_min - mean absolute value in [c0,c0+cn) >= vcrit
- abs_val_avg_max - mean absolute value in [c0,c0+cn) <= vcrit
- sqr_avg_min - mean squared value in [c0,c0+cn) >= vcrit
- sqr_avg_max - mean squared value in [c0,c0+cn) <= vcrit
- var_min - variance in [c0,c0+cn) >= vcrit
- var_max - variance in [c0,c0+cn) <= vcrit
- val_min - all values in [c0,c0+cn) >= vcrit
- val_max - all values in [c0,c0+cn) <= vcrit
- varcond_c0 - (sta_start) start of time window for varcond (c0 and cn must define a window that fits in the sta window)
- varcond_cn - (sta_period) duration of time window for varcond
- varcond_vcrit - (0.0) numerical criterion for varcond
- ft_sta_flag - (0) compute FT of STA.
- ft_phase_deg_flag - (1) phase in degrees (1) or radians (0).
- ft_win_center - (-50) Center for Gaussian window to multiply data.
- ft_win_sigma - (0.0) SD for Gaussian window to multiply data.
Example:
# # STA # sta all plot ...