-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathconfigFile_KiloSort25.m
52 lines (38 loc) · 2.07 KB
/
configFile_KiloSort25.m
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
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
% sample rate
ops.fs = 20000;
% frequency for high pass filtering (150)
ops.fshigh = 300;
% threshold on projections (like in Kilosort1, can be different for last pass like [10 4])
ops.Th = [10 4];
% how important is the amplitude penalty (like in Kilosort1, 0 means not used, 10 is average, 50 is a lot)
ops.lam = 10;
% splitting a cluster at the end requires at least this much isolation for each sub-cluster (max = 1)
ops.AUCsplit = 0.9;
% minimum spike rate (Hz), if a cluster falls below this for too long it gets removed
ops.minFR = 1/50;
% number of samples to average over (annealed from first to second value)
ops.momentum = [20 400];
% spatial constant in um for computing residual variance of spike
ops.sigmaMask = 30;
% threshold crossings for pre-clustering (in PCA projection space)
ops.ThPre = 8;
% spatial scale for datashift kernel
ops.sig = 20;
% type of data shifting (0 = none, 1 = rigid, 2 = nonrigid)
ops.nblocks = 0;
%% danger, changing these settings can lead to fatal errors
% options for determining PCs
ops.spkTh = -6; % spike threshold in standard deviations (-6)
ops.reorder = 1; % whether to reorder batches for drift correction.
ops.nskip = 25; % how many batches to skip for determining spike PCs
ops.GPU = 1; % has to be 1, no CPU version yet, sorry
% ops.Nfilt = 1024; % max number of clusters
ops.nfilt_factor = 4; % max number of clusters per good channel (even temporary ones)
ops.ntbuff = 64; % samples of symmetrical buffer for whitening and spike detection
ops.NT = 64*1024+ ops.ntbuff; % must be multiple of 32 + ntbuff. This is the batch size (try decreasing if out of memory).
ops.whiteningRange = 32; % number of channels to use for whitening each channel
ops.nSkipCov = 25; % compute whitening matrix from every N-th batch
ops.scaleproc = 200; % int16 scaling of whitened data
ops.nPCs = 3; % how many PCs to project the spikes into
ops.useRAM = 0; % not yet available
%%