Demo data with known answers

An 8-channel LabChart-style export, 300 s at 1000 Hz. Channel 6 is the stimulus: 5 V pulses of 5 s every 30 s from t = 30 s (9 pulses). Channel 8 is the LDF signal:

LDF(t) = 120 + 5·sin(2πt/400) + 3·sin(2π·0.13·t) + 1.5·sin(2π·6·t) + noise(SD 2)
         + Σ_onsets 30·g((t − onset)/4),   g(x) = x³·e^(3(1−x)) for x > 0, else 0

That is a ~120 PU baseline, a slow drift (period 400 s), vasomotion at 0.13 Hz, a 6 Hz cardiac ripple and white noise, plus a response of +30 PU peaking 4 s after every onset. The response shape is a gamma-variate function, a common model of slow haemodynamic responses^[12].

Cropped file

20–280 s of the same recording (stim, LDF, t, Fs); onsets at 10, 40, … 250 s.

Trials file

The 8 onsets whose −5 to +20 s window fits, block-averaged to 10 Hz (segmentedLDF, segmentedTime, Fs).

Ground truth

onsets, stimulus duration 5 s, response peak delay 4 s, amplitude 30 PU, baseline 120 PU.

30 s, 8 channels 100 µm apart, whisker stimulus pulses of 20 ms every 2 s from 1 s (15 stimuli). The raw channels are at the TDT rate 24414.0625 Hz; the LFP and stimulus are at 1017.25 Hz (24× decimation).

LFP

ERP(x) = −120 µV·exp(−(x − 15 ms)²/(2·(5 ms)²)) + 60 µV·exp(−(x − 40 ms)²/(2·(12 ms)²))
LFP_c(t) = profile_c·Σ ERP(t − onset) + shared(t) + local_c(t) + white noise (1 µV)
profile_c = exp(−(depth_c − 300 µm)²/(2·(150 µm)²))   → largest on channel 4

The background shared(t) is the same on every channel (an AR(1) process with coefficient 0.98, scale 3 µV), so it cancels in the CSD; local_c(t) is small independent AR(1) noise (0.9, 0.3 µV). Ground truth: onsets, sink channel 4, N1 at 15 ms, P2 at 40 ms, spacing 100 µm.

MUA

Spike trains are Poisson with a 2 ms refractory period; each unit is seen on neighbouring channels with a gain of exp(−|distance|/0.8) (skipped below 0.05). Noise is white, 10 µV. The MUA file holds channels 3–5. Ground truth: every unit's spike times, rates and home channel.

Oscillation demo (core/demo/demoLFPOscillations.m)

The same LFP plus a 6 Hz theta rhythm (40 µV, identical on all channels, frequency wandering with SD 0.5 Hz, so not phase-locked to the stimuli) and a 40 Hz burst (10 µV, 50–250 ms after each stimulus with 25 ms raised-cosine ramps, channels 3–5, same phase at every stimulus).

Basic stack (DemoData.imagingStack)

96 × 96 px, 150 frames at 10 Hz. A dark vertical vessel at x = 60 whose diameter is 12 + 3·sin(2π·0.2·t) px (9–15 px, period 5 s); a bright red blood cell (Gaussian spot, sigma 2 px) moving down the upper part of the vessel at 2 px per frame (it never reaches the diameter line at y = 70); a cell at (24, 30), radius 6 px, with calcium transients of ΔF/F = 1 at 3, 7 and 11 s (rise 0.1 s, decay 0.8 s) added on top of the tissue background; sensor noise SD 0.02. roiMask marks the cell.

Advanced stack (core/demo/demoImagingAdvanced.m)

• Every frame translated by a sub-pixel shift that follows a smooth random walk, zero mean, at most 3 px in each direction (the true shifts are stored).
• Three cells: (22, 24) r = 6 px, events at 4.0, 8.5, 13.0 s; (26, 78) r = 5 px, events at 5.5, 10.5 s; (82, 30) r = 7 px, events at 7.0, 12.0 s (ΔF/F amplitudes 1.0, 0.8, 1.2).
• The vessel at x = 60 (12 ± 3 px, period 5 s) and a red blood cell moving 3 px per frame down the whole image, crossing the diameter line (35, 64)–(85, 64) about every 3 s.

• Scalp: an oddball study of 8 participants, 32 channels at 250 Hz, trials from −200 to 800 ms, already cleaned; Standard (40), Target (15) and Novel (15) trials, 5 rejected per participant. P300 at 350 ms, largest at Pz: Target 10 > Novel 6 > Standard 2 µV; P1 +2 µV at 60 ms (Oz), N1 −5 µV at 100 ms (Cz), 10 Hz alpha in random phase.
• Rodent: a continuous 60 s recording from 4 skull screws at 1000 Hz with 30 light flashes; a visual evoked potential over V1 (−40 µV at 50 ms, +25 µV at 100 ms), 30% of that over M1.
• Both are written in every EEG format the window reads (EEGLAB, FieldTrip, BrainVision, plain .mat). Expected results are on the EEG page.

24 LDF trial files: the same 8 animals in three conditions (Control, Stimulated, Drug), 8 trials each (−5 to 20 s at 10 Hz), response peaking 4 s after onset on a 110–130 PU baseline.

amplitude(animal, condition) = trueMean(condition) + animal effect (SD 3 PU)
                              + animal × condition scatter (SD 2.5 PU)
trueMean = 18, 30, 24 PU (Control, Stimulated, Drug);  each trial adds jitter (SD 2 PU)

Population effects: Stimulated − Control = 12 PU, d_z = 3.27 (paired), Cohen's d = 3.02 (unpaired), ANOVA η² = 0.60. The realised per-file amplitudes are stored as well.