That's a good point. You should start your script with:
spm('defaults', 'eeg');
Sorry I forgot to include that.
Vladimir
On Mon, Apr 19, 2010 at 10:18 PM, Jeremy Nelson <[log in to unmask]> wrote:
> Hi Vladimir,
>
>> OK. That's much clearer now. What you can do is press '3D source
>> reconstruction', load the example file, delete the source
>> reconstruction data that is there and save. Then when you load it to
>> DCM you'll be asked at some point to specify the head model parameters
>> from scratch. That should solve this problem.
>
> This worked - thank you for the help!
>
>> You should not put all the ERPs from all subjects in the same dataset.
>> You can either compute the grand average and do DCM of that or do DCM
>> of each subject individually and then do statistics on the parameters
>> across subjects. Look at the published DCM paper of Marta Garrido et
>> al. for some examples. In any case what you should input to DCM are
>> files with 3 trials, one per condition. You don't need to specify the
>> number of repetitions that went into each average for DCM. You should
>> specify the condition labels by changing the line:
>>
>> D = conditions(D, 1:Ntrials, 'Condition 1'); % Sets the condition label
>>
>> from the script and changing into something like:
>>
>> D = conditions(D, 1, 'Pre');
>> D = conditions(D, 2, 'Post 1');
>> D = conditions(D, 3, 'Post 2');
>
> Yes, that is the approach I am using - each subject's pre, post 1, and
> post 2 data are being compiled into the same dataset (one dataset per
> subject, each of which will be run through DCM). After making the
> changes you suggested to the script (which I have attached), here are
> the errors I am getting when I try to run it:
>
> checkmeeg: no channel type, assigning default
> checkmeeg: no units, assigning default
> checkmeeg: transform type missing, assigning default
> checkmeeg: data scale missing, assigning default
> checkmeeg: data type is missing or incorrect, assigning default
> ??? Undefined function or method 'fileio_chantype' for input arguments of
> type 'cell'.
>
> Error in ==> spm_eeg_prep at 60
> type = fileio_chantype(D.chanlabels);
>
> Error in ==> spm_eeg_ft2spm at 145
> D = spm_eeg_prep(S1);
>
> Error in ==> spm_eeg_convert_arbitrary_data at 127
> D = spm_eeg_ft2spm(ftdata, fname);
>
> Error in ==> run at 57
> evalin('caller', [s ';']);
>
>
> Do you have any suggestions?
>
> Thanks again,
>
> Jeremy
>
|