A cortical sheet mesoscopic model for investigating focal seizure onset dynamics (Wang et al. 2014)

 Download zip file 
Help downloading and running models
Accession:155565
The model uses realistically coupled, discretised, Wilson-Cowan units to describe the spatio-temporal activity of a cortical sheet. This model has been used the investigate the dynamic onset mechanisms of focal seizures.
Reference:
1 . Wang Y, Goodfellow M, Taylor PN, Baier G (2014) Dynamic mechanisms of neocortical focal seizure onset PLoS Computational Biology 10(8):e1003787 [PubMed]
Model Information (Click on a link to find other models with that property)
Model Type: Neural mass;
Brain Region(s)/Organism:
Cell Type(s):
Channel(s):
Gap Junctions:
Receptor(s):
Gene(s):
Transmitter(s):
Simulation Environment: MATLAB;
Model Concept(s): Oscillations; Spatio-temporal Activity Patterns; Epilepsy; Delay; Brain Rhythms; Bifurcation;
Implementer(s): Wang, Yujiang [yujiang.wang at newcastle.ac.uk];
/
WebPublication
lib
ConnLocGaussian.m *
ConnPatchyRemOverlap.m *
Conns_n150.mat
convolve2.m *
distSheet.m *
distTorus.m *
exindex.m *
FilterEEG.m
Gaussian.m *
GaussianLocConnFunc.m
generatePatchesOverlap.m *
getDelayMatrix.m
getNoise.m
getParam.m *
getParamDelay.m
makeCellCluster.m *
makeCellClusterToroidal.m *
MayColourMap.mat *
meanMacroCol.m *
runSheet.m *
runSheetDelay.m *
runSheetPRamp.m *
Sigm.m *
                            
function PatchCoordinates = generatePatchesOverlap(numPatches,realNumPatch,cl,neighbourFullIndex,xF,yF,n,patchPercent,nOverlap,PatchCoordinates,CellLocFunc)
%subfunction used by ConnPatchyRemOverlap

        for l=nOverlap+1:numPatches
            xpatch=xF(neighbourFullIndex(randi(cl,1)));
            ypatch=yF(neighbourFullIndex(randi(cl,1)));
            [~,CellLocsVector] = CellLocFunc(1.5, patchPercent, n, ypatch, xpatch);%in this script the coordinate system is where x and y are actually in correct (x=horizontal, y = vertical) directions
            PatchCoordinates(:,l)=CellLocsVector;
        end
        
        
                
        alreadydeleted=find(PatchCoordinates(1,:)==-1);

        if numPatches>realNumPatch+length(alreadydeleted)
            numtodelete=numPatches-(realNumPatch+length(alreadydeleted));
            todelete=randperm(numPatches);
            if isempty(alreadydeleted)
                alreadydeleted=0;
            end
            todelete=todelete(~ismember(todelete,alreadydeleted));
            todelete=todelete(1:numtodelete);
            PatchCoordinates(:,todelete)=-1;
        end
        
end

Loading data, please wait...