Scripts for statistical processing of patient and simulated datas for the article Transition to seizure in focal epilepsy: from SEEG phenomenology to underlying mechanisms" for publication in Epilepsia

Execute in order of cells.

For all questions, contact me on github, or by mail at mehmet-alihan dot kayabas at inserm dot fr or the corresponding author "Fabrice Wendling"

Inputs for patients:

• for each patient, three lists of h² values for each comparisons of the three selected channels (i.e. channel 1 vs 2, channel 1 vs 3 and channel 2 vs 3) => 'patients/','*.dat' files
• timestamps of start and end of preictal, fast-onset and ictal periods => 'patients/timestamps_patient.csv'
• power values at each subbands (theta, alpha, beta and gamma) => 'patients/psd_subbands_data_out.csv'
• Names of the three SEEG electrode channels selected => 'patients/montage_data/' Simulations are performed on Coalia software. h² values and power values are calculated on Amadeus software.

Inputs for simulations:

• for each of the ten simulations, two lists of h² values (channel 1 vs 2 and channel 2 vs 1). In the script we take the average of the $h^2_{1,2}$ and $h^2_{2,1}$ as the focus of the study was not directionality. => 'simulations/h2_simus/'
• power values at each subbands (thet,a alpha, beta and gamma) => 'simulations/power_simu_data.csv'
• firing rates of the pyramidal cells and PV interneurons => 'simulations/PSPout_EXC.txt'
  • (for figure 5) corresponding simulated timeseries for the firing rates => 'simulations/phase_value.dat'

NB1: In the scripts, FOA (fast-onset activity) period is interchangeably written as LVF (low-voltage fast).

NB2: Some simulation results might vary due to stochasticity (external input, the seed of RNG). The general h² synchronization pattern remains valid: Correlation during preictal phase, followed by a decorrelation at FOA followed by a correlation during the ictal phase.