Question: Read_ctf_rc usw..

Matsuhashi, Masao (NIH/NINDS) matsuham at NINDS.NIH.GOV
Fri Oct 8 03:46:19 CEST 2004


Dear Jan-Mathijs,

Thank you very much for the description. Now I am trying...

Actually the reference signal is not included in the dataset, so I did as
follows:

========================
% load data, in 5 sec segments
cfg.dataset='X:\raid\DOSTPIDM_track_20040826_01.ds'                     % 6+
min, 600Hz data, 219600 time points
cfg.trl=[[1:3000:219600-3000]',[3000:3000:219600]',zeros(fix(219600/3000),1)
]
cfg.channel='MEG'
cfg.detrend='yes'
data=preprocessing(cfg)       %take some time, so do it only neccesory.
% combine with ref signal  w{}
for loop = 1:length(data.trial)
  data.trial{loop}=[data.trial{loop}; w{1}([1:3000]+3000*(loop-1),1)'];
%w is a cell array of ref signals (219600 x 2 matrix)
end
data.label=[data.label;{'REF'}];
% freq analysis
cfg=[];
cfg.method = 'fft';
cfg.sgncmb = channelcombination({'MEG' 'REF'},data.label);
cfg.output     = 'powandcsd'
cfg.foilim = [0,20]
cfg.tapsmofrq=0.5;
cfg.pad='maxperlen';
[freq] = freqanalysis(cfg, data)
freq=freqdescriptives(cfg, freq)
==========================

And I encountered an error:

==========================
....
processing trial 71, nfft: 3000 samples, taper length: 3000 samples, 4
tapers
processing trial 72, nfft: 3000 samples, taper length: 3000 samples, 4
tapers
processing trial 73, nfft: 3000 samples, taper length: 3000 samples, 4
tapers

freq =

        label: {275x1 cell}
          sgn: {275x1 cell}
       dimord: 'sgncmb_frq'
    powspctrm: [275x101 double]
       sgncmb: {274x2 cell}
     labelcmb: {274x2 cell}
    crsspctrm: [274x101 double]
          foi: [1x101 double]
         grad: [1x1 struct]
          cfg: [1x1 struct]

??? Reference to non-existent field 'pow'.

Error in ==> C:\MATLAB\fieldtrip\freqdescriptives.m
On line 207  ==> frqjck.powspctrm = frqjck.pow;

>>
===========================

I think I am doing something wrong......

------------------------------------------------
Masao Matsuhashi, MD, DMS
NINDS, NIH
Tel: 1-301-594-9139
E-mail: MatsuhaM at ninds.nih.gov
------------------------------------------------



More information about the fieldtrip mailing list