Code repository for all of the methods described and the analysis performed in Switching Functional Network Models of Oscillatory Brain Dynamics,The 56th Asilomar Conference on Signals, Systems & Computer, IEEE, 2022.
The detailed explantion of the code is at the top of each Matlab script file. Below is a brief description of all scripts in this repo.
-
core_functions
- skf.m -- switching Kalman filter
- smoother.m -- switching RTS smoother
- em_B.m -- EM on B matrices
-
single_rhythm_model
- single_rhythm_model.m -- propofol example of implementing the common oscillator model with a single rhythm (alpha)
-
multiple_rhythms_model
- multi_rhythms_model.m -- propofol example of implementing the common oscillator model with multiple rhythms (alpha + slow wave)
- std_kf.m -- standard Kalman filter (used for EM on B during awake & unconscious periods)
- std_smth.m -- standard RTS smoother (used for EM on B during awake & unconscious periods)
- kf_em_B.m -- EM on B matrices without switching components
-
plotting
- plt_B_single.m -- plot the estimated B matrices for the single rhythm model
- plt_B_multi.m -- plot the estimated B matrices for the multiple rhythms model
- plt_dosage_bhvr_SW.m -- plot the propofol dosage, behavioral responses, estimated switching states
- helper_functions
- plt_funcs.m -- create the scalp layout of the electrodes
- adjust_sw.m -- adjust the estimated switching states for removed segments