GetSpikeFeatures

GetSpikeFeatures - Get spike features.

  USAGE

    features = GetSpikeFeatures(group)

    group          electrode group

  OUTPUT

    features       (timestamp,group,cluster,features...) t-uples.

  NOTE

    Contrary to e.g. GetSpikes, this function retrieves information for all clusters
    of a given electrode group, rather than from a list of selected units. This is
    because .fet files are usually very large and it would be very inefficient to
    repeatedly read the same file for each cluster. It is preferable to load the whole
    file only once, and then select individual clusters if appropriate.

  SEE

    See also LoadSpikeFeatures, IsolationQuality.

0001 function [data,indices] = GetSpikeFeatures(group)
0002 
0003 %GetSpikeFeatures - Get spike features.
0004 %
0005 %  USAGE
0006 %
0007 %    features = GetSpikeFeatures(group)
0008 %
0009 %    group          electrode group
0010 %
0011 %  OUTPUT
0012 %
0013 %    features       (timestamp,group,cluster,features...) t-uples.
0014 %
0015 %  NOTE
0016 %
0017 %    Contrary to e.g. GetSpikes, this function retrieves information for all clusters
0018 %    of a given electrode group, rather than from a list of selected units. This is
0019 %    because .fet files are usually very large and it would be very inefficient to
0020 %    repeatedly read the same file for each cluster. It is preferable to load the whole
0021 %    file only once, and then select individual clusters if appropriate.
0022 %
0023 %  SEE
0024 %
0025 %    See also LoadSpikeFeatures, IsolationQuality.
0026 
0027 % Copyright (C) 2004-2011 by Michaël Zugaro
0028 %
0029 % This program is free software; you can redistribute it and/or modify
0030 % it under the terms of the GNU General Public License as published by
0031 % the Free Software Foundation; either version 3 of the License, or
0032 % (at your option) any later version.
0033 
0034 global DATA;
0035 if isempty(DATA),
0036     error('No session defined (did you forget to call SetCurrentSession? Type ''help <a href="matlab:help Data">Data</a>'' for details).');
0037 end
0038 
0039 if nargin < 1,
0040   error('Incorrect number of parameters (type ''help <a href="matlab:help GetSpikeFeatures">GetSpikeFeatures</a>'' for details).');
0041 end
0042 
0043 if ~isiscalar(group,'>=0'),
0044     error('Incorrect unit (type ''help <a href="matlab:help GetSpikeFeatures">GetSpikeFeatures</a>'' for details).');
0045 end
0046 
0047 % Load all waveforms for this group
0048 filename = [DATA.session.path '/' DATA.session.basename '.fet.' int2str(group)];
0049 rate = DATA.rates.wideband;
0050 data = LoadSpikeFeatures(filename,rate);