Meg information based connectivity coding session

From MEG Core
Jump to navigation Jump to search

We will be meeting virtually over Zoom from 1pm-3pm with breakout rooms for questions and collaboration
There will be a quick discussion and some example data on biowulf.

FRITES toolbox

Club MEG Presentation Recording

 Getting Started Scripts
 Example Scripts

Install in new environment (RECOMMENDED!):

 conda create -n frites conda-forge::mne conda-forge::jupyterlab pip -y 
 conda activate frites
 pip install frites
 git clone
 cd CookingFrites
 jupyter lab
 #Once jupyter is open, the notebooks are in the notebooks folder
 #For MEG data to use with this --- Just test data, there are many optimizations that have not been performed
 pip install git+
 cp /vf/users/MEGmodules/modules/ ./  
 #In python, frites_output_root=$(pwd)/frites_dset
 from frites_connectivity_tutorial.load_data import load_dataset
 dt = load_dataset(<<FRITES_OUTPUT_ROOT>>)

Install in current environment (Better to use New environment as listed above!)-

 #Download getting started tutorials
 git clone
 pip install frites

Use on biowulf

 sinteractive --mem=12G --cpus-per-task=8 
 module use --append /data/MEGmodules/modulefiles
 module load mne_spyder  #loads mne/mne-bids/frites/frites_tutorial

SAM Symoblic Transfer Entropy Connectivity


Club MEG SAM Connectivity Presentation Recording


 This process is a trigger-free process that estimates the symbolic transfer entropy (connectivity assessed by delayed propagation of signals) between two ROIs 
 70-185Hz has been tested and is recommended as a filter band
 Use an atlas based input for the bivariate connectivity assessment (AAL atlas, Desikan K., Destreux...)
 SAMwts or Patch_wts to create the beamformer
 Output is a text file
 Must determine an embedding dimension (4 as default) - the higher the embedding the more data needed
 Minimal requirements (>=1Khz sampling rate and >=300s of data)
 The inputs do not need to be brain regions (for example can compare STE between audio track to brain region)

Use on Biowulf

 #Log into biowulf
 sinteractive --mem=6G --cpus-per-task=4  #You may need more mem/cpus
 #Make sure modules are loadable -- this can also be added to ~/.bashrc file
 module use --append /data/MEGmodules/modulefiles
 module load SAMsrcDev
 #Make the SAM beamformers and apply the symbolic transfer entropy
 sam_cov ...
 sam_wts ...
 STEdelay ....  

 Usage: STEdelay [options]	Version 5.0 (64-bit) rev-1, Aug 10 2022
 If a parameter begins with '--', it is allowed on the command line,
 otherwise it is only allowed in a parameter file (see -m).
 All times are in seconds.
   -h, --help               show this help
   -r DSNAME, --DataSet DSNAME
                            MEG dataset name [env_var = ds]
   -m PFILE, --param PFILE  parameter file name (optionally ending
                            in ".param") [env_var = param]
   -v, --verbose            verbose output
   --CovType GLOBAL|SUM|ALL which covariance matrix to use for the analysis
   --CovBand LO HI          covariance bandwidth limits in Hz
   --ImageBand LO HI        imaging covariance bandwidth limits in Hz
   --OrientBand LO HI       band to use for orientation instead of Global (Hz)
   --SmoothBand LO HI       used to smooth Hilbert envelope, kurtosis, or RVE (Hz)
   --Notch                  apply a powerline notch
                            (including harmonics)
   --Hz HZ                  frequency for powerline (Hz, default 60)
   ImageMetric METRICSPEC   imaging metric
   --Extent DIST            radial extent from ROI voxels or centroid (mm)
                            file containing target coordinates
   -i MRIDIR, --MRIDirectory MRIDIR
                            MRI directory root [env_var = mridir]
   -o IMAGEDIR, --ImageDirectory IMAGEDIR
                            image output directory [env_var = imagedir]
   --PrefixLength N         DataSet prefix may be specified as a number
                            or as a prefix delimiter (default "_")