Oops! I didn't see the tilde. Also, I was using the image written out
to the same directory as SPM.mat - the ResMS.img file. If I use
SPM.VResMS, I get the same answer as SPM graph. So then, what is
ResMS.img really?
First, SMP.xX.Bcov isn't really Bcov, then ResMS.img really isn't
ResMS!!! This nomenclature is extremely confusing!!
Allison Nugent
MRI Physicist
SNMAD/MIB/NIMH/NIH
Office: (301)451-8863
Mobile: (301)408-8560
[log in to unmask]
-----Original Message-----
From: Ged Ridgway [mailto:[log in to unmask]]
Sent: Saturday, January 20, 2007 7:25 AM
To: Nugent, Allison C. (NIH/NIMH) [E]
Cc: [log in to unmask]; Wang, Hongye (NIH/NIMH) [C]; Drevets, Wayne
(NIH/NIMH) [E]
Subject: Re: [SPM] Beta SE
Nugent, Allison C. (NIH/NIMH) [E] wrote:
> Hello,
> SO, the problem wasn't ResMS. Examining the code for
> spm_graph.m, around one 130, you can see that the Bcov reported to the
> workspace is ResMS*SPM.xX.Bcov *only* when xSPM.STAT is a 'P'
Other way round actually:
line 142: if xSPM.STAT ~= 'P'
translates as "if xSPM.STAT does not equal 'P'"
So line 149: Bcov = ResMS*SPM.xX.Bcov;
does get run for both T and F statistic images.
> covariance calculated by SPM.xX.Bcov*ResMS is an order of magnitude
> greater than the workspace Bcov following the plot function.
Are you sure you're using the same ResMS as spm_graph, i.e.
line 148: ResMS = spm_get_data(SPM.VResMS,XYZ);
and at the same voxel? (maybe put a debugging break-point after lines
95-97 to check this)
If so, then I can't see what could possibly be different, Bcov is not
modified in spm_graph after line 149, so I expect spm_graph might be
looking at a different voxel to you, maybe?
Hope this helps,
Ged.
|