Mismatch between traces and preview #161
Description
If you have a usage or support question, come join our Gitter and ask your question there.
If you found a bug, please provide all of the following:
Output of jrc version
JRCLUST v4.0.0-beta-p1 "Edward"
The command you used
jrc preview D:\Dario\tmp\1097843a_g0_t0.exportedt.imec.ap.prm
jrc traces D:\Dario\tmp\1097843a_g0_t0.exportedt.imec.ap.prm
Hello,
I am using JRC4 to spikesort my data. I found a mismatch between the voltage traces plotted by the command preview and traces (look at sites 112 for an example). Shouldn't they be identical since I am using the same parameter file on the same data
?
Thank you
I have pasted below the output of the traces and preview command and the parameter file I am using.
% JRCLUST parameters
% For a description of these parameters, including legal options, see https://jrclust.readthedocs.io/en/latest/parameters/index.html
% USAGE PARAMETERS
outputDir = 'D:\Dario\tmp'; % Directory in which to place output files (Will output to the same directory as this file if empty)
batchMode = 1; % Suppress message boxes in favor of console messages
verbose = 1; % (formerly fVerbose) Be chatty when processing
% EXECUTION PARAMETERS
gpuLoadFactor = 5; % GPU memory usage factor (Use 1/gpuLoadFactor amount of GPU memory)
maxBytesLoad = []; % (formerly MAX_BYTES_LOAD) Maximum number of bytes to load into memory
maxSecLoad = []; % (formerly MAX_LOAD_SEC) Maximum sample duration (in s) to load into memory (Overrides maxBytesLoad if nonempty)
nThreadsGPU = 128; % (formerly nThreads) Number of GPU threads to use for clustering
ramToGPUFactor = 8; % (formerly nLoads_gpu) Ratio of RAM to GPU memory
randomSeed = 0; % Seed for the random number generator
useGPU = 1; % (formerly fGpu) Use GPU where appropriate
useParfor = 1; % (formerly fParfor) Use parfor where appropriate
% PROBE PARAMETERS
probePad = [12, 12]; % (formerly vrSiteHW) Recording contact pad size (in μm) (Height x width)
shankMap = [1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]; % (formerly viShank_site) Shank ID of each site
siteLoc = [16, 0; 48, 0; 0, 20; 32, 20; 16, 40; 48, 40; 0, 60; 32, 60; 16, 80; 48, 80; 0, 100; 32, 100; 16, 120; 48, 120; 0, 140; 32, 140; 16, 160; 48, 160; 0, 180; 32, 180; 16, 200; 48, 200; 0, 220; 32, 220; 16, 240; 48, 240; 0, 260; 32, 260; 16, 280; 48, 280; 0, 300; 32, 300; 16, 320; 48, 320; 0, 340; 32, 340; 48, 360; 0, 380; 32, 380; 16, 400; 48, 400; 0, 420; 32, 420; 16, 440; 48, 440; 0, 460; 32, 460; 16, 480; 48, 480; 0, 500; 32, 500; 16, 520; 48, 520; 0, 540; 32, 540; 16, 560; 48, 560; 0, 580; 32, 580; 16, 600; 48, 600; 0, 620; 32, 620; 16, 640; 48, 640; 0, 660; 32, 660; 16, 680; 48, 680; 0, 700; 32, 700; 16, 720; 48, 720; 0, 740; 16, 760; 48, 760; 0, 780; 32, 780; 16, 800; 48, 800; 0, 820; 32, 820; 16, 840; 48, 840; 0, 860; 32, 860; 16, 880; 48, 880; 0, 900; 32, 900; 16, 920; 48, 920; 0, 940; 32, 940; 16, 960; 48, 960; 0, 980; 32, 980; 16, 1000; 48, 1000; 0, 1020; 32, 1020; 16, 1040; 48, 1040; 0, 1060; 32, 1060; 16, 1080; 48, 1080; 0, 1100; 32, 1100; 48, 1120; 0, 1140; 32, 1140; 16, 1160; 48, 1160; 0, 1180; 32, 1180; 16, 1200; 48, 1200; 0, 1220; 32, 1220; 16, 1240; 48, 1240; 0, 1260; 32, 1260; 16, 1280; 48, 1280; 0, 1300; 32, 1300; 16, 1320; 48, 1320; 0, 1340; 32, 1340; 16, 1360; 48, 1360; 0, 1380; 32, 1380; 16, 1400; 48, 1400; 0, 1420; 32, 1420; 16, 1440; 48, 1440; 0, 1460; 32, 1460; 16, 1480; 48, 1480; 0, 1500; 16, 1520; 48, 1520; 0, 1540; 32, 1540; 16, 1560; 48, 1560; 0, 1580; 32, 1580; 16, 1600; 48, 1600; 0, 1620; 32, 1620; 16, 1640; 48, 1640; 0, 1660; 32, 1660; 16, 1680; 48, 1680; 0, 1700; 32, 1700; 16, 1720; 48, 1720; 0, 1740; 32, 1740; 16, 1760; 48, 1760; 0, 1780; 32, 1780; 16, 1800; 48, 1800; 0, 1820; 32, 1820; 16, 1840; 48, 1840; 0, 1860; 32, 1860; 48, 1880; 0, 1900; 32, 1900; 16, 1920; 48, 1920; 0, 1940; 32, 1940; 16, 1960; 48, 1960; 0, 1980; 32, 1980; 16, 2000; 48, 2000; 0, 2020; 32, 2020; 16, 2040; 48, 2040; 0, 2060; 32, 2060; 16, 2080; 48, 2080; 0, 2100; 32, 2100; 16, 2120; 48, 2120; 0, 2140; 32, 2140; 16, 2160; 48, 2160; 0, 2180; 32, 2180; 16, 2200; 48, 2200; 0, 2220; 32, 2220; 16, 2240; 48, 2240; 0, 2260; 16, 2280; 48, 2280; 0, 2300; 32, 2300; 16, 2320; 48, 2320; 0, 2340; 32, 2340; 16, 2360; 48, 2360; 0, 2380; 32, 2380; 16, 2400; 48, 2400; 0, 2420; 32, 2420; 16, 2440; 48, 2440; 0, 2460; 32, 2460; 16, 2480; 48, 2480; 0, 2500; 32, 2500; 16, 2520; 48, 2520; 0, 2540; 32, 2540; 16, 2560; 48, 2560; 0, 2580; 32, 2580; 16, 2600; 48, 2600; 0, 2620; 32, 2620; 48, 2640; 0, 2660; 32, 2660; 16, 2680; 48, 2680; 0, 2700; 32, 2700; 16, 2720; 48, 2720; 0, 2740; 32, 2740; 16, 2760; 48, 2760; 0, 2780; 32, 2780; 16, 2800; 48, 2800; 0, 2820; 32, 2820; 16, 2840; 48, 2840; 0, 2860; 32, 2860; 16, 2880; 48, 2880; 0, 2900; 32, 2900; 16, 2920; 48, 2920; 0, 2940; 32, 2940; 16, 2960; 48, 2960; 0, 2980; 32, 2980; 16, 3000; 48, 3000; 0, 3020; 16, 3040; 48, 3040; 0, 3060; 32, 3060; 16, 3080; 48, 3080; 0, 3100; 32, 3100; 16, 3120; 48, 3120; 0, 3140; 32, 3140; 16, 3160; 48, 3160; 0, 3180; 32, 3180; 16, 3200; 48, 3200; 0, 3220; 32, 3220; 16, 3240; 48, 3240; 0, 3260; 32, 3260; 16, 3280; 48, 3280; 0, 3300; 32, 3300; 16, 3320; 48, 3320; 0, 3340; 32, 3340; 16, 3360; 48, 3360; 0, 3380; 32, 3380; 48, 3400; 0, 3420; 32, 3420; 16, 3440; 48, 3440; 0, 3460; 32, 3460; 16, 3480; 48, 3480; 0, 3500; 32, 3500; 16, 3520; 48, 3520; 0, 3540; 32, 3540; 16, 3560; 48, 3560; 0, 3580; 32, 3580; 16, 3600; 48, 3600; 0, 3620; 32, 3620; 16, 3640; 48, 3640; 0, 3660; 32, 3660; 16, 3680; 48, 3680; 0, 3700; 32, 3700; 16, 3720; 48, 3720; 0, 3740; 32, 3740; 16, 3760; 48, 3760; 0, 3780; 16, 3800; 48, 3800; 0, 3820; 32, 3820]; % (formerly mrSiteXY) Site locations (in μm) (x values in the first column, y values in the second column)
siteMap = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96, 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 107, 108, 109, 110, 111, 112, 114, 115, 116, 117, 118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143, 144, 145, 146, 147, 148, 149, 150, 151, 153, 154, 155, 156, 157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182, 183, 184, 185, 186, 187, 188, 190, 191, 192, 193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207, 208, 209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221, 222, 223, 224, 225, 226, 227, 229, 230, 231, 232, 233, 234, 235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247, 248, 249, 250, 251, 252, 253, 254, 255, 256, 257, 258, 259, 260, 261, 262, 263, 264, 266, 267, 268, 269, 270, 271, 272, 273, 274, 275, 276, 277, 278, 279, 280, 281, 282, 283, 284, 285, 286, 287, 288, 289, 290, 291, 292, 293, 294, 295, 296, 297, 298, 299, 300, 301, 302, 303, 305, 306, 307, 308, 309, 310, 311, 312, 313, 314, 315, 316, 317, 318, 319, 320, 321, 322, 323, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 336, 337, 338, 339, 340, 342, 343, 344, 345, 346, 347, 348, 349, 350, 351, 352, 353, 354, 355, 356, 357, 358, 359, 360, 361, 362, 363, 364, 365, 366, 367, 368, 369, 370, 371, 372, 373, 374, 375, 376, 377, 378, 379, 381, 382, 383, 384]; % (formerly viSite2Chan) Map of channel index to site ID (The mapping siteMap(i) = j corresponds to the statement 'site i is stored as channel j in the recording')
% RECORDING FILE PARAMETERS
bitScaling = 2.3438; % (formerly uV_per_bit) ADC bit scaling factor (Conversion factor for ADC bit values to μV)
dataType = 'int16'; % (formerly vcDataType) Format of raw recordings
headerOffset = 0; % (formerly header_offset) Recording file header offset (in bytes)
nChans = 385; % Number of channels stored in recording file (Distinct from the number of AP sites)
rawRecordings = {'D:\Dario\tmp\1097843a_g0_t0.exportedt.imec.ap.bin'}; % Path or paths to raw recordings to sort
recordingFormat = 'SpikeGLX'; % Format of raw recording
sampleRate = 30000; % (formerly sRateHz) Sampling rate (in Hz) of raw recording
tallSkinny = 1; % (formerly fTranspose_bin) Recording will be interpreted as nChannels x nSamples if true
% PREPROCESSING PARAMETERS
blankThresh = 60; % (formerly blank_thresh) Threshold (in MADs) above which to reject samples exceeding channel median after filtering
filterType = 'ndiff'; % (formerly vcFilter) Type of filter to use on raw data
filtOrder = 3; % Bandpass filter order
freqLimBP = [300, 3000]; % (formerly freqLim) Frequency cutoffs for bandpass filter
ignoreChans = []; % (formerly viChanZero) Channel numbers to ignore manually
ignoreSites = []; % (formerly viSiteZero) Site IDs to ignore manually
fftThresh = 0; % (formerly fft_thresh) Threshold (in MADs of power-frequency product) above which to remove frequency outliers
freqLimNotch = []; % Frequency ranges to exclude for notch filter
freqLimStop = []; % Frequency range to exclude for band-stop filter
gainBoost = 1; % (formerly gain_boost) Scale factor to boost gain in raw recording (Used in filtering operation)
loadTimeLimits = []; % (formerly tlim_load) Time range (in s) of samples to load at once (All samples are loaded if empty)
nDiffOrder = 2; % (formerly nDiff_filt) Order for differentiator filter (Used if and only if filterType is 'sgdiff' or 'ndiff')
nSamplesPad = 100; % (formerly nPad_filt) Number of samples to overlap between chunks in large files
useElliptic = 1; % (formerly fEllip) Use elliptic (bandpass) filter if true (Uses Butterworth filter if false)
userFiltKernel = []; % (formerly vnFilter_user) User-specified filter kernel (Ignored unless filterType is 'user')
% SPIKE DETECTION PARAMETERS
blankPeriod = 0.3; % (formerly blank_period_ms) Duration of blanking period (in ms) when the common mean exceeds blankThresh
CARMode = 'median'; % (formerly vcCommonRef) The meaning of 'average' in 'common average reference'
evtDetectRad = 50; % (formerly maxDist_site_um) Maximum distance (in μm) to search over for duplicate peaks
evtWindow = [-0.25, 0.75]; % (formerly spkLim_ms) Time range (in ms) of filtered spike waveforms, centered at the peak
nSiteDir = 6.5; % (formerly maxSite) Number of neighboring sites to group in either direction (nSitesEvt is set to 1 + 2nSiteDir - nSitesExcl)
nSitesExcl = 6; % (formerly nSites_ref) Number of sites to exclude from the spike waveform group for feature extraction
qqFactor = 5; % Spike detection threshold factor (Thr = qqFactormed(abs(x-med(x)))/0.6745)
refracInt = 0.25; % (formerly spkRefrac_ms) Spike refractory period (in ms)
detectBipolar = 0; % (formerly fDetectBipolar) Detect positive as well as negative peaks
evtManualThresh = []; % (formerly spkThresh_uV) Manually-set spike detection threshold (in μV)
evtWindowRaw = [-0.5, 1.5]; % (formerly spkLim_raw_ms) Time range (in ms) of raw spike waveforms, centered at the peak
groupShank = 0; % (formerly fGroup_shank) Group all sites on the same shank if true
minNeighborsDetect = 0; % (formerly nneigh_min_detect) Minimum number of sample neighbors exceeding threshold for a sample to be considered a peak
realignTraces = 0; % (formerly fRealign_spk) Realign spike traces after subtracting local CAR (Realign if 1, perform subpixel interpolation if 2)
spikeThreshMax = []; % (formerly spkThresh_max_uV) Maximum absolute amplitude (in μV) permitted for spikes
threshFile = ''; % (formerly vcFile_thresh) Path to .mat file storing the spike detection threshold (Created by preview GUI)
% FEATURE EXTRACTION PARAMETERS
clusterFeature = 'pca'; % (formerly vcFet) The feature to extract from your spike waveforms in order to cluster them
evtGroupRad = 75; % (formerly maxDist_site_spk_um) Maximum distance (in μm) for extracting spike waveforms
nPCsPerSite = 1; % (formerly nPcPerChan) Number of principal components to compute per site
extractAfterDetect = 0; % Extract features only after detecting all spikes across all files if true (This is effectively set to true if you specify clusterFeature='gpca')
interpPC = 1; % (formerly fInterp_fet) Interpolate 1st principal vector to maximize projection of spikes if true
nPeaksFeatures = 2; % (formerly nFet_use) Number of potential peaks to use when computing features
% CLUSTERING PARAMETERS
autoMergeBy = 'pearson'; % (formerly autoMergeCriterion) Metric to use for automerging clusters based on average waveform
distCut = 2; % (formerly dc_percent) Percentile of pairwise distances between spikes on a site to use as a cutoff distance
evtMergeRad = 35; % (formerly maxDist_site_merge_um) Maximum distance (in μm) to consider for merging spike waveforms
log10DeltaCut = 0.6; % (formerly delta1_cut) Log10 of delta cutoff (Spikes with delta values below this cutoff will not be considered as cluster centers)
log10RhoCut = -2.5; % (formerly rho_cut) Log10 of rho cutoff (Spikes with rho values below this cutoff will not be considered as cluster centers)
maxUnitSim = 0.98; % (formerly maxWavCor) Threshold for merging two units having similar spike waveforms (Units with a similiarity score above this value will be merged)
minClusterSize = 30; % (formerly min_count) Minimum number of spikes per cluster (Automatically set to the maximum of this value and twice the number of features)
nClusterIntervals = 4; % (formerly nTime_clu) Number of intervals to divide the recording into around a spike (When clustering, take the 1/nClusterIntervals fraction of all spikes around a spiking event to compute distance)
RDDetrendMode = 'global'; % (formerly vcDetrend_postclu) Detrending mode to apply to rho-delta values in order to determine cluster centers
driftMerge = 1; % (formerly fDrift_merge) Compute multiple waveforms at three drift locations based on the spike position if true
evtWindowMergeFactor = 1; % (formerly spkLim_factor_merge) Ratio of samples to take when computing correlation
maxClustersSite = 20; % (formerly maxCluPerSite) Maximum number of cluster centers computed per site (Used if RDDetrendMode is 'local')
meanInterpFactor = 1; % (formerly nInterp_merge) Interpolation factor for mean unit waveforms (Set to 1 to disable)
minSitesWeightFeatures = 5; % (formerly min_sites_mask) Minimum number of sites to have if using weightFeatures (Ignored if weightFeatures is false)
nPassesMerge = 10; % (formerly nRepeat_merge) Number of times to repeat automatic waveform-based merging
outlierThresh = 7.5; % (formerly thresh_mad_clu) Threshold (in MADs) to remove outlier spikes for each cluster
useGlobalDistCut = 0; % (formerly fDc_global) Use a global distance cutoff for all sites if true
weightFeatures = 0; % (formerly fSpatialMask_clu) Weight display features by distance from site if true
% CURATION PARAMETERS
figList = {'FigCorr', 'FigHist', 'FigISI', 'FigMap', 'FigPos', 'FigProj', 'FigRD', 'FigSim', 'FigTime', 'FigWav'}; % List of tags of figures to display in feature view
frFilterShape = 'triangle'; % (formerly filter_shape_rate) Kernel shape for temporal averaging (Used in estimation of the firing rate of a given unit)
frPeriod = 2; % (formerly filter_sec_rate) Time period (in s) over which to determine firing rate (Used in estimation of the firing rate of a given unit)
frSampleRate = 1000; % (formerly sRateHz_rate) Resampling rate (in Hz) for estimating the firing rate (Used in estimation of the firing rate of a given unit)
% DISPLAY PARAMETERS
dispTimeLimits = [0, 0.2]; % (formerly tlim) Time range (in ms) to display
nSpikesFigProj = 500; % (formerly nShow_proj) Maximum number of spikes per cluster to display in the feature projection view
nSpikesFigWav = 30; % (formerly nSpk_show) Maximum number of spikes per cluster to display generally
colorMap = [0.5, 0.5, 0.5; 0, 0, 0; 1, 0, 0]; % (formerly mrColor_proj) RGB color map for background, primary selected, and secondary selected spikes (The first three values are the R values, the next three are the G values, and the last three are the B values.)
corrRange = [0.9, 1]; % (formerly corrLim) Correlation score range to distinguish by color map
dispFeature = 'vpp'; % (formerly vcFet_show) Feature to display in the feature projection plot
dispFilter = 'ndiff'; % (formerly vcFilter_show) Filter to apply in traces plot
maxAmp = 250; % Amplitude scale (in μV)
nSitesFigProj = 5; % Number of sites to show in feature projection view
nSpikesFigISI = 200; % Maximum number of spikes to show in ISI view
pcPair = [1, 2]; % Pair of PCs to display
projTimeLimits = []; % (formerly tLimFigProj) Time range (in s) to display in feature projection view
showRaw = 0; % (formerly fWav_raw_show) Show raw traces in waveform view if true
showSpikeCount = 1; % (formerly fText) Show spike count per unit in waveform plot
umPerPix = 20; % (formerly um_per_pix) Vertical site center-to-center spacing
% TRIAL PARAMETERS
psthTimeLimits = [-0.025, 0.05]; % (formerly tlim_psth) Time range (in s) over which to display PSTH
trialFile = ''; % (formerly vcFile_trial) Path to file containing trial data (Can be .mat or .csv, must contain timestamps of trials in units of s)
psthTimeBin = 0.001; % (formerly tbin_psth) Time bin (in s) for PSTH view
psthXTick = 0.005; % (formerly xtick_psth) PSTH time tick mark spacing
% VALIDATION PARAMETERS
gtFile = ''; % (formerly vcFile_gt) Path to file containing ground-truth data
% PREVIEW PARAMETERS
nLoadsMaxPreview = 30; % (formerly nLoads_max_preview) Number of time segments to load in preview
nSecsLoadPreview = 1; % (formerly sec_per_load_preview) Number of seconds to load in preview
siteCorrThresh = 0; % (formerly thresh_corr_bad_site) Threshold to reject bad sites based on maximum correlation with neighboring sites (Set to 0 to disable)
% TRACES PARAMETERS
nSkip = 1; % (formerly nSkip_show) Show every nSkip samples when plotting traces
nSegmentsTraces = 1; % (formerly nTime_traces) Number of time segments to display in traces view (A value of 1 shows one continuous time segment)
% LFP PARAMETERS
lfpSampleRate = 2500; % (formerly sRateHz_lfp) Sampling rate for LFP channels
% AUX CHANNEL PARAMETERS
auxChan = []; % (formerly iChan_aux) Auxiliary channel index
auxFile = ''; % (formerly vcFile_aux) Path to file containing auxiliary channel
auxLabel = 'Aux channel'; % (formerly vcLabel_aux) Label for auxiliary channel data
auxSampleRate = []; % (formerly sRateHz_aux) Sample rate for auxiliary file
auxScale = 1; % (formerly vrScale_aux) Scale factor for aux data
nClustersShowAux = 10; % (formerly nClu_show_aux) Number of clusters to show in the aux vs. firing rate correlation