Dear Meng:
I think you need to do this via the matlab command line. The code below
outlines the steps, and is taken from spm_dcm_ui.
% select and load a DCM file.
load(spm_get(1,'DCM*','Select a DCM file'));
% number of inputs
m = DCM.M.m;
% number of outputs
l = DCM.M.l;
T = XX; % enter whatever threshold you want
% calculate the posterior probabilities
pp = 1 - spm_Ncdf(T,abs(DCM.Ep),diag(DCM.Cp));
% reshape into 3 matrices.
[pA pB pC] = spm_dcm_reshape(pp,m,l,1);
% pA contains the probabilities for matrix A, pB for matrix B, etc. at the
selected threshold.
Darren
At 10:16 PM 4/8/2005, Meng Liang wrote:
>Dear SPMers,
>
>A previous email titled "Re: [SPM] DCM output conflicts" has mentioned
>that "DCM is calculated the probabilities are stored using a threshold of
>log(2)/4, while with displaying DCMs you may have used a threshold of 0 or
>other threshold."
>But if i want to get the probabilities of connectivities larger than a
>threshold of 0, what should I do?
>To my knowledge, one can not get the values of the probabilities using the
>graphics interface of SPM2. (Maybe I'm wrong?)
>
>Any suggestions are appreciated very much!
>
>
>¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡Meng Liang
>¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡[log in to unmask]
>¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡¡2005-04-09
-------------------------------------------------------------------------
Darren R. Gitelman, M.D.
Cognitive Neurology and Alzheimer¹s Disease Center
Northwestern Univ., 320 E. Superior St., Searle 11-470, Chicago, IL 60611
Voice: (312) 908-9023 Fax: (312) 908-8789
-------------------------------------------------------------------------
|