post_merge_()

After the initial clustering using DPCLUS, units with similar spike waveforms (averaged raw) are merged if their correlation score exceeds P.maxWavCor. Since the probe drift changes the spike waveforms, averaged waveforms are computed at three depth ranges per unit. The spike depth is estimated by the center of mass (CoM) using the site positions and the peak spike amplitude at each site. Since each unit has three average waveforms, maximum correlation score is taken from the 3x3 correlation matrix for each unit pair. Time shift is applied between a pair of unit (P.spkRefrac_ms = 0.25 ms) to find the maximum cross-correlation value.

Syntax

• [S_clu, S0] = post_merge_(S_clu, P, fPostCluster)

Input

• S_clu
• P
• fPostCluster: flag to call postCluster_() function. Defauult: 1. This function assigns the cluster membership

Output

• S_clu: Cluster index S_clu.viClu gets updated, as well as cluster quality metrics, mean unit waveforms and waveform correlation matrix (S_clu.mrWavCor).
• S0

Operations

• postCluster_(): Assign spike memberships based on Rho and Delta values (DPCLUS)
• post_merge_wav_(): Apply automated merging based on the waveform similarity
• S_clu_position_(): Compute cluster positions by computing the center of mass
• S_clu_quality_(): Compute the cluster quality scores