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