Dear Chou-Ching,
On Thu, Apr 15, 2010 at 2:50 PM, Chou-Ching K. Lin <[log in to unmask]> wrote:
> Hi, Dear SPM friends, now I am probing and running DCM for M/EEG (induced response). I found some small problems.
> First, it stops at 64 rounds of calculation every time before convergence. Is there any way to change this magic number, 64?
64 iterations should be sufficient to get close to the optimum even
with no convergence but if you'd like to increase this number you can
only do it by changing the code (spm_nlsi_N around line 278).
Also we recently found a bug in how the data for DCM-IR was prepared.
Put the attached function in toolbox\dcm_meeg to fix it.
> Second, when I reload the saved DCM file, all the parameters are recovered, but those dots for B in the lower box are missing (every time).
This might be a GUI problem. I haven't seen it before. You can send me
one example DCM file and I'll check.
> Third, is there any way to examine the calculated parameters (A and B)? When I open the DCM result file in Matlab command line, it shows it is a matrix or structure. But it does not give me the values.
The values are in the structure. I suggest that you put a breakpoint
in spm_dcm_ind_results when you do the results display, find the
relevant part if the code and see how what you are interested in is
extracted from the structure. Some things are not actually stored but
are computed on the fly from things that are stored. If you are just
interested in coupling matrices in Hz you can export them as images
using the option in the 'Results' part.
Best,
Vladimir
|