load('mri.mat')
load('segmentedmri.mat')
load('headmodel_struct.mat')
load('raw_filtered.mat')
load('clean_segmentedData.mat', 'clean_segmentedData')
load('erp.mat')
clean_segmentedData.grad = raw_data.grad;
dataAll = clean_segmentedData;
cfg = [];
cfg.covariance = 'yes';
cfg.covariancewindow = 'all';
cfg.vartrllength = 2;
timelockAll = ft_timelockanalysis(cfg, dataAll);
erp = rmfield(erp,'grad');
cfg=[];
cfg.layout = '4D248.lay';
cfg.interactive = 'yes';
figure;ft_singleplotER(cfg,erp)
% 1st time - symmetric
cfg.latency = [0.115 0.135];
cfg.numdipoles = 2;
cfg.symmetry = 'y';
cfg.grid.resolution = 1;
cfg.grid.unit = 'cm';
cfg.gridsearch = 'yes';
cfg.vol = headmodel_struct;
cfg.channel = 'MEG';
source_planar = ft_dipolefitting(cfg, timelockAll);
% 2nd time - nonsymmetric
load('mri.mat')
load('segmentedmri.mat')
load('headmodel_struct.mat')
load('raw_filtered.mat')
load('clean_segmentedData.mat', 'clean_segmentedData')
load('erp.mat')
clean_segmentedData.grad = raw_data.grad;
dataAll = clean_segmentedData;
cfg = [];
cfg.covariance = 'yes';
cfg.covariancewindow = 'all';
cfg.vartrllength = 2;
timelockAll = ft_timelockanalysis(cfg, dataAll);
erp = rmfield(erp,'grad');
cfg=[];
cfg.layout = '4D248.lay';
cfg.interactive = 'yes';
figure;ft_singleplotER(cfg,erp)
% 1st time - symmetric
cfg = [];
cfg.latency = [0.115 0.135];
cfg.numdipoles = 2;
cfg.symmetry = 'y';
cfg.grid.resolution = 1;
cfg.grid.unit = 'cm';
cfg.gridsearch = 'yes';
cfg.vol = headmodel_struct;
cfg.channel = 'MEG';
source_planar = ft_dipolefitting(cfg, timelockAll);
% 2nd time - nonsymmetric
cfg = [];
cfg.latency = [0.115 0.135];